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A  study  is  presented  of  some  minimax  location  problems  in  the 
Euclidean  plane  under  the  assumption  that  distance  between  points  is 
the  rectilinear  distance.   Each  problem  is  to  determine  the  location 
of  new  facilities  with  respect  to  existing  facilities  so  as  to  min- 
imize the  maximum  cost,  where  costs  are  specified  functions  of  the 
rectilinear  distance  between  new  and  existing  facilities  or  among  new 
facilities.   For  all  but  one  of  the  problems  considered,  upper  bound 
constraints  are  assumed  on  the  rectilinear  distances  between  new  and 
existing  facilities. 

A  decomposition  into  independent  one-dimensional  minimax 
subproblems  is  obtained  for  all  but  one  of  the  problems.   The  optimal 
solutions  to  the  subproblems  determine  the  optimal  solutions  to  the 
original  problem;  solution  procedures  are  developed  for  the  subprob- 
lems, based  on  their  structure.   The  structure  of  some  of  the  sub- 
problems  gives  an  explicit  expression  for  their  unique  solution  in 
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terms  of  parameters  of  the  subproblem.   The  structure  of  other  subprob- 
lems  enables  an  equivalent  network  flow  problem  to  be  formulated  which 
is  efficiently  solved  by  determining  a  simple  chain  of  maximum  cost 
to  weight  ratio  in  the  network. 

For  the  remaining  problem  considered,  properties  of  the  recti- 
linear distance  enable  the  specification  of  a  finite  iterative  proce- 
dure, where  each  iteration  solves  a  set  covering  problem,  to  determine 
an  optimal  location  of  facilities. 


IX 


CHAPTER  1 
INTRODUCTION  AND  OVERVIEW 

1. 1  Introduction 

This  dissertation  contains  a  study  of  some  minimax  location 
problems  in  the  Euclidean  plane  under  the  assumption  that  distance 
between  points  is  the  rectilinear  distance.   These  problems  may  be 
characterized  as  determining  the  locations  of  new  facilities  with 
respect  to  existing  facilities  so  as  to  minimize  the  maximum  cost, 
where  costs  are  specified  functions  of  the  rectilinear  distance  between 
new  and  existing  facilities  or  among  new  facilities.   An  example  of 
such  a  problem  is  to  determine  the  location  of  one  point  in  the  plane 
that  minimizes  the  maximum  rectilinear  distance  from  it  to  any  of  a 
finite  number  of  fixed  points  in  the  plane. 

1.2  Applications 

Minimax  location  problems  considered  in  this  dissertation  may 
have  applications  in  determining  the  location  of  emergency  service  facil- 
ities such  as  police  stations,  fire  stations  or  ambulance  centers  that 
must  respond  to  emergencies  at  fixed  locations.   For  example,  it  might 
be  desired  to  minimize  the  maximum  distance  traveled  from  an  ambulance 
center  to  an  emergency  at  any  fixed  location.   Similar  applications 
might  be  found  in  locating  emergency  repair  centers  that  must  respond 


quickly  to  breakdowns  at  fixed  facilities.   Applications  also  may  be 
found  in  locating  facilities  that  emphasize  accessibility  or  conven- 
ience to  customers,  such  as  public  libraries,  self-service  post  offices 
or  banks. 

The  rectilinear  distance  assumption  is  made  to  approximate  the 
distance  between  two  points  when  travel  is  restricted  to  a  system  of 
orthogonal  paths.   Applications  using  the  rectilinear  distance  may 
arise  in  locating  facilities  within  a  plant  or  warehouse  where  travel 
occurs  on  a  system  of  orthogonal  aisles.   Applications  may  also  be 
found  in  locating  facilities  in  an  urban  area  where  travel  is  restricted 
to  a  system  of  orthogonal  streets. 

There  are  recognized  limitations  in  applying  the  problems 
studied  here  to  locate  facilities  in  many  actual  settings.   All  the 
factors  that  require  consideration  for  the  location  of  some  facility 
might  not  be  included  in  the  assumptions  of  the  problems  studied. 
However,  good  methodology  and  solution  procedures  for  these  problems 
form  a  beginning  for  the  study  of  problems  containing  more  realistic 
considerations.   Also,  solutions  to  some  of  the  problems  in  the  disser- 
tation may  be  used  as  guidelines  for  locating  a  facility  where  there 
are  additional  factors  to  consider. 

1.3  Related  Work 

Most  of  the  problems  considered  in  this  dissertation  are  anal- 
ogues to  "total  cost"  location  problems.  The  objective  of  these  prob- 
lems is  to  minimize  the  sum  of  costs,  where  costs  are  functions  of  the 
distance  between  the  new  and  existing  facilities.   Perhaps  the  earliest 


total  cost  location  problem  was  posed  in  the  seventeenth  century  by 
Fermat  using  Euclidean  distance.   Since  then,  variations  and  general- 
izations of  the  problem  have  received  considerable  attention  in  the 
literature.   Recent  contributions  to  the  total  cost  problem  using 
Euclidean  distances  are  reported  by  Kuhn  [20,21],  Kuhn  and  Kuenne  [22], 
Love  [23],  and  Francis  and  Cabot  [14],  and  for  general  distances  by 
Wendell  [28].   These  papers  also  contain  extensive  references  to  loca- 
tion problems.   Two  versions  of  the  total  cost  location  problem  assum- 
ing rectilinear  distances  are  considered  by  Francis  [10,11],  Wesolowsky 
and  Love  [30],  and  Cabot,  Francis  and  Stary  [2].   The  last  paper  also 
contains  an  extensive  list  of  references  for  location  problems. 

A  minimax  location  problem  using  Euclidean  distance  was  posed 
by  Sylvester  [26]  in  1857.   Recent  contributions  to  two  versions  of 
this  problem  are  reported  by  Elzinga  and  Hearn  [6,7].   Properties  of 
a  minimax  location  problem  with  general  distances  are  developed  by 
Francis  T13]. 

Literature  concerned  with  minimax  location  problems  explicitly 
using  the  rectilinear  distance  has  appeared  only  since  1971.   Francis 
[15]  and  Elzinga  and  Hearn  [7]  have  solved  a  special  case  of  one  of 
the  problems  considered  in  this  dissertation.   In  [16] ,  Francis  formu- 
lates special  cases  of  other  problems  considered  here  as  linear 
programs.   Wesolowsky  [29]  also  considers  some  special  cases  of  problems 
considered  here  and  suggests  that  they  may  be  solved  by  parametric  lin- 
ear programming.   This  literature  will  be  considered  in  more  detail  as 
it  applies  to  the  specific  minimax  problem  developed  in  each  of  the 
following  chapters. 


1.4   Overview 

Chapters  2,  3,  4  and  5  each  contains  a  study  of  a  different 
minimax  location  problem.   Each  chapter  includes  a  formulation  of  the 
problem  under  consideration,  a  discussion  of  some  possible  applications 
and  a  review  of  the  literature  associated  with  the  problem.   Then  prop- 
erties of  the  problem  are  developed  and  a  solution  procedure  is  given. 

An  important  aspect  of  the  presentation  in  this  dissertation 
is  a  transformation  of  the  problems  under  consideration  into  equiva- 
lent minimax  location  problems  in  another  space.   This  transformation 
was  first  applied  to  a  minimax  location  problem  by  Francis  ("15].   Here 
the  transformation  is  shown  to  apply  to  the  more  general  problems  that 
are  considered.   For  each  problem  discussed  in  Chapters  2,  3  and  5,  the 
transformed  problem  decomposes  into  independent  subproblems  which  are 
minimax  location  problems  in  one  dimension.   The  solutions  to  the  sub- 
problems  provide  the  solutions  to  the  original  problems.   For  the  prob- 
lem studied  in  Chapter  4,  the  transformed  problem  does  not  decompose 
but  the  transformation  is  used  in  the  solution  procedure.   Chapter  2 
includes  a  development  of  this  transformation. 

There  is  also  discussed  in  Chapter  2  the  one  facility  minimax 
location  problem.   This  problem  is  to  locate  a  single  facility  in  the 
plane  with  respect  to  a  finite  number  of  existing  facilities.   Upper 
bound  constraints  on  the  rectilinear  distance  between  the  new  facility 
and  the  existing  facilities  are  assumed.   Properties  and  conditions  at 
optimality  are  developed  and  a  solution  procedure  based  on  these  condi- 
tions is  presented.   A  lower  bound  and  a  procedure  for  determining 
level  sets  for  different  values  of  the  objective  functions  are  also 
presented. 


Considered  in  Chapter  3  is  a  multif acility  minimax  location 
problem  for  locating  n  new  facilities  with  respect  to  existing  facil- 
ities.  The  costs  for  this  problem  are  proportional  to  the  rectilinear 
distances  among  facilities,  that  is,  between  new  and  existing  facil- 
ities and  among  new  facilities.   Upper  bound  constraints  on  the  recti- 
linear distance  among  facilities  are  also  assumed.   The  subproblems 
obtained  from  the  decomposition  of  the  transformed  problem  are  formu- 
lated as  equivalent  network  weighted  flow  problems  for  which  a  solution 
procedure  is  given. 

The  problem  considered  in  Chapter  4  is  a  multif acility  minimax 
location  problem  referred  to  as  the  p-center  problem.   This  problem  is 
to  locate  p  new  facilities  with  respect  to  existing  facilities  so  as 
to  minimize  the  maximum  cost.   In  the  p-center  problem,  there  is  a  cost 
for  each  existing  facility  which  is  the  minimum,  over  the  new  facil- 
ities, of  a  linear  function  of  the  rectilinear  distance  between  the 
existing  facility  and  each  new  facility.   From  the  transformed  problem 
there  is  determined  a  sequence  of  set  covering  problems.   The  solutions 
to  these  set  covering  problems  determine  the  optimal  p-center  for  all 
possible  values  of  p. 

Discussed  in  Chapter  5  is  a  one  facility  minimax  location 
problem  for  locating  a  single  facility  with  respect  to  existing  facil- 
ities, but  where  the  existing  facilities  are  points  of  a  compact 
subset  of  the  plane.   Properties  developed  by  Francis  [13]  for  minimax 
location  problems  with  respect  to  a  finite  set  of  existing  facilities 
are  extended  to  a  compact  set  of  existing  facilities.   These  properties 


together  with  the  results  of  Chapter  2  are  used  to  develop  properties 
and  solution  procedures  for  the  problem. 

Chapter  6  contains  a  summary  of  the  dissertation  and  suggests 
some  areas  for  further  research. 


CHAPTER  2 
A  ONE  FACILITY  PROBLEM 

2. 1  Introduction 

The  problem  considered  in  this  chapter  is  to  locate  one  new 
facility  with  respect  to  a  finite  number  of  existing  facilities  so  as 
to  minimize  the  maximum  cost,  where  costs  are  linear  functions  of  the 
rectilinear  distance  between  the  new  and  existing  facilities.   The 
problem  also  includes  a  constraint  for  each  existing  facility  which 
is  an  upper  bound  on  the  rectilinear  distance  from  that  facility  to 
the  new  facility. 

A  transformation,  introduced  in  [15],  is  shown  to  apply  to 
this  problem  and  to  give  a  decomposition  into  two  independent  one- 
dimensional  minimax  subproblems.   All  optimal  solutions  to  the  orig- 
inal one  facility  problem  are  obtained  from  the  solutions  to  the  sub- 
problems.    Solution  procedures  are  developed  for  each  subproblem  in 
terms  of  the  maximum  of  a  finite  number  of  simple  functions  of  the 
subproblem' s  parameters. 

2.2  Formulation 

The  one  facility  minimax  location  problem  is  now  formulated 
and  notation  to  be  used  throughout  the  dissertation  is  defined. 

A  notational  convenience  adopted  here  is  that  "l  <  i  <  m" 
denotes  all  the  integers  and  only  the  integers  from  1  to  m. 


Let  X  =  (x,y)  and  P  =  (a,b)  by  any  two  points  in  the  Euclidean 
plane,  E_.  Denote  the  rectilinear  distance  between  X  and  P  as  ||  X  —  Pjf 
so  that 


II  X  -  Pl^  =  ||  (x,y)  -  (a.b))^  =  |x  -  a|  +  |y  -  b|  . 

Let  X  denote  the  point  in  E  to  be  determined  for  the  location 
of  the  new  facility.   Let  S  be  a  finite  set  in  E„  of  m  points  P  , 
1  <  i  <  m,  representing  existing  facilities.   For  each  P.  let  g  be 
any  real  constant  and  let  w.  be  any  real  positive  constant. 

Define  a  function  f(X)  by 

Each  expression,  w  ||x-r|L  +  g.,  for  1  <  i  <  m,  is  a  linear 
function  of  the  rectilinear  distance  between  X  and  P. ,  referred  to 
previously  as  a  cost.   Results  in  [12]  extend  to  show  that  each  of 
these  expressions  is  a  continuous,  convex  function  in  X.   Thus  the 
function  f(X)  is  continuous  and  convex  since  it  is  the  pointwise  max- 
imum of  a  finite  number  of  continuous,  convex  functions  in  X. 

For  each  point  P  let  c.  be  a  positive  constant  for  1  <  i  <  m 
11 

which  is  an  upper  bound  on  the  rectilinear  distance  between  P  and  X. 
With  the  above  terminology,  the  one  facility  minimax  location 
problem  to  be  studied  in  this  chapter  can  be  stated  as 


mm    max 


i{-i»»-'i«i*«J 


X£E0   l<i<m 


s.t.    I  X  -  P. ||.  <  c.    1  <  i  <  m. 


(2.2.1) 


i"  1 


Observe  that  the  unconstrained  one  facility  minimax  location 
problem 

|w.||  X  -  P.||1  +  gj  (2.2.2) 


mm    max 
X6E   l<i<m 


is  a  special  case  of  problem  (2.2.1).   This  follows  by  taking  the  c. 
in  problem  (2.2.1)  sufficiently  large  so  that  the  constraints  are  not 
binding. 

A  generalization  of  problem  (2.2.1)  considered  in  Section  2.8 
allows  the  w.  to  be  nonnegative.   However,  if  all  w.  =  0,  the  objec- 
tive function  value  is  independent  of  the  location  of  X,  and  the  prob- 
lem is  not  a  well-formulated  location  problem.   Therefore,  it  will  be 
assumed  that  some  of  the  w.  are  positive  when  considering  this  general- 
ization. 

2.3   Applications 

Possible  applications  for  the  one  facility  minimax  location 
problem  are  included  in  those  uses  suggested  in  Chapter  1  for  more 
general  minimax  location  problems. 

As  a  specific  example,  [16,19],  suppose  it  is  desired  to  locate 
an  ambulance  center  in  a  city  to  respond  to  m  existing  facilities  in 
the  city.   Each  response  involves  traveling  to  one  of  the  existing 
facilities  and  then  to  a  hospital.   The  ambulance  center  is  to  be 
located  so  that  the  maximum  distance  traveled  to  respond  to  an  emergency 
is  minimized.   In  addition,  there  are  requirements  that  the  ambulance 
center  be  within  a  specified  distance  from  each  existing  facility. 

Let  P. ,  for  1  <  i  <  m,  be  the  location  of  an  existing  facility. 

Take  w.  =  1  and  g.  to  be  the  distance  from  the  facility  at  P   to  the 
li  '  i 


10 


hospital.   Take  c.  to  be  the  upper  bound  on  the  distance  that  the 
ambulance  center  must  be  located  from  the  facility  at  P..   Then  this 
problem  may  be  formulated  as  problem  (2.2.1). 

The  numbers  w  may  have  several  interpretations  as  constants 
per  unit  of  distance  between  X  and  P.  with  corresponding  interpreta- 
tions for  the  constants  g..   For  example,  in  locating  a  service  repair 
center  at  some  point  X  in  a  plant,  it  might  be  desired  to  minimize 

the  maximum  cost  of  travel  from  X  to  the  fixed  locations  P..   In  this 

1 

application  the  w.  may  be  interpreted  as  costs  per  unit  distance  and 

the  g.  as  constant  costs.   In  another  application  of  locating  an 

emergency  center  at  some  point  X,  it  might  be  desired  to  minimize  the 

maximum  time  to  travel  from  the  fixed  locations  P.  to  X.   Then  w. 

l  l 

would  be  interpreted  as  time  per  unit  distance  and  the  g.  as  fixed 
times  associated  with  each  point  P  . 

The  w.  may  also  be  interpreted  as  a  weight  or  a  measure  of 
relative  importance  given  to  each  point  P..   For  example,  in  locating 
a  bank  in  a  city,  customers  located  at  the  points  P.  may  be  given  dif- 
ferent weights  according  to  their  frequency  of  travel  to,  or  level  of 
business  with,  the  bank.   Then  locating  the  bank  so  as  to  minimize  the 
maximum  weighted  distance  that  customers  have  to  travel  may  be  inves- 
tigated using  the  problem  formulated  as  (2.2.1). 

Another  interpretation  of  the  points  P.  is  obtained  by  ideal- 
izing an  area  as  a  point.   Then  the  one  facility  problem  may  find 
applications  in  locating  a  facility  with  respect  to  given  areas  or 
regions  in  an  urban  or  industrial  setting. 


11 


2.4   Related  Work 

Francis  [12]  considered  properties  of  a  one  facility  minimax 
location  problem  for  general  distances  in  1966.   However,  the  one 
facility  problem  using  the  rectilinear  distance  appears  to  have 
received  explicit  consideration  only  since  1971.   Elzinga  and  Hearn 
[7]  consider  a  geometrically  based  procedure  that  gives  an  explicit 
solution  to  a  special  case  of  problem  (2.2.2)  obtained  by  setting  w.  =  1 
and  g.  ^  0  for  1  <  i  <  m.   Francis  [15]  considers  this  problem  with 
g.  =0,  and  introduces  a  transformation  of  the  problem  that  gives  an 
explicit  solution  procedure.   This  transformation  is  developed  in 
Section  2.5  for  application  to  problem  (2.2.1)  and  to  the  other  minimax 
location  problems  considered  in  the  sequel. 

In  [16],  Francis  considers  some  special  cases  of  problem  (2.2.2) 
and  formulates  them  as  equivalent  linear  programs.   For  the  special 
case  with  w.  =  1  and  g.  =  0,  he  obtains  an  explicit  solution  based  on 
the  structure  of  the  linear  program. 

Using  the  development  in  [16],  an  equivalent  linear  program 
to  (2.2.1)  is  obtained.   A  new  variable  z  is  introduced  and  the  con- 
strained problem 
min  z 
s.t.    w.j|  X  -  Pi||1  +  g.  <  z 


X  "  Pi"l      *%, 


1  <  i  <  m  (2.4.  1) 


is  written  which  is  equivalent  to  problem  (2.2.1).   Each  constraint  of 
(2.4.1)  may  be  written  equivalently  as  four  linear  inequality  constraints, 
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giving  a  linear  programming  problem.   The  linear  programming  problem  is 
not  used  directly  in  the  solution  procedure  developed  in  the  sequel; 
however,  this  formulation  is  useful  in  proving  some  results  that  follow 
and  in  providing  insight  into  the  one  facility  problem. 

Wesolowsky  [29]  considers  the  problem  (2.2.2)  with  w.  >  0 
and  g.  =0,  l<i<m.   He  suggests  a  graphical  solution  procedure 
but  notes  "when  m  is  large  this  graphical  trial  and  error  method  is 
unwieldy."   He  also  reports  that  the  problem  can  be  written  as  a  linear 
program  but  erroneously  states  "the  approach  is  only  workable  in  the 
one-dimensional  case." 

Thus,  while  the  literature  reports  that  some  special  cases  of 
problem  (2.2.1)  are  equivalent  to  linear  programs,  an  explicit  solu- 
tion procedure  for  the  problem  has  not  been  developed.   The  notion  of 
the  transformation  does  not  appear  to  have  been  applied  to  minimax 
location  problems  previous  to  [15]  and  has  not  been  developed  as  a  gen- 
eral tool  for  application  to  the  minimax  location  problems  studied  here. 

2.5  A  Transformation  and  Decomposition 

This  section  provides  an  analytical  development  of  a  transfor- 
mation that  is  applied  to  the  one  facility  problem  and  allows  a  decom- 
position into  simpler  subproblems.   A  geometrical  development  is  given 
for  this  transformation  in  [15]  where  it  is  applied  to  a  special  case 
of  problem  (2.2.2).   The  geometrical  development  is  based  on  the  obser- 
vation that  the  locus  of  all  points  in  E  of  equal  rectilinear  distance 
from  a  point  X  in  E„  is  a  "diamond"  with  sides  of  equal  length  and 
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slope  ±  1.   Rotating  this  diamond  by  45  degrees  gives  a  square  with 

sides  parallel  to  the  axes,  and  a  square  is  the  locus  of  all  points 
in  E„  of  equal  Chebyshev  distance  from  a  point  X  in  E  . 

The  Chebyshev  distance  between  two  points  X  and  P  in  E„, 
denoted  by  ||  X  -  p|j   ,  is  defined  by 

II  X  -  v\\a    =  ||  (x,y)  -  (a,b)||B  =  max  (|x-a|,  |y-b|)  . 

The  transformation,  denoted  by  T,  is  a  counterclockwise  rota- 
tion of  the  axes  through  45  degrees  and  is  given  by  the  nonsingular 
matrix 

T    -i-  U  1 

The  first  property  relates  the  rectilinear  distance  in  E   to 
the  Chebyshev  distance  in  T(E  ).   This  result  is  stated  without  proof 
by  Beckenbach  and  Bellman  [1]  and  by  Francis  [15]  ,  but  it  is  given  and 
proved  here  for  completeness. 

Property  2.5.1.   [|  X  -  p|{   =  J2    ||  T(X)  -  TCP)^  for  any  two  points 
X  and  P  in  E  . 

Proof:  J2    \\   T(x,y)-T(a  ,b)  J|  a  =  slZ   max  jl/v^f  x+y-a-b  |  ,  l/v^|-x+y+a-b|} 

=  max  |(x+y-a-b) ,  (-x-y+a+b) ,  (-x+y+a-b),  (x-y-a+b)| 

=  max  |(x-a+y-b)  ,  (x-a-y+b)  ,  (-x+a+y-b)  ,  (-x+a-y+b)j- 
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=  |x  -  a|  +  Jy  -  b| 


x  -p\\1 


Now  let  P'  =  T(P)  be  the  image  of  a  point  P  under  T.   Define 

S'  =  T(S)  =  JP!  :  P!  =  T(P  ) ,  P   e  s\    . 
I  1    1      i    i    J 

Let  g^  =  %jj2   and  c'.  =  e  AjS     for  1  <  i  <  m.   Consider  the  problem 

fjlx-  -  p:||„  ♦  g^ 


mm     max 
X'  eE   l<i<m 


(2.5.1) 
s.t.     ||  X'  -  Pi  |L  <  c\  1  <  i  <  m. 

Problem  (2.5.1)  will  be  referred  to  as  the  rotated  one  facility 
problem. 

Using  Property  2.5.1,  the  equivalence  of  problems  (2.2.1)  and 
(2.5.1)  is  established  in  the  next  property. 

Property  2.5.2.   X   is  an  optimal  solution  to  problem  (2.2.1)  if  and 
only  if  T(X  )  is  an  optimal  solution  to  problem  (2.5.1).   Furthermore, 
if  z   and  z'  are  the  optimal  objective  function  values  for  problems 
(2.2.1)  and  (2.5.1),  respectively,  then 


z  =  J2   z*  .  (2.5.2) 

o       o 


Proof:   The  proof  is  an  application  of  Property  2.5.1. 
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Writing  problem  (2.5.1)  using  the  definition  of  the  Chebyshev 
distance  gives 


{w.    [max   (|x'    -   a^|,     |y'    -   b!|)]+   gl} 


min  max 

X'  eE9      l<i<m 


s.t.  max   (  |x'    -   a! | ,    |y'    -  b^|   1  <  c^        1  <  i  <  m 


Each  constraint  of  the  above  problem  may  be  written  as  two  constraints 
and  the  order  of  the  maximizations  may  be  interchanged  in  the  objective 
function  to  give 


CS.Ki",-»ii  +  «J'S.{-ii',.--ii*'i}] 


mm 

„.  _   max 
X'€E2 


s.t.       x'  -  a'  <  c! 
i'    i 


(2.5.3) 


1  <  i  <  m  . 


y'  -  b!  <  c! 
i1    i 


From  the  above  problem,  two  subproblems  are  defined. 


mm  ma 


L  k'*'  -  ai'  +  gi) 


x'  6E        ]<i<m 


s.t.  x*    -   a!      <  c!         1  <   i  <  m 

1  l '  l 


(2.5.4) 


and 


min    max  /   i  ,      i     ,  ~) 
i<mlWi'y  "  V  +  sl! 


y'  eE^^  1^ 


s.t.         |y'  -  b! |  <  c!    1  <  i  < 


(2.5.5) 


m. 


Each  subproblem  is  a  one  facility  minimax  location  problem  in  one 
dimension. 
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The  following  property  relates  the  optimal  solutions  of  the 
subproblems  to  the  rotated  problem,  (2.5.1). 

Property  2.5.3.   If  x'  is  an  optimal  solution  to  subproblem  (2.5.4) 

with  objective  function  value  z'  and  if  y'  is  an  optimal  solution  to 

1         o 

subproblem  (2.5.5)  with  objective  function  value  z'   then  (x',y')  is 

2.  o  o 

an  optimal  solution  to  problem  (2.5.1)  with  optimal  objective  function 

value  z'  =  max  (z'.z*). 

o         1  Z 

Properties  2.5.2  and  2.5.3  show  that  the  optimal  solutions  to 
the  one  facility  problem  (2.2.1)  are  obtained  from  the  independent 
solutions  of  the  subproblems  (2.5.4)  and  (2.5.5).   In  the  next  section 
some  properties  of  the  subproblems  are  developed. 

Notice  that  the  transformation  developed  in  this  section  for 
problem  (2.2.1)  in  E   does  not  extend  directly  to  analogous  problems 
in  E„.   The  locus  of  all  points  of  equal  rectilinear  distance  from 
a  point  X  in  E   is  the  surface  of  an  octahedron;  whereas  the  locus  of 
all  points  of  equal  Chebyshev  distance  from  a  point  X  in  E„  is  the  sur- 
face of  a  cube.   Neither  of  these  figures  is  the  image  of  the  other 
under  a  linear  transformation. 

2.6   Subproblem  Development 

This  section  develops  properties  of  the  subproblems  that  form 
a  basis  for  their  solution. 

Each  subproblem  is  first  rewritten  with  simplified  constraints. 

Each  constraint  of  problem  (2.5.4)  may  be  written  equivalently  as  two 

inequality  constraints:  x'  s  a'.  -  c'  and  x'  <  a!  +  c!  for  1  <  i  <  m. 
^     J  li  ii 
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Let 


„,  max  ,    <  <  \ 

C\    =  ^  <r       (a   "  c!  ) 

1    1<i<m   1     l 


mm     ,     ,  . 

2      1<1<TT1      1        1 


Then  problem  (2.5.4)  is  equivalent  to  the  problem 

{w.jx'  -  all  +  g«} 


mm   max 


x'  SE   l<i<m 


s.t.   C^  <  x'  <  c; 


Similarly  for  problem  (2.5.5),  let 


.     max 
D'  =  ^  ^  (b!  -  c!) 
1    l<i<m    l     l 


mm     .     ,  . 
Dl  =  __.  (b  +  c!) 

2    l<i<m   i     l 


and  problem  (2.5.5)  is  equivalent  to 


min    max   f   i  ,      |     ,"l 

|w.|y  -  b!|  +  gW 


y'eE   l<i<m 


s.t.   D'  <  y'  <  D* 
1    J  2 


(2.6.1) 


(2.6.2) 


The  subproblems  (2.6.1)  and  (2.6.2)  provide  a  test  for  the 

feasibility  of  the  original  problem  (2.2.1).   That  is,  problem  (2.2.1) 

has  a  feasible  solution  if  and  only  if  each  subproblem,  (2.6.1)  and 

(2.6.2),  has  a  feasible  solution.   Either  subproblem  (2.6.1)  or  (2.6.2) 

is  inf easible  if  C  <  C  ,  or  D  <  D  ,  respectively. 
&         1      —I 

Problems  (2.6.1)  and  (2.6.2)  have  the  same  structure  so  that 
for  any  property  of  (2.6.1)  there  is  an  analogous  property  for  (2.6.2). 
Thus  it  is  sufficient  to  develop  the  following  properties  only  for 
problem  (2.6.1). 
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Problem  (2.6.1)  without  constraints  is  written  as 


mm    max 


If  x'  is  an  optimal  solution  to  problem  (2.6.3)  and  x'  is 
o  o 

feasible  to  problem  (2.6.1),  then  x'  is  an  optimal  solution  to  (2.6.1). 
The  case  where  x'  is  not  feasible  to  problem  (2.6.1)  is  considered  by 
the  next  property. 


Property  2.6.1.   If  x'  is  an  optimal  solution  to  (2.6.3)  and  x1  is 

: O  O 

infeasible  to  (2.6.1),  then  the  optimal  solution  to  (2.6.1)  is  either 

C1  or  Cl.   Specifically,  if  x'  <  C  ,  C  is  the  unique  optimal  solution 
1     2    F         '      oil 

to  (2.6.1),  or  if  x'  >  C'  ,  C'  is  the  unique  optimal  solution  to  (2.6.1) 

The  proof  of  Property  2.6.1  is  given  after  the  proof  of 
Property  2.6.2. 


Property  2.6.1  shows  that  it  is  sufficient  to  obtain  solutions 
only  to  the  unconstrained  problem,  (2.6.3),  when  solving  problem  (2.6.1) 
The  following  properties  are  developed  for  the  unconstrained  problem. 

Property  2.6.2.   Problem  (2.6.3)  has  a  unique  optimal  solution. 

Proof:   The  linear  programming  problem  for  (2.6.3)  is 
min  z' 


s.t.   x1  +  z'/w.  £  a'.  +  g! /w. 
l     ill 


-x'  +  z'  -  W  2:  -a!  +  g'  w. 
1     111 


x'  ,  z'  unrestricted. 


1  <  i  <  m        (2.6.4) 
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Suppose  that  with  optimal  objective  function  value  z' ,  problem  (2.6.3) 
has  alternative  optima  x'  and  x'  with  x'  4   x'    Then,  x'  and  x'  must 

X        £>  X      ^  1        Jt 

satisfy  some  common  constraint,  that  is,  both  points  lie  on  some  edge 
of  the  feasible  set.   Then,  either  there  exists  an  i  such  that 


x'  +  z'/w.  =  a!  +  g!/w.    and   x'  +  z'/w.  =  a!  +  g'/w 
1     li     ill  2     li     ill 

or   there   exists   a  j    such   that 

-x'    +    z'/w.    =   -a'.    +   g'. /w.    and      -x'    +    z'/w     =    -a'    +   g'/w 
llj  J  J      J  2  1      j  J  j      j 

In   either   case,    eliminating   z'    yields    x'    =   x' ,    a    contradiction. 

X  X      £ 

Proof  of  Property  2.6.1:   Write  the  objective  function  of  problem  (2.6.3) 


as 


V"-^^ -•!!♦«!} 


Then,  f  (x' )  is  a  continuous,  convex  function  in  x'  with  a  unique  min- 
imum, say  x! .  These  properties  imply  that  f  (x' )  is  strictly  increas- 
ing in  the  interval  Tx'  ,°°)  and  strictly  decreasing  in  the  interval 

(-=°,x  ].   Then,  if  x'  <  C'  <  x'  <  C'   f(C')  <  f „  (x1 )  for  all  x'  in 
o  o     1         2    111 

(C' , C' ) .   If   C'  <  x'  <  C'  <  x'  ,  then  f  (x' )  >  f , (C' )  for  all  x' 
x^  1  J*  o  X         12 

in  (C^.CJj). 

Next  a  lower  bound  is  given  for  problem  (2.6.3).   A  bound  was 
obtained  by  Francis  [12]  for  problem  (2.6.3)  with  g!  =  0  and  is  extended 
here  to  the  problem  with  g!  4   0.   For  each  i  and  j  with  1  <  i  <  j  <  m, 


define  the  number 


b  .  =  max 


I w .  w  .  a '.  -  a '.  I  +  w  .  g !  +  w .  g '   ,  g '  ,  g '  \ 

Lw±+w  L  i  3      i    J1     J  i     isjj   &i'  BjJ 
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and  define 


max        f, 

b    *    =    -*-<r~z      1b-  -r    •  (2.6.5) 

st        3<i<T<in    L   ijJ 


n  n   ~   -.        *         ^-     niin        max      f      i    ,  ,  i  ,~1 

Property   2.6.3.       b^  <  ^  ^    ^^  {w. |x«    -   a«  J  +   glj    . 

The  proof   of   Property   2.6.3    is    a   direct  extension  of   that   given 
in    [12]    and   is   omitted.      The  next   property    is    a  new  result    and   estab- 
lishes  that    this    lower  bound    is    always    attained  at   optimality   for 
problem    (2.6. 3) . 

Property    2.6.1.       The  optimal    objective   function  value  of   problem    (2.6.3) 
equals    its    respective   lower   bound. 

Proof:      The   dual    to    the   linear   program,     (2.6.4),    is 


max    (a'+g'    w   )u     +  ..  +  (a'+g'    w   )u    +  (-a'+g'    w   ) v„  +  .  .  +  (-a'+g' .  w   ) v 
1111  mm     mm  1111  mm     mm 

s.  t .  (1)      u„  +  u_  +  .  . .+   u     -  v     -  v     -  ...   -    v      =0 

I  J  m        1        1  m 

(2)  u      w       +  .  .  .  +  u     w     +    v,    w    +  .     .  -i.     v      w    =1 

II  mmll  mm 

(3)  u    ,    v      s  0  1  <   i  <  m. 

i         i 

The  constraints  imply  that  any  basis  matrix  to  the  dual  must  be  of  the 
form 

1     -1 
1  w .   1/w  . , 

for  some  i  and  i.   Constraints  1  and  3  imply  that  u.  =  v.  so  that 

i    J 

constraint  2  gives 
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u./w.  +  v./w.  =  u. (1/w.  +  1/w.)  =  1 
i   i     J  3  i     i      J 


and 


u.  = 


w  w 
i  J 


1    w . +w . 

1  J 


Thus  the  optimal  objective  function  value  has  the  form 

l<i<j<m  lii   Bi   i    j    J   J  J 

=  «****■**■      l w.w.(al  -  a'.)  +  g'.w.  +  g'.w.   ,  g'.,g'.r 

l<i<j<m  Lw.+w.  L  i  J   i    J     i  J    J  iJ    i  JJ 

which  is  the  lower  bound. 

The  next  property  is  an  extension  of  a  similar  result  in  [12] 


Property  2.6.5.   Suppose  b    is  the  lower  bound  for  (2.6.3),  deter- 

st 

mined  by  a'  and  a'    Define 
s      t 


=  — ^—  |w  a'  +  w  a'  +  g!  -  g' I  .  (2.6.6) 

w  +w   |_  S  S      t  t     t      SJ 


X 

St 

s   t 


Then  x   is  the  optimal  solution  to  (2.6.3). 

Proof:   For  a'  <  a'   consider  the  basis  matrix  B  for  the  dual  problem 
s     t 

and  its  inverse  B 

1     -1  I  w  w   /  1/w   1 

b"1  =  -JL*. 

1/w-   1/w  /  Ws+Ut   -1/w^   1 

t      s  /  \     t 

Then  the  optimal  solution  to  the  primal  is  given  by 

(a!  +  gi/w,,  -  a'  +  g'/w  )B_1  =  (xT ,z«) 
ttt     sss         ol 
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which  gives 


'  =  w  w  (a'  -  a'  )  +  g'w  +  g' w   =  b 

1   w  +w  |_stt    s     ts     st.J    st 


s   t 


x'  =  w  a'  +  w  a'  +  g'  -  g' 

o   w  +w  |_  t  t     s  s     t     sj 


s   t 

Since  the  optimal  solution  is  unique,  it  is  given  by  (2.6.6).   Note 

that  if  a'  =  a'  ,  equation  (2.6.6)  gives  x  __  =  a'  =  a'  ,  and  in  this  case 
s    t  st    s    t 

the  bound  is  given  by  the  max  (g' ,  g' )  =  b   or  b   .   Then  for  the 

case  b   ,  for  example,  the  dual  basis  gives  x'  =  a'   and   z'  =  g' . 
ss  o    s        1    s 

2.7   Solution  Procedure 

The  properties  developed  in  the  preceding  section  for  the 
subproblem  (2.6.1)  and  the  analogous  properties  that  exist  for  sub- 
problem  (2.6.2)  are  used  to  determine  a  solution  procedure  for  both 
subproblems.   Using  the  optimal  solutions  to  the  subproblems,  a  proce- 
dure is  given  for  finding  all  alternative  solutions  to  the  one  facil- 
ity problem  in  two  dimensions,  (2.2.1). 

Subproblem  (2.6.1)  is  solved  by  computing  the  lower  bound, 

(2.6.5),  and  determining  the  value  x'  using  (2.6.6).   If  x'  is  feasible 

o  o 

to  problem  (2.6.1),  then  the  unique  optimal  solution  is  x'  =  x'  ,  and 

z'  =  b     If  it  is  not  feasible,  the  unique  optimal  solution  is 
1    st 

x1  =  C  or  x'  =  C'  as  determined  by  Property  2.6.1,  and 
o     1     o    2 


,     max 
Zl  "  l<i<m 


{wjc-   -  all  +  gj 


for  j   -    1  or  2 ,  respectively. 
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Similarly,  subproblem  (2.6.2)  is  solved  by  computing  the  lower 
bound 

b  „  =    ™^L  {— —  rw.wjb!  -  b'.|  +  w.g!  +  w.g'.l  ,  g!  ,g'\    . 

st    3<i<j<m  Lw.+w  L  i  j   i    J      J  *     i  jJ  i   J-1 


Then  compute  y   using 


=  — - —  w  b'  +  w  b*  +  g!  -  g' 
t   w  +w  L  s  s     t  t    &t     sj 


ys- 

s   t 

If  v'   is  feasible  to  (2.6.2),  the  unique  optimal  solution  to  (2.6.2)  is 
st 

y'  =  y'    and   z'  =  b   ;  if  not,  the  unique  optimal  solution  is  y'  =  D1 
3o        Jst        1    st'  Jo     1 


or  y'  =  D'   and 
O     2 


,     max 
Z2  ;:  ]<i<m 


{wJd-  -  b«j  +  g:} 


for  j  =  1  or  2 ,  respectively. 

By  Property  2.5.3,  the  optimal  objective  function  value,  z' , 

o 

to  the  rotated  problem,  (2.5.1),  is  given  by  z'  =  max  (z'  z')  where 

O  i.   ^ 

z'  and  z'  are  the  optimal  objective -function  values  to  (2.6.1)  and 

(2.6.2),  respectively. 

The  point  (x' ,y')  is  one  optimal  solution  to  (2.5.1)  by 
o   o 

Property  2.5.3.   Three  cases  need  to  be  considered  to  determine  all 

the  alternative  solutions. 

Case  1:   z'  =  z'  =  z'.   Then  (x*  ,y' )  is  the  unique  optimal 
o     1    2  o  ^o 

solution  to  (2.5.1). 


Case  2:   z'  =  z'  >  z' .   Then  for  all  y'  that  satisfy  the 
o     1    2  J  J 


inequalities 
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*.|y'  -  b'l  +  g'  <  z'      1  <  i  <  m 

1         1       1      o 

and  (2.7.1) 

D'  <  y'  <  Dl 

1    J  2 

the  points  (x' ,y')  are  also  optimal  to  (2.5.1).   For  1  <  i  <  m.  the  i 
o  ' 

inequality  in  (2.7.1)  may  be  written  as  two  linear  inequalities  so  that 
(2.7.1)  is  equivalent  to 


b!  -  (z'  -g' )  /w.  <  y'  <  b:  +  (z'  -g'.  )/w.      1  <  i  <  m 

1        Oil  1        Oil 

D^  <  y'  <  D;,. 
These  inequalities  determine  an  interval,  (y1  ,y' )  ,  of  alternative  y' 

J.     •- 


where 


y'    =   max  Td'         ™**        |b'    -    (z1    -   g')/w  }  "I 
Jl  L   1      l<i<m      I   i  o        bi        iJ    J 


(2.7.2) 


y2  =  mi 


■fc.,^fo  ♦<■;-«!>**}] 


Then  the  points  (x'.y')  on  the  line  segment  determined  by  (x',y')  and 

(x' .y^)  ai"e  the  optimal  solutions  to  (2.5.1). 
o  2 

Case  3:   z1  =  z'  >  z1 .   Then  determine 
o    2    l 

x'  =  max  Tc*  ,  ™&*     |a*  -  (z'  -  g')/w  ]■  1 
1       L  1   l^i^m  I  i     o   Bi    ii  J 


(2.7.3) 


X2 


2  =  ^[c2'^4k  +  <z;-^/wi}] 


Then  the  points  (x*  ,y' )  on  the  line  segment  determined  by 
(x^,yM  and  (x^,y^)  are  all  optimal  solutions  to  problem  (2.5.1). 
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By  Property  2.5.2  the  optimal  solution  to  the  problem  (2.2.1) 

is  obtained  by  taking  the  inverse  rotation,  T   ,  of  the  solutions  to 

(2.5.1).   In  Case  1,  (x  ,y  )  =  T   (x' ,y' )  is  the  unique  solution  to 

o   o         o   o 

(2.2.1).   For  Cases  2  and  3  with  a  line  segment  of  solutions,  the  image 

of  a  line  segment  under  a  nonsingular  transformation  of  E  onto  E~  is 

a  line  segment.   Thus  for  Case  2,  the  points  (x  ,y  )  on  the  line  seg- 

o   o 

ment  determined  by  (x  ,y  )  =  T   (x' ,y')  and  (x  y  )  =  T   (x' ,y' )  are 

11         ol        zz         o   z 

all  optimal  solutions  to  (2.2.1).   Similarly  for  Case  3,  the  optimal 

solutions  are  on  the  line  segment  determined  by  the  points 

(x  y  )  =  T_1(x'  y')   and   (x  y  )  =  T_1(x'  y'). 
11  lo  z   z  zo 

2.8  Additional  Properties 

Using  the  above  developments,  additional  properties  of  the  one 
facility  problem  are  given. 

A  lower  bound  is  determined  for  the  problem  (2.2.1)  and  it  is 
shown  that  for  the  unconstrained  problem,  (2.2.2),  this  lower  bound  is 
always  attained  at  optimality.   A  similar  bound  is  obtained  in  [12] 
and  is  extended  here  to  the  case  with  g.  /  0.   Define  the  numbers  B  . 
for  each  i  and  j  such  that  1  <  i  <  j  <  m  by 

B.  .    =   max  \    w.w  .  ||  P.  -  P  .  j!       +  w.g.    +    w  g        ,  g    ,g    f 

ij  I  wi+w     L   i   j"     i        JM1         Ji  i&jJ'6i'sjJ 


and  define 


max  f 

st  "'    l<i<j<m    I   ij 


J 
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Property  2.8.1.   If  z   is  the  optimal  objective  function  value  of 

problem  (2.2.1),  then  B  .  <  z  . 

st    o 

Proof:   By  a  direct  extension  of  the  proof  given  [12],  for  any  X  in  E„, 


Therefore  inequality  (2.8.1)  holds  for  all  feasible  points,  X,  to 
problem  (2.2.1)  and  for  the  minimum  point,  X  ,  in  the  feasible  region. 


Hence, 

B   <   maX 
st  ~  l<i<m 


fil|Xo-Pi'fl  +  «i}  =  V 


Property  2.8.2.   The  optimal  objective  function  value  of  the  uncon- 
strained problem,  (2.2.2),  equals  the  lower  bound  B  . 

Proof:   Let  z'  be  the  optimal  objective  function  value  of  the  rotated 
o 

problem,  (2.5.1),  without  constraints.   Then  by  Property  2.5.3, 

z;  =  raax  {2i'z2; 

where  z'  and  z'  are  the  respective  lower  bounds  of  the  subproblems 
(2.6.1)  and  (2.6.2)  without  constraints.   Then 


,  max 

z1  =  max  /  ,^-.^^- 
o       |  l^i<j<m 


{w"Tw-[}'iwjlal  "  aj'  +  Vj  +  wjg0'gi'gJ  • 


max 
l<i<j<m 


{— —  Iw.w.lb'.    -   b'.|    +   w.g'.    +  w.g!    |  ,    g!.g,}> 
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Property  2.5.2  implies  that  z'  =  z  ,V2  .  Applying  the  transformation 
T_1,  where  a!  =  (a.  +  b.)A/2  and  b!  =  (-a.  +  b.)/«/2,  and  multiplying 
through  by  ^2  gives 

f       max         f — 1 —  r  I        +b.-a.-b.|+w.g.+w.g.       ,    g    ,g    \    , 

1  l<i<j<m   Iw.+w.  L  i   3      i        i        J        J  x   J        J    XJ  x      JJ 

V  i      J 

j — - —     w  w.l -a.  +b.  +  a  .  -b.|  +w.g.  +wg.      ,  g    ,g        / 


z     =   max 
o 


max 
]<i<j<m 


i      3 


■\ ww.   max        a.+b.-a.-b., 

i<j<m    U.+w.  L   i   J  L      i        i        3        3 

| -a.  +  b.  +  a  .  -b.|]  +   w.g.  +  w  g         ,    g. ,g.f 


l£l<£3         i 


max 
l<i<j<m 


i — 1 —  [ww    f|a    -a!+|b.-b.|~|+w.g.+w.g.       ,g.  ,g.f 


i      J 


max        ;_L_[W        ||p     _p    (I      +w.g.+w.g.l,g.,g.} 
i<j<m    Lw.+w.   L  i   J         l        jHl         l   J        J    iJ        i      3-> 


l<i<j<m 


max   J   1 
l<i<j<m  I  ijJ 


The  above  results  are  now  used  to  develop  a  procedure  to 
determine  level  sets  and  contour  sets  for  problem  (2.2.1). 

For  the  function  f(X)  and  for  a  real  positive  number  k,  a 
level  set  with  respect  to  k,  denoted  by  S  ,  is  defined  by 


Sk  =  {x  e  E2  :  f (X)  <  kj  . 
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A  contour  set,  for  the  function  f(X),  with  respect  to  k  is  denoted  by 


T  ,  and  is  defined  by 


Tk  =  {x  e  E9  :  f(X)  =  kj 


Denote  the  set  of  feasible  points  to  problem  (2.2.1)  by  F. 
Then  a  level  set  with  respect  to  k  of  problem  (2.2.1)  is  F  fl  S   and 
a  contour  set  of  problem  (2.2.1)  with  respect  to  k  is  F   T  . 

The  level  set,  F  fl  S.,  is  the  set  of  points  X  such  that 

w   ||  X-P   ||      +g±  <  k  l<i<m  (2.8.2) 

II  X-P    II  <  c.  1  <   i  <  m.  (2.8.3) 

"  l     1  l 

The  points  X  satisfying  (2.8.2)  and  (2.8.3)  may  be  determined 
using  the  preceding  development.   That  is,  transform  the  inequalities 
by  T,  decompose  the  resulting  inequalities  and  determine  the  points 
(Xl.y*)  (x'  y')  (x*  y')  (xl, yl),  using  expressions  (2.7.2)  and  2.7.3) 
with  k  replacing  z' .   Then  the  corners  of  the  set  of  solutions  are 
determined  by  taking  the  inverse  transformation  of  these  points. 

Alternatively,  the  points  X  in  the  level  set  F  M  S  may  be 

determined  directly  from  the  expressions  (2.8.2)  and  (2.8.3).   For 

each  i,  the  i    inequality  in  (2.8.2)  and  (2.8.3)  is  equivalent  to  four 

linear  inequalities.   Define  G.  (k)  =  (k-g.)/w..   Then  the  equivalent 

linear  inequalities  are 

x  +  y  S  a.  +  b.  -  G.  (k) 
•'111 

x  -  y  S  a.  -  b.  -  G.  (k) 
ill 

-x  +  y  ^  -a.  +  b.  -  G.  (k) 


1  <  i  <  m 


-x  -  y  >  -a.  -  b.  -  G.  (k) 

'      ill 
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and 


1  <  i  <  m 


x+y^a      +   b.    -   c 

J  1  i  l 

x   -   y    >   a      -b      -c 

1  i  l  i 

-x  +  y^-a     +b      -c 
l  l  i 

-x-yS-a      -b      -c. 
l  i  i 

All    but    four    inequalities   from  the   above   are  redundant.       To   elim- 
inate  the  redundant    inequalities   define   the  numbers,  dependent   on  k, 

/i  •,  max      •  ,_  „    ,,  ■>  ^         max 

e1  (k)    =    max  (a.    +   b.    -   G.  (k))  ,    ,__,      (a.    +   b.    -   c.) 

1  Ll^i^m        ill  l<i<m        i  i  i  J 

e.(k)    =   maxl"^      (a.    -   b.    -   G.(k)),      ™a*      (a.    -  b.    -   c.)"l 

2  Ll^i<m        ill  l<i<m        l  i  i  J 

e    (k)    =   max  I"™1*      (-a.    +   b.    -  G.(k)),      ™a*      (-a  +    b.    -   c  1 

3  Ll^i^m  ill  l<i<m  i  ij 


(2.8.4) 


e    (k)    =   maxf™^      (-a.    -  b.    -  G.(k)),      ™*      (-a.    -   b.    -   c.)~] 
4  LlSiSra  ill  l<i<m  l  i  i  J 


Then   the   four    inequalities 

x  +   y   >   ex(k) 

x  -  y  >  e  (k) 
-x  +  y  >  e  (k) 
-x  -  y  S:   e4(k) 

are  equivalent  to  the  inequalities  (2.8.2)  and  (2.8.3),  and  the  set  of 
points   X  =  (x,y)  satisfying  (2.8.5)  are  the  points  in  the  level  set 

f  n  s, . 

k 


(2.8.5) 
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It  can  be  shown  that  the  contour  set,  F  I  T ,  is  the  boundary 

K 

of  the  level  set  F  fl  S,  .   Hence  to  determine  the  contour  set,  use  the 

k 

expressions  (2.8.4)  to  determine  e  (k) ,  e  (k) ,  e  (k)  and  e  (k) .   Then 
write  the  expressions  (2.8.5)  as  equalities,  as  follows: 

x  +  y  =  e  (k) 

x  -  y  =  e  (k) 

-x  +  y  =  e  (k) 

-X  -  y  =  e„(k)  . 
4 

Then  the  set  of  all  X  =  (x,y)  satisfying  these  equations  are 

the  points  in  the  contour  set  F  11  T.  .   Notice  that  if  k  is  chosen  as 

K 

the  optimal  objective  function  value,  F   T   is  the  set  of  optimal 
solutions  to  problem  (2.2.1). 

In  the  previous  development  it  was  assumed  that  w.  >  0  for 
all  P.  in  S.   The  following  procedure  considers  the  case  where  w.  =  0 

for  some  P.  in  S.   Let 

l 

I  s  |i  :  P  e  S,  w.  >  oj  . 

Then  the  problem  (2.2.1)  may  be  written  as 

min      ("max  /   n  „   „  ||  max  /  1 

XeE.  maX  Lew  Wl  X  "  Pi"l  +  SJ    '  l<i^n  lgi/  J 

iew 


s.  t.      II  X  -  P.  II,  <  c.     1  <  i  <  m  . 
i"  1     i 

The   subproblem 

min     max  /      n  n  r 

XeE2      iCwK'IX-PjIj    +   gJ 

s.t.  !|  X   -  P.  ||,    <  c.  1  <   i  <  m 

'  l"  1  l 
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ay  be  solved,  using  the  procedure  determined  above,  for  its  optimal 


m 


objective  function  value  z  .   Let 


max   f  1 


go  =:  l<i<m 
i£W 


If  z   s  g  ,  then  the  solution  to  problem  (2.2.1)  is  given  by  the  solu- 
o    o 

tion  to  the  subproblem  over  the  set  W.   If  z  <  g  ,  than  all  X  such 

o    o 

that 

w.l!  X  -  P. ||,  +  g.  <  g     for  all     i  in  W  (2.8.6) 

1         1   1      1     o 

and       \\   X  -  P.JL       ^  c     for   all    1  <  i  <  m  (2.8.7) 

l  1  l 


are  solutions  to  the  original  problem  with  optimal  objective  function 

value  g  . 
o 

Because  of  the  similarity  between  (2.8.6),  (2.8.7)  and  (2.8.2), 

(2.8.3),  the  solutions  may  be  determined  using  the  procedure  developed 

for  the  level  sets.   That  is,  in  (2.8.4)  set  g   =  k  and  take  the  max- 

o 

imum  of  the  first  terms  over  the  index  set  W  to  determine  e  (g  ) , 

1   o 

e„(g  )»  e„(g  )  and  e  (g  ).   Then  the  points  X  =  (x,y)  satisfying  the 
inequalities 

x  +  y  2:  e1(gQ) 

x  -  y  £  e2(gQ) 

-x  +  y  >  e  (g  ) 
o   o 

-x  -  y  >  e4(gQ) 

are  solutions  to  problems  being  considered  with  g  >  z  . 

o   ■  o 
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A  lower  bound  for  problem  (2.2.1)  assuming  that  w  =  0  for 
some  i  is  also  given  by  this  development.   Let  B   be  the  lower  bound 
for  the  subproblem  (2.8.2).   Then  a  lower  bound  for  problem  (2.2.1) 

with  some  w.  =  0  is 

i 


[Bsf  iSmkl] 


B  =  max 

L  st   J 

iew 

From  Property  2.8.1  and  the  above  discussion,  B  is  a  lower 
bound  for  problem  (2.2.1). 


CHAPTER  3 
A  MULT I FACILITY  PROBLEM 

3.1   Introduction 

The  problem  considered  in  this  chapter  is  to  locate  several 
new  facilities  with  respect  to  a  finite  number  of  existing  facilities 
so  as  to  minimize  the  maximum  cost,  where  the  costs  are  proportional 
to  the  rectilinear  distances  between  new  and  existing  facilities  or 
proportional  to  the  rectilinear  distances  among  new  facilities. 
The  problem  also  includes  upper  bound  constraints  on  the  Tectilinear 
distances  between  new  and  existing  facilities  and  among  new  facilities. 

The  transformation  T,  developed  in  Chapter  2,  is  applied  to 
this  problem  and  gives  a  decomposition  into  two  independent  one- 
dimensional  minimax  subproblems. 

If  there  are  costs  or  constraints  among  new  facilities,  each 
subproblem  is  formulated  as  an  equivalent  network  weighted  flow  problem 
and  is  solved  by  determining  a  simple  chain  of  maximum  cost  to  weight 
ratio  in  the  network.   If  there  are  no  costs  and  no  constraints  among 
new  facilities,  each  subproblem  decomposes  into  one  facility  subprob- 
lems which  are  solved  using  the  procedures  of  Chapter  2. 

All  optimal  solutions  to  the  original  multif acility  problem  are 
obtained  from  the  solutions  to  the  subproblems. 
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3.2   Formulation 

The  multif acility  minimax  location  problem  is  now  formulated 
using  the  notation  developed  in  Chapter  2. 

Let  X.,  1  <  j  <  n,  denote  the  n  points  in  E   to  be  determined 

J  - 

for  the  location  of  n  new  facilities.   Let  S  be  a  finite  set  in  E9  of 

m  points  P. ,  1  <  i  <  m,  representing  the  location  of  existing  facilities. 

For  each  point  P   and  X   let  d..  be  a  positive  constant  which  is 
i      J      Ji 

an  upper  bound  on  the  rectilinear  distance  between  P.  and  X..   Also 

let  c   be  a  positive  constant  which  is  an  upper  bound  on  the  rectilin- 
jk 

ear  distance  between  the  pair  of  new  facility  locations  X.  and  X   for 
1  <  j  <  k  <  n. 

The  objective  function  is  a  function  of  the  n  variables  X.  in 
E  ,  1  <  j  <  n,  and  is  the  maximum  of  all  expressions  which  are  propor- 
tional to  the  rectilinear  distance  between  each  pair  of  new  facility 

locations  X   and  X   and  between  each  new  facility  location  X.   and  each 
j      k  J 

existing  facility  location  P..   Let  u.   be  a  real  nonnegative  constant 

of  proportionality  for  the  rectilinear  distance  between  each  point  X. 

and  P   for  1  <  j  <  n,  and  1  <  i  <  m.   Let  v   be  a  real  nonnegative 
i  Jk 

constant  of  proportionality  for  the  rectilinear  distance  between  each 

pair  of  points  X.  and  X,  for  1  <  j  <  k  <  n. 
J      k 

Then  the  multif acility  minimax  location  problem  to  be  studied 
in  this  chapter  can  be  stated  as 
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Xl,.  m,nXneE2   maX  [l^kk  fjk"  Xj  -Xk"l}'    ^n  (w.Jl  X.  -P   ||     }  1 
1  n      2  ]<i<m    I   ji        J         i     U    J 

s.t.  ||  X.    -   X    ||       <   c.,  1  <  j    <  k  <  n 

1      J  k'l  jk 

(3.2.1) 
I!  Xj    -  P±||1    <dj±  1<  j  <  n 

1  <  i  <  ra. 

As  before,  the  unconstrained  multif acility  problem,  which  will 

be  referred  to  as  problem  (3.2.2),  is  a  special  case  of  problem  (3.2.1) 

obtained  by  letting  the  upper  bounds  c .,  and  d..  be  large. 

Jk      ji 

Thus  each  of  the  n  new  facilities  is  to  be  located  with  respect 

to  the  m  existing  facilities  and  also  with  respect  to  the  other  new 

facilities.   The  location  of  the  point  X.  may  depend  on  the  location  of 

some  point  X  because  of  the  terms  involving  v    in  the  objective  func- 
k  j  k 

tion  or  because  of  the  constraint  with  upper  bound  c .,  .   This  depend- 

Jk 

ence  among  new  facility  locations  is  made  more  explicit  by  the  follow- 
ing definitions  which  were  developed  in  [14]  for  an  unconstrained  total 
cost  multif acility  location  problem  and  are  extended  here  to  problem 
(3.2.1). 

New  facility  locations  X.  and  X,  are  said  to  be  linked  if  v  , 

J      k  jk 

is  positive  or  if  there  exists  a  constraint  for  X.  and  X,  with  c  , 

J      k       jk 

finite,  and  are  said  to  be  not  linked  if  v.,  is  zero  and  c    is  infin- 

Jk  Jk 

ite,  that  is,  there  is  no  constraint  for  X.  and  X  . 

J      k 


The  index  set,  V(j),  of  positive  v    for  each  1  <  j  <  n  -  1  is 

Jk 


defined  by 


V(j)  =  [k  :  1  <  j  <  k  <  n,  v   >  0} 

J1'* 
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The  index  set,  C(j),  for  which  there  exists  finite  c    for 

Jk 

1  <  J  <  n  -  1  is  defined  by 


C(j)  =  {k  :  1  <  j  <  k  <  n,    0  <  c   <  -}  . 

JK 

It  is  assumed  that  each  new  facility  location  X.  is  linked  with 

J 

at  least  one  other  new  facility  location.   If  this  assumption  is  not 
made,  then  the  location  of  X.  could  be  determined  independently  of  the 
other  new  facility  locations  by  considering  a  separate  problem,  and 
problem  (3.2.1)  would  not  be  formulated  properly. 

In  addition,  a  new  facility  location  X.  and  an  existing  facil- 
ity location  P   are  said  to  be  linked  if  w..  is  positive  or  if  there  is 
J  i  Ji 

a  constraint  for  X   and  P   with  d.   finite,  and  are  said  to  be  not 

j       i       Ji  

linked  otherwise.   The  index  set,  W(j),  of  positive  w   lor  each 
1  <  j  <  n  is  defined  by 

W(j)  =  [i  :  1  <  i  <  m,  w   >  0}  . 

The  index  set,  D(j),  for  which  there  exist  finite  d..  for  1  <  j  <  n  is 


defined  by 


D(j)  =  {i  :  1  <  i  <  m,         0  <  d   <  «} 


The  new  facility  location  Xj   is  said  to  be  chained  if  there 

exists  a  sequence  of  new  facility  locations  Xj  ,Xj0, .  .  .  ,Xj  ,  such  that 

X-i   and  X.:     in  this  sequence  are  linked  for  1  <  k  <  p  -  1,  and  the 
J  k      J  k-i- 1 

new  facility  ar  X-:   is  linked  with  an  existing  facility  at  P.  for 

'      Jp  i 

some  i.   It  is  assumed  that  all  new  facility  locations  are  chained, 
for  otherwise  the  unchained  locations  could  be  placed  at  the  same 
point  in  the  plane. 
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The  chaining  assumptions  are  used  in  [14]  to  prove  that  there 
exists  an  optimal  solution  to  the  total  cost  raultif acility  problem. 
In  an  analogous  fashion,  the  above  assumptions  imply  that  if  there  is 
a  feasible  solution  to  problem  (3.2.1),  there  is  an  optimal  solution. 

With  these  definitions  the  multif acility  minimax  location 
problem  (3.2.1)  becomes 

Xl,.°lnx  £b2  - LSU  hJi  vMJ  ■  £L   {'J  VillJ  1 

kev(j)  iew(j) 

(3.2.3) 

s.t.     \\  X.    -   xj^  <  cjR     1  <  j  <  n  -  l,  k  e  C(j) 

"  xj  ■  PA  -  dji      *  - j  - n>      i  e  D(j)- 

In  order  to  simplify  the  notation  and  presentation  in  the 

remainder  of  this  chapter,  the  development  is  given  for  problem  (3.2.1) 

with  the  assumption  that  the  constants  v.,  are  positive  and  the  c    are 

Jk  jk 

finite  for  all  1  <  j  <  k  <  n,  and  that  the  constants  w   are  positive 

Ji 

and  the  d..  are  finite  for  all  1  <  i  .<  n,  1  <  i  <  m.   Differences 
Ji 

between  this  development  and  a  development  for  the  more  general  problem 
(3.2.3)  will  be  pointed  out  in  the  sequel. 

Notice  that  for  n  =  1,  the  multif acility  problem  (3.2.1)  reduces 
to  a  one  facility  problem. 

A  generalization  of  problem  (3.2.1)  may  be  formulated  to  include 

constants  g .  .  corresponding  to  X   and  P   for  1  <  i  <  n  and  1  <  i  <  m: 
Ji  j       i 

and  h    corresponding  to  X   and  X   for  1  <  j  <  k  <  n: 
Jk  j      k 
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^hJv^W-^kiiiv'A^i}] 


min  ma; 

„  max   „  _.  .  _  . 

l<j<n 


s.t.     \\X.     -   Xk\\1    <  cjk      1<  j  <  k<  n 

(3.2.4) 
\\x.     -   Pi|l1  <  d       1  <  i  <  m,   1  <  j  <  n. 

Most  of  the  development  in  this  chapter  applies  to  problem  (3.2.4)  with 

specific  assumptions  for  the  g..  and  h  .,  .   A  discussion  of  problem 

Ji      jk 

(3.2.4)  in  terms  of  the  properties  and  procedures  developed  for  prob- 
lem (3.2.1)   is  included  in  a  final  section  of  this  chapter. 

3.3   Applications 

Applications  of  multif acil ity  location  problems  occur  in  the 
same  contexts  as  those  discussed  in  Chapters  1  and  2. 

As  a  specific  example,  suppose  that  n  new  machines  are  to  be 

located  in  a  plant  with  respect  to  m  existing  machines  located  at 

points  P. ,  1  <  i  <  m.   Travel  on  a  system  of  orthogonal  aisles  is  to 

take  place  between  the  new  machines  and  existing  machines,  and  among 

the  new  machines.   Let  w..  be  an  estimate  of  the  number  of  trips  per 

Ji 

month  between  new  machine  j  and  existing  machine  i.   Let  v   be  an 

estimate  of  the  number  of  trips  per  month  between  new  machines  j  and  k. 

It  is  desired  to  locate  the  n  new  machines  so  that  the  maximum  estimated 

travel  distance  per  month  is  minimized.   This  problem  may  be  formulated 

as  problem  (3.2.2). 

In  other  applications  the  number  w..  and  v.,  may  find  different 

Ji      Jk 

interpretations  as  constants  per  unit  of  distance  between  facilities. 
For  example,  it  may  be  desired  to  minimize  the  maximum  cost  of  travel 
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between  new  and  existing  facilities  and  among  new  facilities;  then  the 
w   and  v   may  be  interpreted  as  costs  per  unit  of  distance.   These 
numbers  may  also  find  interpretations,  other  than  in  the  example,  as 
weights  or  levels  of  interaction  between  facilities. 

3.4  Related  Work 

First  it  is  noted  that  the  analogous  total  cost  multif acility 
location  problem  has  received  attention  previously  for  both  Euclidean 
distances  [14,30]  and  rectilinear  distances  [11,2].   Also,  the  Euclidean 
distance  multif acility  minimax  problem  is  the  subject  of  a  paper  by 
Elzinga  and  Hearn  [6] ,  and  Love  and  Wesolowsky  [24] . 

There  appear  to  be  only  two  reports  extant  concerning  the 
multif acility  problem  explicitly  using  the  rectilinear  distance.   The 
unconstrained  problem,  (3.2.2),  was  formulated  independently  by  Francis 
[16]  and  Wesolowsky  [30].   In  [16]  it  is  shown  that  the  unconstrained 
problem  may  be  written  as  an  equivalent  linear  program.   Wesolowsky  [30] 
also  makes  that  observation,  but  only  for  the  unconstrained  problem  in 
one  dimension.   He  then  suggests  a  solution  procedure  using  parametric 
linear  programming. 

3.5  The  Transformation  and  Decomposition 

The   transformation  T  developed   in  Chapter   2    is   now  applied  to 
problem    (3.2.1)    to  obtain   the  rotated   problem. 


mm 
v.  vi  eF     maX 

X1'""VE2 


Jv^U^faJIV'lU] 


max 
.l<j<k<r 

l<i<m 


■-t-  IIX'     -*JL<c'k        1<  j    <k<n  (351) 

||x'     -P'll       <  d'  1  <  j  <  n,    1  <   i  <  m 

II     j  ill  „  JX 


where   X'.  -  T(X  ■»      „■  ,  FT      ^      .,.  r— 

J  J}»    cjk  =   cjk/v2  and     dji  =   djiA/^ 
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The  following  property  is  an  extension  of  Property  2.5.2  and 
shows  that  problems  (3.2.1)  and  (3.5.1)  are  equivalent. 

Property  3.5.1.   The  set  {X  ,  ...,X  }   is  an  optimal  solution  to  probl 

(3.2.1)  if  and  only  if  [T(X  ) , . . . ,T(X  )}  is  an  optimal  solution  to 

problem  (3.5.1).   Furthermore,  if  z   and  z'  are  the  optimal  objective 

o      o 

function  values  of  problem  (3.2.1)  and  (3.5.1),  respectively,  then 

z   =  z'  </2    . 
o    o 

Analogous  to  the  development  of  the  one  facility  problem,  two 
subproblems  are  defined  from  the  rotated  problem  (3.5.1)  as  follows: 


em 


mm  max 

x'    .  ...x'ez.,  maXLl^J<k<n    L'j] 
l  n      l 


P     max        f 


s.t. 


x'.    -   x' 

J  k 


x'.    -   a' 

J  i 


and 


min 


max 


y^ y;^       _i<j<k<n  I  jk 


r  max  (v 

_l<j<k<n   I  jl 


s.  t. 


y     -  y 


k  - b: 


x- -k:|},"ax   {w..|x-  -a:|}"| 

j        k '  J       l<j^n    I  j  i  '    i        i'J    J 


Jk 


l<i<m 

1  <  j   <  k  <  n 


<  d'.  1  <  j  <  n,       1  <   i  <  m 

Ji 


i-'il}.S{'jil'i-bil}] 


max 
l<j<n 
l<i<m 


-  c  , 
Jk 

<  d'.  . 
Ji 


1  <  j  <  k  <  n 

1  <  j  <  n,   1  <  i  <  m. 


(3.5.2) 


(3.5.3) 


Each  subproblem  is  a  one-dimensional  multif acility  minimax  location 
problem.   The  following  property  is  an  extension  of  Property  2.5.3  and 
relates  the  optimal  solution  of  the  subproblems  to  the  rotated  problem 
(3.5.1). 
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Property  3.5.2.   If  [x   ,...,x   }  is  an  optimal  solution  to  subproblem 

(3.5.2)  with  optimal  objective  function  value  z'  and  if  (y   ,...,y   } 

1        '1       n 

is  an  optimal  solution  to  subproblem  (3.5.3)  with  optimal  objective 

function   value  z ' ,  then  {(x,  ,y   ),..., (x   ,y   ) }  is  an  optimal  solu- 
2  11         n   n 

tion  to  problem  (3.5.1)  with  optimal  objective  function  value  z'  = 
max  (z',z^). 

Thus  Properties  3.5.1  and  3.5.2  show  that  the  optimal  solutions 
to  the  multif acility  problem  (3.2.1)  are  obtained  from  the  independent 
solutions  to  the  subproblems  (3.5.2)  and  (3.5.3).   In  the  next  section, 
the  subproblems  are  formulated  as  network  problems  and  properties  are 
developed. 

3. 6   A  Network  Formulation 

In  this  section  it  is  shown  that  each  subproblem,  (3.5.2)  and 
(3.5.3),  may  be  formulated  as  equivalent  network  weighted  flow  problems. 
Properties  of  the  network  problems  are  developed  that  form  a  basis  for 
a  solution  procedure.   The  subproblems  have  the  same  structure  so  that 
subsequent  properties  are  developed  only  for  problem  (3.5.2),  with 
analogous  properties  existing  for  (3.5.3). 

The  constraints  between  new  and  existing  facilities  in  subprob- 
lem (3.5.2)  are  simplified  by  using  the  development  in  Section  2.6  and 
defining 

,,    max 
d.  =  ,  -.  -   (-d*.  .  +  a'. ) 

D\   =   -5;   Cd»   +  a!) 

j        l<i<m   ji     i 

for  1  <  j  <  n.   Then  the  subproblem  (3.5.2)  is  equivalent  to  the  problem 
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min  max 

max 


kJv-tD.^ki'-s  "I"}] 


x-    ...  ,x'€El    ""*"  |_l<j<k<n 
1  n      1  l<i<m 


s.t.  x-    -   x,'      <  c'  1  <  j   <  k  <  n 

1    J  V  jk 


(3.6.1) 


d'    <   x*   <  D'.  1  <  j  <  n. 

J  J  J 


From  the  chaining  assumptions  made  in  Section  3.2,  if  problem 
(3.6.1)  has  a  feasible  solution.it  has  an  optimal  solution. 

Based  on  the  discussion  in  Section  2.4,  problem  (3.6.1)  is 
equivalent  to  the  following  linear  program: 

min  z' 


s.t.    (la)    x-  -  x'  *   -c'R 

(lb)   -x-  t  x;  *   -cjfc 

(2a)    x*  ^   d' 

J  J 

(2b)   -x'.  £   -D'. 

J  J 

(3a)    x'  -  x,'  +  z*  /v.,     ^   0 
J     k       jk 

(3b)   -x'  +  x*  +  z'/v.,     2:   0 
J     k       Jk 

(4a)    x'  +  z'  ,v.  .         2   a'. 
3  Ji  * 

(4b)   -x*  +  z'  w  ^  -a' 

J       Ji  i 


1  <  j  <  k  <  n 


1  <  j  <  n 


1  <  j  <  k  <  n 


1  <  j  <  n,   1  <  i  <  m 


z'  ,  x'  unrestricted  1  <  j  <  n.        (3.6.2) 

j 
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A  dual  problem  to  (3.6.2)  is  written  as  problem  (3.6.3) 
The  dual  variables  in  (3.6.3)  are  defined  in  Table  3.1. 


ma 


m  n 


m 


-)    c'F+d'.F.     tValf1     -)    c'F       -D'.F    .-Va'f1 
^     jk   jk        j    jt     L     i   jt     L     jk  kj        j    sj     L     i    sj 


k>j 


i=l  k>j 


i=l 


n 


n  m 


n 


n  m 


s.'t. 


dj)  y  f  +f.,  +  Vf.l  +  y f1  -  Vf, .-f  . -  y f . - y f1. =o 

L,        Jk        jt       U     jk       L     Jt       L     kj         sj        ^     kj        L      sj 


k=l 

k^j 


n       n 


k=l  i=l  k=i 

k^j  k^j 


k=l  i=l 

k^j 


n       m 


n        n 


1  <  j  <  n 

n        m 


(2)  y  yf../v..  +  y  Y**/*.. +  y  y*7y.1  +  y  y^vw  =1 

La       U     jk     jk       L       L     jt      ji       L       L     kj      jk       L       L     sj      ji 


j=l  k>j 


j=l   i=l 


j=l  k>j 


j=l   i=l 


(3)   all  variables  nonnegative 


(3.6.3) 


n 


n   m 


(4) 


Zj    sj     L    L    sj 
j=i        j=i  i=i 


n 


n   m 


(5) 


-IFjt-  1  Zfit 

j=l      j=l  i=l 


=  -v. 


The  constraints  (Ij),  for  1  <  j  <  n,  (2)  and  (3)  in  problem  (3.6.3) 
are  the  dual  constraints  obtained  from  problem  (3.6.2).   The  additional 
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constraints  (4)  and  (5)  are  appended  to  the  original  constraints  (lj) , 
(2)  and  (3) ,  and  the  new  variable  v  is  introduced  without  changing  the 
problem.   The  additional  constraints  are  obtained  from  the  sum  of  the 
constraints  (lj),  for  1  <  j  <  n,  so  that  they  are  redundant  in  problem 
(3.6.3). 

TABLE  3 . 1 
DUAL  VARIABLES 


Corresponding 
Variable  Name  Range  of  Indices         Constraint 


F.,  1  <  j  <  k  <  n  (la) 

Jk 

F,  .  1  <  j  <  k  <  n  (lb) 

kj 

F  1  <  j  <  n  (2a) 

J  * 

F  1  <  j  <  n  (2b) 

sj 

f  1  <  j  <  k  <  n  (3a) 

jk 

f  1  <  j  <  k  <  n  (3b) 

kj 

f1  1  <  j  <  n,   1  <  i  <  m         (4a) 

f1.  l<j<n,   l<i<m         (4b) 

sj 


The  coefficient  matrix  of  constraints  (lj),  1  <  j  <  n,  (4)  and 
(5)  is  a  node-arc  incidence  matrix  of  a  network,  denoted  by  N,  with 
n +  2  nodes.   The  n  constraints  (lj),  for  1  <  j  <  n,  correspond  to 
a  complete  subnetwork  of  n  nodes  j ,  called  intermediate  nodes,  such 
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that  between  each  pair  of  nodes  j  and  k  there  are  two  pairs  of  directed 
arcs.   For  any  two  intermediate  nodes,  arcs  in  a  pair  have  the  same 
direction  but  have  opposite  direction  to  the  other  pair  of  arcs;  the 
arcs  in  a  pair  have  costs  -c'   and  zero,  respectively. 

Node  n+1  corresponds  to  constraint  (4)  and  is  denoted  as  the 
source,  s.   Node  n+2  corresponds  to  constraint  (5)  and  is  denoted  as 
the  sink,  t.   The  source  and  sink  are  connected  to  the  other  n  inter- 
mediate nodes  in  the  following  manner.   For  each  intermediate  node  j  of 
the  complete  network  there  are  m+1  directed  arcs  from  the  source  s  to 

node  j.   Each  of  the  first  m  arcs  has  cost  -a'.  ,  1  <  i  <  m,  and  the 

st 
(m+1)  '  arc  has  cost  -D'. .   There  are  also  m+1  directed  arcs  from  the 

J 

node  j    to  the  sink  t.   Each  of  the  first  m  arcs  has  cost  a! ,  1  <  i  <  m, 

st 

and  the  (m+1)    arc  has  cost  d'. . 

J 

The  variables  in  problem  (3.6.3)  will  be  interpreted  as  flow 
variables  on  the  corresponding  arcs  in  the  network  N.   The  objective 
of  (3.6.3)  is  to  determine  the  flow  of  maximum  total  cost  in  the  net- 
work N. 

Constraint  (2)  of  problem  (3.6.3)  implies  that  there  are  weights 

associated  with  some  arcs  in  N  and  that  the  sum  of  the  weighted  flow  in 

these  arcs  must  equal  one.   The  arcs  that  do  not  have  weights  associated 

with  them  are  the  arcs  from  node  s  to  node  j  with  cost  -D'.  and  the  arcs 

from  node  j  to  node  t  with  cost  d'.  for  1  <  j  <  n.   Also  the  arcs  between 

J 

pairs  of  intermediate  nodes  j  and  k  with  cost  -c1   for  all  1  <  j  <  k  <  n, 

jk 

do  not  have  weights.   All  such  arcs  will  be  called  unweighted  arcs. 
The  remaining  arcs  will  be  called  weighted  arcs. 
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Then  problem  (3.6.3)  is  to  determine  the  flow  in  N  of  maximum 
total  cost  subject  to  the  restriction  that  the  total  of  the  weighted 
flow  in  the  weighted  arcs  must  equal  one. 

For  n=3  and  m  -  2 ,  the  network  N  is  illustrated  in  Figure  3.1. 
The  number  corresponding  to  each  arc  is  the  arc  cost.   Broken  lines 
represent  unweighted  arcs,  and  unbroken  lines  represent  weighted  arcs. 

The  unweighted  arcs  in  the  network  N  correspond  to  constraints 
in  the  linear  program   (3.6.2)  that  are  determined  from  the  constraints 
of  problem  (3.6.1).   Therefore,  for  the  problem  (3.2.3),  where  not  all 
constraints  were  assumed  to  exist,  unweighted  arcs  in  the  network  N 
would  only  exist  for  these  constraints  that  are  included  in  (3.2.3). 
The  network  N  associated  with  the  unconstrained  problem  (3.2.2)  would 
have  no  unweighted  arcs. 

Similarly,  weighted  arcs  in  the  network  N"  correspond  to  con- 
straints in  the  linear  program  (3.6.2)  that  are  determined  from  the 
objective  function  of  problem  (3.6.1).   For  the  problem  (3.2.3), 
weighted  arcs  in  the  network  N  would  exist  only  for  those  terms  in  the 
objective  function  with  positive  v.,  or  w... 

For  the  one  facility  problem,  obtained  by  letting  n=l,  the 
network  N  has  only  one  intermediate  node,  but  still  has  a  source  and 
sink  and  connecting  arcs  as  in  the  discussion  above. 

3. 7   Properties  of  the  Network  Problem 

After  establishing  some  definitions,  properties  of  solutions 
to  problem  (3.6.3)  are  developed.   Define  a  chain  [ST   from  node  i 
to  node  i   in  a  network  to  be  a  set  of  directed  arcs 

q 
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Figure   3.1      A  Network  N 
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{(i  ,i  ),(i,  ,i0), (i,  ,i  )}  where  (i.,i.  n)  is  a  directed  arc  from 

pi    12        kqJ         J   J+l 

node  i   to  node  i    in  the  network.   A  cycle  is  a  chain  with  i  =i  , 
j  J  +  1  P   q 

and  a  simple  chain  is  a  chain  that  contains  no  cycles,  that  is,  the 

i   are  distinct.   For  any  chain  C=  f (s , j  ) , ( j  , j0) (j„>t)]  inN 

J  l    l  ^        p 

from  s  to  t,  denote  the  total  of  the  arc  costs  in  C  by  T  and  denote  the 
total  of  the  arc  weights  in  C  by  W.   Also,  let  F  represent  the  flow  in 
each  arc  of  the  chain  C. 


Lemma  3.7.1.   There  exists  a  positive  feasible  flow  in  N  with  positive 
objective  function  value. 

Proof:   Of  all  the  costs  -a' ,  1  <  i  <  m,  for  arcs  between  node  s  and 
node  1,  let  -a'  be  the  largest.   Of  all  the  costs  a',  for  arcs  between 

1  and  node  t,  let  a'  be  the  largest.   Consider  the  chain  C =  { (s , 1) , (1 , t) } 

m 

with  arc  costs  -a'  and  a' ,  respectively.   Now  a'  *  a'   for  otherwise 
1      m  in    j. 

all  the  a'  would  be  equal  and  problem  (3.2.1)  would  have  the  trivial 
i 

solution  x'  =  a*   for  all  1  <  j  <  n.   Therefore,  T  =  a'  -  a'  >  0  and 
j    i  mi 

1    m   1  . 
W=  w   +w   .   A  feasible  solution  is  obtained  by  taking   fe1  =  1n-  =  w 
11   lm  SJ.J.LH 

and  taking  all  other  flow  variables  zero.   This  solution  has  objective 
function  value  T/W  >  0. 

Lemma  3.7.2.   Any  cycle  in  N  has  nonpositive  total  arc  cost  and  consists 
only  of  arcs  between  intermediate  nodes. 

Proof:   All  coefficients  in  constraint  (4)  of  problem  (3.6.3)  are  posi- 
tive so  that  there  are  no  arcs  directed  into  node  s.   Similarly,  all 
coefficients  in  constraint  (5)  are  negative  so  that  there  are  no  arcs 
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directed  out  of  the  sink  t.   Hence  any  cycle  must  occur  among  inter- 
mediate nodes.   From  the  objective  function,  all  arc  costs  for  arcs 
among  intermediate  nodes  are  nonpositive.   Thus  any  cycle  will  have 
nonpositive  total  arc  cost. 

The  next  lemma  characterizes  a  feasible  flow  in  N  for  use  in 
the  proof  of  subsequent  properties. 

Lemma  3.7.3.   Any  feasible  flow  in  N  consists  of  either  a  cycle  of  flow, 
a  chain  of  flow  from  s  to  t  of  unweighted  arcs,  a  chain  of  flow  from  s 
to  t  of  at  least  one  weighted  arc,  or  some  combination  of  these. 

If  there  is  an  unweighted  chain  in  N  of  zero  total  arc  cost, 
then  any  feasible  solution  to  problem  (3.6.3)  can  include  any  amount 
of  flow  on  this  chain  without  affecting  the  objective  function  value. 
To  exclude  from  further  consideration  those  feasible  solutions  that 
contain  flow  on  such  chains  allows  a  more  explicit  presentation  and 
causes  no  loss  of  generality.   Therefore,  it  is  assumed  for  the  remainder 
of  the  development  that  any  feasible  solution  has  zero  flow  on  unweighted 
chains  of  zero  total  arc  cost. 

Property  3.7.1.   Problem  (3.6.1)  has  no  feasible  solution  if  and  only 
if  there  exists  a  chain  of  unweighted  arcs  in  N  which  has  a  positive 
total  arc  cost. 


Proof:   Suppose  there  exists  a  chain  C  in  N  of  unweighted  arcs  with 

total  arc  cost  T,  >  0.   By  Lemma  3.7.1  there  exists  a  feasible  flow  on 
k 

some  chain  C.  with  T.  >  0  and  W.  >  0.   Then  putting  a  flow  of  1/W.  on 
11  l  ^  l 

each  arc  of  C.  and  putting  infinite  positive  flow  on  C  is  a  feasible 


50 


solution  to  problem  (3.6.3)  with  unbounded  objective  function  value. 
Since  problem  (3.6.2)  is  the  linear  programming  dual  of  problem  (3.6.3), 
problem  (3.6.2)  has  no  feasible  solution. 

A  contrapositive  proof  is  given  for  the  converse  statement. 
Suppose  all  chains  from  s  to  t  of  unweighted  arcs  in  N  have  nonpositive 
total  costs.   Since  problem  (3.6.3)  is  a  maximization  problem,  Lemma 
3.7.1  and  the  above  assumption  imply  that  any  optimal  solution  will 
have  zero  flow  on  all  such  chains.   There  will  also  be  zero  flow  on  all 
cycles,  since  cycles  have  nonpositive  total  arc  cost  by  Lemma  3.7.2. 
By  Lemma  3.7.3  all  remaining  feasible  solutions  will  correspond  to 
chains  from  s  to  t  with  at  least  one  weighted  arc  in  each  chain.   Then 
constraint  (2)  implies  that  flow  on  these  chains  is  bounded.   Thus  the 
existence  of  a  bounded  feasible  solution  to  problem  (3.6.3)  implies 
that  there  is  an  optimal  solution  to  problem  (3.6.3).   Hence  there 
exists  a  feasible  solution  to  problem  (3.6.2). 

Property  3.7.2.   If  a  feasible  solution  to  problem  (3.6.3)  is  optimal 
(with  bounded  objective  function  value),  then  the  optimal  solution 
consists  of  positive  flow  on  simple  chains  from  s  to  t  with  positive 
total  arc  weight. 

Proof:   An  optimal  solution  to  problem  (3.6.3)  with  bounded  objective 
function  value  implies  that  all  chains  of  unweighted  arcs  in  N  have 
nonpositive  total  arc  cost.   If  such  a  chain  has  zero  total  arc  cost, 
by  assumption   there  is  zero  flow  on  that  chain.   If  such  a  chain  has 
negative  total  arc  cost,  a  positive  flow  on  the  chain  will  not  be 
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optimal.   Thus,  flow  on  unweighted  chains  of  nonpositive  total  arc 
cost  is  zero  in  the  optimal  solution. 

Also,  all  cycles  have  nonpositive  total  cost  by  Lemma  3.7.2 
and  will  have  zero  flow  in  an  optimal  solution. 

By  Lemma  3.7.1  an  optimal  solution  will  have  positive  flow  on 
some  arcs.   Thus  by  Lemma  3.7.3  an  optimal  solution  must  correspond  to 
positive  flow  on  chains  from  s  to  t,  each  of  which  contains  at  least  one 
weighted  arc.   These  chains  do  not  contain  cycles  since  all  cycles  have 
zero  flow.   Thus  an  optimal  solution  consists  of  positive  flow  on 
simple  chains  from  s  to  t  that  contain  at  least  one  weighted  arc. 

Property  3.7.3.   Suppose  problem  (3.6.1)  has  a  feasible  solution  and 

let  C, ,C  ,  ...,C  be  all  the  simple  chains  from  s  to  t  in  N  with  chain  C. , 
1  2,  r  l 

having  total  arc  cost  T  and  total  arc  weight  W  >  0.   Let  chain  C 

l  l  k 


have  the  property  that 


T  /W,  £  T./W. 


for  1  <  i  <  r.   Then  putting  a  flow  of  value  F  =  1/W  on  all  arcs  in 

k      k 

C  and  zero  flow  elsewhere  is  an  optimal  solution  to  problem  (3.6.3). 

Proof:  Since  problem  (3.6.1)  has  a  feasible  solution,  problem  (3.6.3) 
has  a  bounded  optimal  solution  and  ay  Property  3.7.2  any  optimal  solu- 
tion will  be  determined  by  simple  chains  with  positive  total  arc  weights. 

Suppose  that  positive  flow  on  the  simple  chains  Ci ,C-j   is 

Jl      Jp 

an  optimal  solution  to  problem  (3.6.3)  with  objective  function  value 
and  total  weighted  flow  given  by 


p 


I  FJi  TJi  and  I  FJi  wJi  =  1. 


i=l  i=l 

respectively.       But 

P  P 


I   FJi  TJi  I   FJi  wJi  ^T 


i=l  i=l 

P 


T 
V 

k 


~       w~     L    FJi  wJi 


i=l 

Tk 

—  -FT 
W,     ~    Vk    ' 
k 


Thus    the   flow   F     on   chain  C      is   also   an  optimal    solution  to 
k  k 

problem    (3.6.3). 

Property  3.7.3  provides  a  basis  for  a  solution  procedure  to 

problem  (3.6.3).   That  is,  if  problem  (3.6.3)  has  a  bounded  objective 

function  value,  an  optimal  solution  is  determined  by  the  simple  chain, 

C,  ,  from  s  to  t  with  T,  >  0  and  maximum  ratio  T,  W,  ,  referred  to  as  the 
k  k  k   k 

cost  to  weight  ratio.   A  procedure  for  determining   a   chain  in  N  of 
maximum  cost  to  weight  ratio  is  developed  next.   This  procedure  also 
determines  if  there  is  an  unbounded  solution  to  problem  (3.6.3),  that 
is,  if  there  is  a  chain  of  unweighted  arcs  with  positive  total  arc 
cost. 

Similar  optimization  problems  concerning  cycles  in  general 
networks  have  received  previous  attention  in  the  literature  and  may 
be  considered  as  special  cases  of  hyperbolic  programming  problems. 
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Dantzig,  Blattner  and  Rao  [5]  consider  a  ship  routing  problem  that  is 
formulated  as  a  problem  of  determining  a  cycle  in  a  network  of  minimum 
cost  to  time  ratio.   Their  solution  procedure  is  a  variation  of  the 
simplex  algorithm.   Fox  [9]  considers  the  same  problem  and  solves  it 
using  a  parametric  application  of  the  out-of -kilter  algorithm. 

A  problem  equivalent  to  problem  (3.6.3)  may  be  formulated  and 
interpreted  as  finding  a  cycle  in  a  corresponding  network  of  maximum 
cost  to  weight  ratio  from  those  cycles  with  positive  total  weight. 
Using  such  a  formulation,  some  developments  in  [5]  and  [9]  would  apply. 
However,  an  exposition  in  terms  of  chains  is  chosen  because  of  the 
presence  of  the  unweighted  arcs  and  the  special  structure  of  the  net- 
work N. 

Notice  that  for  the  unconstrained  problem  (3.6.2),  there  are 
no  unweighted  arcs  in  the  network  N  and  the  dual  problem  is  always 
bounded. 

Consider  a  network  N'  with  the  same  node-arc  structure  as  N 
but  with  different  arc  costs.   The  arc  costs  in  N'  are  obtained  from 
the  arc  costs  of  N  by  the  following  procedure.   Subtract  X/w . .  from 
the  corresponding  arc  costs  a!  and  -a',  for  1  <  i  <  m  and  1  <  j  <  n. 
Subtract  X/v   from  the  zero  arc  cost  of  corresponding  intermediate 
arcs.   Then  multiply  all  arc  costs  by  -1.   The  arc  costs  for  arcs  in 
N'  are  defined  as  follows: 
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r a\    +    X/w, 


Ji 


D'. 
J 


c   +  \W   =  / 

pq    pq  \ 


-a'  .  +  x/w . . 

-d' 
J 


X/v„ 


V 


"jk 


P=s,  q  =  j,  1  <  i  <  m" 


1  <  j  <  n 


P  =  s,  q  =  j  ,  l 


m  +  1 , 


P  =  j,  q  =  t,  1  <  i  <  m" 


P  =  j,q  =  t,  i  =  m  +  1 


1  <  j  <  n 


The  network  N'  corresponding  to  the  network  N  in  Figure  3.1  is 
given  in  Figure  3.2. 

Every  simple  chain  from  s  to  t  in  X  corresponds  to  a  unique 
simple  chain  from  s  to  t  in  N"  having  the  same  node-arc  progression 
but  different  arc  costs.   Denote  a  simple  chain  in  N  by  C.  and  the  corre- 
sponding simple  chain  in  N1  by  C .   Any  chain  C'  in  N'  has  total  arc 
cost  denoted  by  T'  and  defined  by 


T'. 

l 


I 


(c   +  \W   ) 

pq    pq 


(P.q)ec: 


If  the  corresponding  chain  C.  has  positive  total  arc  weight, 
the  total  arc  cost  of  C  is  a  strictly  increasing  function  of  \. 
For  each  chain  C'  whose  corresponding  chain  C.  has  positive  total  arc 
weight,  define  the  root  value  of  C' ,  denoted  X. ,  as  the  value  of  X  for 

which  chain  C  has  zero  total  arc  cost.   That  is, 

l 


I 


(c   +  X.w  )  =  0  . 

pq   i  pq 


(p.q)ec^ 
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a  +\/w 


ai+X/wil 


Figure  3.2      A  Network   N' 
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Solving  for  X.  gives 


\.  = 


^       pq     f  L  pq    W. 

(p,q)6C^        /  (p,q)CC^         * 


by  the  definition  of  c   and  W   in  terms  of  the  chain  C  in  N. 

pq     pq  i 

Therefore,  finding  the  chain  in  N  with  positive  total  arc  weight 
and  maximum  cost  to  weight  ratio  is  equivalent  to  finding  the  correspond- 
ing chain  in  N"  having  the  largest  root  value.  Denote  the  largest  root 
value  by  X*. 

Property  3.7.4.   If  problem  (3.6.3)  has  an  optimal  solution  with  bounded 
objective  function  value,  say  z' ,  then  X*  =  z' 

Proof:   Let  chain  C  correspond  to  X*,  that  is, 
o 


y  -c    /  y   w 

L  pq       L  p 


X*  = 

pq  /     La  pq 

(p.q)eC'  /  (p,q)eC 

o  '        o 


The  corresponding  chain  C   in  N  has  total  arc  weight  W  and  total  arc  cost 

o  o 

T  with 
o 

T    T 

X*  =  -2.  >  -L 

WW 

O      1 


for  all  chains  C.  in  X.   By  Property  3.7.3,  the  optimal  flow  on  C  is 

F  =  1/W  .   Then  the  optimal  objective  function  value  is  F  T  = T  /W  =  X*. 
o      o  o  o   o  o 
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Consideration  is  now  given  to  chains  in  N  of  unweighted  arcs. 

If  a  chain  C  in  N'  corresponds  to  a  chain  C.  in  N  of  unweighted  arcs, 
i  i 


then  C  has  total  arc  cost 
i 


X 


c 

pq 


(P.q)ec^ 


which  is  a  constant  with  respect  to  X,  since  all  W   are  zero  in  C . 

If  the  constant  total  arc  cost  of  C  is  negative,  then  C.  in  N  is 

a  .chain  of  unweighted  arcs  with  positive  total  arc  cost.   By  Property  3. 7.  1, 

problem  (3.6.1)  has  no  feasible  solution. 

3. 8  A  Solution  Procedure 

Using  the  properties  of  the  previous  section,  a  solution  proce- 
dure is  developed  for  problem  (3.6.3).   The  following  algorithm  deter- 
mines a  chain  C  in  N' ,  having  the  largest  root  value  of  all  chains  in 
N' ,  or  it  determines  if  there  is  a  chain  in  N'  of  constant  total  arc 
cost  which  is  negative.   Therefore  the  algorithm  either  determines  the 
chain  in  N  corresponding  to  the  optimal  solution,  or  it  determines  that 
there  is  no  bounded  optimal  solution. 

Algorithm 

1.  Set  i  =  0,  X.  =  0. 

l 

2.  Using  the  current  value  of  \. ,  find  the  chain  from  s  to  t 
in  N'  of  minimum  total  arc  cost.   This  may  be  done  by  using  any  short- 
est chain  algorithm  that  handles  both  positive  and  negative  arc  costs, 
for  example   the  algorithm  given  in  Ford  and  Fulkerson  [8].   Use  the 
arc  costs  determined  by  the  current  value  of  X.  as  the  arc  lengths  in 
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the  shortest  chain  algorithm.   Denote  the  chain  of  minimum  total  arc 

cost  by  C!  „  and  its  total  arc  costs  by  T'  ..   If  T'    =0,  stop: 
l+l  l+l       i+1 

C  is  optimal  and  X.  is  the  largest  root  value.   If  T'  <  0  and 
11  i 

E      W   >  0,  go  to  step  3.    If  T!  <  0  and     I      W   =  0, 

(P,q)ec:+1    pq  x  <P'^eci+i    Pq 

stop;  problem  (3.6.1)  is  infeasible. 

3.   Set   T!  ,  =  0,  compute  X.  ,  by 
l+l  l+l 


£ 


i  +  1        L,  pq  /     L  pq 

(P,q)ec;+ 

Set  i  +  1  =  i ,  go  to  step  2. 

The  following  remarks  show  that  the  algorithm  converges  finitely 
to  the  desired  chain. 

Remark  3.8.1.   X.  ,  >  \  . 

l  +  l     i 

Proof:   If   i  =  0,  then  \  =  0,  and  by  Lemma  3.7.1  there  exists  some 

o 

chain  in  N  with  positive  total  arc  cost  and  positive  total  arc  weight. 
Then  the  corresponding  chain  in  N"  has  negative  total  cost.   There- 
fore, with  X  =  0,  the  minimum  arc  cost  chain,  denoted  by  C' ,  has 
o  l 


T^  <  0. 


Next,  suppose  that  i  £  0,  and,  using  the  current  value  of  X.  , 


i 


chain  C  „  has  T*    <  0  and   .    £  „,    W   >  0.   Then  by  the  defin- 
i+1       l+l  (p,q)eC.     pq 

ition  of  Tl  .,  and  X.  .. 
l+l      l+l 


'..  <     V     -c    /    )       W 
i       L>  pq  /     L,  ] 

(P,q)eci      /  (p.q)eC' 

1+1       I  1-rl 


pq     l-  1 
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and  the  assertion  holds.   If  T'    <  0  and     *■      W   =  0, 

i+l         (P.q)£Ci+1   pq 

then  by  step  2  the  algorithm  stops  without  specifying  X.  . 

If  T'    =0,  then  the  algorithm  stops  at  optimality  without  specify- 
i+1 

ing  X.   . 
i+l 

Remark  3.8.2.   The  algorithm  is  finite. 


Proof:   At  each  iteration  the  algorithm  either  stops  or  a  new  root 

value  X    >  X  is  determined.   Since  the  number  of  chains  is  finite, 
i+l     i 

and  any  chain  is  considered  at  most  once,  the  largest  X.  will  be  deter- 
mined in  a  finite  number  of  iterations. 

Remark  3.8.3.  The  chain  identified  at  the  termination  of  the  algorithm 
corresponds  to  a  chain  in  N  that  is  either  an  unweighted  chain  of  posi- 
tive total  arc  cost  or  a  weighted  chain  of  maximum  cost  to  weight  ratio. 

Proof:   Suppose  the  algorithm  stops  at  step  2  by  determining  a  chain 

C  .  with  T'  .  -   0.   Then  i  ^  1  and  the  previous  chain  C!  has,  for 
i+l       i+l  i 

the  current  root  value  X. ,  total  arc  cost  T!  -   0  and     £  .  w   >  0. 

i  i         (p,q)eC^  pq 

Since  T!  is  the  minimum  arc  cost,  T!  <  T'.  for  all  other  chains  C.  in 
i  i    J  J 

N'  ,  using  the  value  X..   Therefore,  in  the  network  N,  T./W.  2:  T./W.  for 

all  chains  C.  in  N,  and  C.  is  a  weighted  chain  of  maximum  cost  to  weight 
J  i 

ratio. 

Suppose  the  algorithm  stops  at  step  2  because  the  chain  C! 

has  T'.  <  0  but  „   £      w   =  0.   Then  the  chain  C    in  N  is  an 
i+l         (p,q)eC^+1   pq  .i+l 

unweighted  chain  of  positive  total  arc  cost. 
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If  problem  (3.6.3)  has  an  optimal  solution  with  bounded  objec- 
tive function  value,  then  the  algorithm  identifies  a  chain  C.  in  N 
that  determines  an  optimal  solution  to  (3.6.3).   However,  it  is  desired 
to  obtain  all  optimal  solutions  to  the  primal  problem  (3.6.1).   If  a 
labeling  algorithm  is  used  in  step  2  to  find  the  minimum  cost  chain, 
one  optimal  solution  to  problem  (3.6.1)  can  be  determined  from  the  node 
labels  (dual  variables).   The  node  labels  also  indicate  a  procedure  to 
determine  alternative  solutions. 

The  dual  of  the  minimum  cost  chain  problem  for  the  network  N' , 
at  iteration  i,  is  written  below  with  the  following  change  of  variables; 
each  variable  is  multiplied  by  -1. 


min   s 

-  t 

s.  t. 

x' 

j 

-    X' 

k 

-x' 

J 

♦* 

x' 

J 

-x' 

J 

x' 

3 

-*i 

-x' 

j 

+    x' 
k 

x'. 

J 

-x'. 

J 

-  t 


+  s 


£ 

-J* 

^ 

"S 

2: 

dJ 

£ 

-D' 
J 

£ 

-X./v.. 

£ 

-X.    v 
i      Jk 

1  <  j  <  k  <  n 


1  <  j  <  n 

(3.8.1) 
1  <  j  <  k  <  n 


-  t  £   a'.  -  X./w.. 
i    i   Ji 

+  s       s  -a!  -  X.  /w  . . 
i   i  ji 


1  <  j  <  n 


At  the  last  iteration,  X.  =  X*  and  by  Property  3.7.4,  X*=z!  , 
the  optimal  objective  function  of  problem  (3.6.1).   Since  the  source  s 
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has  no  arcs  directed  into  it,  the  minimum  cost  chain  algorithm  gives 
s  =  0.   Since  the  minimum  cost  chain  at  the  terminal  iteration  has 
zero  cost,  t  =  0.   This  implies  that,  at  optimality,  the  constraints 
to  (3.8.1)  are  equivalent  to  the  constraints  of  problem  (3.6.2). 
Therefore  the  node  labels  in  appropriate  shortest  chain  algorithms 
determine  the  primal  variables  x'. ,  1  <  j  <  n,  that  form  an  optimal 
solution  to  problem  (3.6.1). 

If  the  minimum  cost  chain  does  not  pass  through  some  node, 
say  j ,  then  flow  into  and  out  of  node  j  is  zero.   This  corresponds  to 
a  degenerate  solution  to  problem  (3.6.3)  so  that  its  dual,  problem 
(3.6.2),  has  alternative  solutions.   The  alternative  solution  values 
for  x'.  may  be  determined  using  those  constraints  in  (3.8.1)  that 
involve  x'.  and  using  the  value  z'  =  X*.   Then  from  the  development  in 
Section  2.7,  the  alternative  values  for  x'.  must  satisfy 


,         max        max   ,  ,     ,  .    ,H 
x'.  2:  max        d*.  ,  ,-_.   <a   -  z'/w.l 
j       l_3<J^n  j    l<j<n  L  i       jiJJ 


l<i<m. 


-       mm   _,     min   -  ,     ,  ,    ,~j 

x*.  <  min   D'.  ,  ,__   |a!  +  z'/w..} 

j    l<j<n  -  i     1   jiJJ 


(3.8.2) 


mm 
'  L_l<X<n  "j  '  l<j<n 
l<j<m 


and 


\x'    -   x'  |  <  min  fc*  ,  z'/v   }   1  <  k  <  n,  j  4   k.        (3.8.3) 

J       K  JK     X    JK 

The  expressions  (3.8.2)  may  be  used  to  determine  an  interval 
for  x'. .   However,  expression  (3.8.3)  is  an  additional  constraint  on 

<J 

alternative  solutions  that  determines  subintervals  depending  on  other 
variables  x"  . 
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Using  the  optimal  solutions  to  the  subproblems  (3.5.2)  and 

(3.5.3),  a  procedure  is  given  based  on  Section  2.7  for  finding  alter- 
native solutions  to  the  original  multif acility  problem  (3.5.1). 

r  O*        O*  ■>      j   r  °'        o'  , 
Solutions  (x   x   J   and   (y   ,...,y   J  to  subproblems 

(3.5.2)  and  (3.5.3),  respectively,  are  obtained  using  the  algorithm. 

Let  z'  and  z'  be  the  respective  optimal  objective  function  value  of  the 

J.         £ 

subproblems.   Then  Property  3.5.2  gives  z1  =  max  (z-l.z')  as  the  optimal 
objective  function  value  of  (3.5.1),  and  gives  a  procedure  for  deter- 
mining one  optimal  solution  to  (3.5.1).   As  in  Section  2.7,  there  are 
three  cases  that  determine  all  alternative  solutions. 

Case  1.   z'  =  z'  =  z'    Then  all  the  solutions  fx'...,x* ]  to 
o    1    2  L  1      nJ 


em 


problem  (3.5.2)  paired  with  all  the  solutions  [yl,...,y'}    to  probl 

(3.5.3)  give  all  the  optimal  solutions  ( (x'y' ),..., (x' ,y' ) }  to 

problem  (3.5.1). 

Case  2.   z'  =  z'  >  z' .   Then  for  all  the  points  y'.  for  1  <  j  <  n, 
o    1    2  J 

that  satisfy  the  expressions 


y'.  £  max  ,_,_   d'.  ,  ,-_.   [b   -  z'/w..j 
1  j       Ll^J^n  J    l^J<n  l  l    o  jiJJ 


l<i<ni 


(3.8.4) 


_       min   ,     min   r.  .     .     ,1 
y'  <  min  ^^,     D'.  ,  ,  . .  ..   fb  +  z*/w. . 1 


ly'  -  y.'l  <  min  [cl,  ,  z'/v.,  ]   1  <  k  <  n,   j  *  k        (3.8.5) 
|Jj   ■'k1       l  jk'   o  jkJ 

the  points  (x°  ,y'.),  1<  j  <  n,  ar  :  alternative  solutions  to  (3.5.1), 

for  all  alternative  solutions  rx,  ,...,x  ]  to  problem  (3.5.3). 

'1       n  J 
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Case  3.  z."0  =   z'  >  z'    Then  for  all  the  points  x'. ,  1  <  j  <  n,  that 
satisfy  expressions  (3.8.2)  and  (3.8.3),  with  z'  replacing  z ' ;  the 
points  (x*.  ,y.),  1  <  j  <  n,  are  alternative  solutions  to  (3.5.1),  for  all 
alternative  solutions  fy   ,...,y   }  to  problem  (3.5.2). 

By  Property  3.5.1,  optimal  solutions  to  the  problem  (3.2.1)  are 
obtained  by  taking  the  inverse  rotation,  T   ,  of  the  solutions  to 
(3.5.1). 

3. 9  A  Special  Case 

A  special  case  of  problem  (3.2.3)  is  considered  in  this  section. 

It  is  obtained  by  assuming  that  none  of  the  new  facilities  are  linked; 

that  is,  all  v.,  are  zero  and  there  are  no  constraints  for  the  new 
Jk 

facilities.   The  problem  is 

111111       m3X   iw   II  X   -  P  II  \ 
X     ,X  €E_  3<j<n  IV1  Xj    Pi"lJ 

1      n  2  ieW(J) 

(3.9.1) 
s.t.     ||  X   -  PiJ|1  <  d     1  <  j  <  n,  i  C  D(j). 


It  is  assumed  that  each  new  facility  X.  is  linked  to  some  exist- 

J 

ing  facility  P.;  otherwise,  X.  could  be  located  independently  of  P.. 
It  is  also  assumed  that  each  existing  facility  P.  is  linked  to  some 
new  facility  X.;  otherwise,  P.  could  be  deleted  from  the  problem. 
These  assumptions  exclude  from  consideration  problems  that  are  not  well 
formulated. 
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Then  problem  (3.9.1)  may  be  treated  as  n  independent  one  facility 
subproblems: 

min     max 


VE2 


uw«>lwjil,xJ  -pi"J 


(3.9.2) 
s.t.    ||  X.  -  P.  |!   <  d..      i  e  D(j) 

for  1  <  j  <  n.   Each  of  these  problems  is  solved  using  the  development 
of  Chapter  2,  and  from  these  solutions,  all  alternative  solutions  are 
obtained. 

The  transformation  T  is  applied  to  each  of  the  problems  (3.9.2) 
and  its  subproblems  are  defined  for  1  <  j  <  n.   Let  z! .  and  z^ .  be  the 

J  "J 

optimal  objective  function  values  for  the  j    subproblems  in  x'.  and 

y*  ,  respectively.   Then  z  =  V^  ,__   fz'  .  ,  z '  }   is  the  optimal  objec- 
j  o      l<j<n  l  lj   2jJ 

tive  function  value  of  problem  (3.9.1).   All  alternative  solutions  can 
be  determined  from  the  following  procedure. 

th 
Points  X   that  are  alternative  solutions  to  the  j    subproblem 
J 

(3.9.2)   are  those  points  satisfying 

w..  ||  X.  -  P.  ||,  <  z       i  €  W(j) 

II  Xj  -  P±||1  <  dj.     i  e  D(j)  . 

Based  on  the  development  in  Section  2.8,  define 

max    ,     .  ,    max  ,   >1 

e,  (z  ,j)  =  max    , ,.  ..  (a.  +  b.  -  z   w..),  ,„_,..  (a.  +  b.  -  d   ) 

1  o  J       Lie^(j)   i    i    o  ji   ieD(j)   i    i    ji  J 

max    ,  i.  \    max    /      V,     A      \    I 

e„(z  ,j)  =  max      ._,.  (a.  -b   -z   w..),  .       (a.  -b.  -d..) 

2  o  J        Lie^'(j)    i     i    o  ji    ieD(j)    i     i    Ji  J 
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e„(z  ,j)  =  max 
J   o 


max  ,     .      .    .    max  1 

,r/.x  (-a.  +b.  -z  /w..),  .„_,.,.  (-a.  +b.  -  d .  .  ) 

LieW(j)    i   1   o  ji  ieD(j)    i   i   ji  J 

max  ,     ,      .    «    max  x  ] 

e  (z  ,j)  =  max    ,  ,  ..  (-a .  -  b.  -  z  /w  .  . )  ,  .  (-a.-b.  -d..) 

4   o'J        LieW(J)     l    l    o   ji  i€D(j)     l    l    ji  J 


Then  the  four  inequalities 


Xj  +  yj  S  el(V  j) 

Xj  "  yj  *  e2(Zo'  J) 

"Xj  +  yj  *  G3(Zo'  J) 

~Xj  "  yj  -  e4(V  j) 


for  1  <  j  <  n,  determine  all  alternative  optimal  solutions  to  problem 
(3.9.1). 

3. 10   An  Extended  Problem 

This  section  contains  a  discussion  of  the  extended  problem, 
(3.2.4).   The  development  of  Sections  3.5  and  3.6  apply  to  problem 
(3.2.4)  without  any  additional  assumptions  on  the  g..  and  h   . 
That  is,  two  subproblems,  corresponding  to  (3.5.2)  and  (3.2.3),  may- 
be defined  and  each  subproblem  has  a  dual  with  a  corresponding  network. 

Let  N"  denote  the  network  of  the  subproblem  in  x'  , . .  .  ,  xr  .   The  network 

In 

N"  has  the  same  node-arc  structure  as  N,  but  with  different  arc  cost. 

Arcs  in  N  incident  to  node  j  with  costs  -a'  and  a',  have  corresponding 

11 

arc  costs  -a'  +  g'  /w    and  a'  +  g'  /w  ,  respectively,  in  N"  for 
i    Ji   Ji       i     Ji   i 

1  <  j  <  n,  1  <  i  <  m.   An  arc  with  a  cost  of  zero  in  N  has  a  correspond- 
ing arc  cost  h   ,  v    in  N"  for  1  <  j  <  k  <  n.   Unweighted  arcs  in  N"  have 
Jk   jk 

the  same  arc  costs  as  corresponding  unweighted  arcs  in  N. 
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Assumptions  must  be  made  on  the  g.,  and  h.   in  order  for  the 

J  -      Jk 

properties  and  lemmas  of  Section  3.7  to  hold  for  the  network  X"  and 

its  corresponding  subproblem. 

If  there  are  positive  h   ,  then  the  weighted  arcs  between 

intermediate  nodes  will  have  positive  cost.   Therefore  the  network 

used  in  the  algorithm  may  have  cycles  of  negative  total  cost  and  the 

shortest  chain  algorithm  cannot  be  used.   Hence  the  above  procedure 

will  not  apply  if  there  exist  positive  h    such  that  a  cycle  in  N"  has 

positive  total  arc  cost.   If  negative  constants  h    exist   such  that 

Jk 

a  cycle  has  negative  total  cost,  there  will  always  be  zero  flow  on  that 
cycle  at  optimal ity.  The  corresponding  arcs  between  intermediate  nodes 
will  not  be  in  the  optimal  solution. 

If  the  g..  are  nonpositive,  then  Lemma  3.7.1  may  not  hold. 
Nevertheless,  the  results  in  [5]  appear  to  extend  so  that  the  optimal 
solution  to  the  network  problem  still  corresponds  to  simple  chains 
in  N".   Additional  research  is  required  to  completely  extend  the  proce- 
dure to  nonpositive  g...   If  assumptions  are  made  on  the  constants  h 

J1  jk 

and  S-.     to  avoid  both  of  the  above  situations,  then  the  properties  and 
solution  procedures  developed  for  problem  (3.2.1)  also  hold  for  problem 
(3.2.4). 


CHAPTER  4 
A  p-CENTER  LOCATION  PROBLEM 

4. 1   Introduction 

Considered  in  this  chapter  is  a  p-center  problem,  which  is  to 
locate  p  new  facilities  in  the  plane,  for  a  given  value  of  p,  with 
respect  to  a  finite  number  of  existing  facilities.   There  is  a  cost 
defined  for  each  existing  facility  which  is  the  minimum,  over  the  new 
facilities,  of  linear  functions  of  the  rectilinear  distance  between 
the  existing  facilities  and  each  new  facility.   The  problem  is  to 
locate  the  new  facilities  so  that  the  maximum  cost,  over  the  existing 
facilities,  is  minimized. 

An  example  of  a  p-center  problem  is  to  locate  p  points  in  the 
plane  so  that  the  maximum  rectilinear  distance  between  each  existing 
point  and  the  nearest  new  point  is  minimized. 

Candidates  for  the  optimal  objective  function  values  of  the 
p-center  problems  for  all  values  of  p  are  determined.   For  each  candi- 
date, a  set  covering  problem  is  formulated  whose  solution  determines 
a  feasible  solution  to  a  p-center  problem  for  some  value  of  p.   By 
solving  a  set  covering  problem  for  each  candidate,  optimal  solutions 
to  the  p-center  problems  for  all  values  of  p  are  determined. 
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4.2   Formulation 

Suppose  p  is  a  fixed  integer.   Let  X.,  1  <  j  <  p,  denote  the 

points  in  E  to  be  determined  for  the  location  of  p  new  facilities 

called  centers  and  let  the  set  X  =  fx ,X  "I  denote  a  p-center. 

-pi      p  

As  in  the  previous  chapters,  S  will  be  a  finite  set  in  E_  of  m  points 
P.  ,  1  <  i  <  m,  representing  the  locations  of  existing  facilities. 
For  each  point  P.,  let  w.  be  any  real  positive  constant  and  let  g.  be 
any  real  constant. 

For  each  i  and  j,    1  <  i  <  m,  1  <  j  <  p,  define  a  continuous 
function  in  the  p  variables  X.  by 


The  p-center  problem,  to  be  studied  in  this  chapter,  may  be 
written  as  follows: 

v  min  v  i^L  i<min  Cf--<xi x  >}•  ^•2-i> 

X  , .  .  .  ,  X   l<i<m  l<j<p  <•  j  l   1      p  i 
1      P 

To  show  that  an  optimal  solution  to  (4.2.1)  exists,  note  that 
the  function 

1      p    l<i<m  l<j<p   ji   1      p 


is  continuous.   Now  consider  the  region  in  E„  determined  by  the  rectangle 
containing  the  points  P.  with  sides  of  minimum  length  and  parallel  to 
the  x  and  y  axes,  respectively.   Any  point  X.  outside  this  region  may 
be  moved,  with  no  increase  in  the  objective  function  value,  in  directions 
parallel  to  the  x  or  y  axes  until  X.  is  contained  in  the  rectangle. 

<J 
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Therefore  it  is  sufficient  to  show  that  there  is  an  optimal  solution 

within  the  rectangle.   Since  the  rectangle  is  compact  in  E  ,  the  set 

in  E_   formed  by  p  cross-products  of  these  rectangles  is  also  compact. 
2p 

Therefore  f(X  , ... ,X  )  attains  its  minimum  value  in  this  compact  set 
1      P 

in  V 

Notice  that  for  p  =  1,  problem  (4.2.1)  is  a  special  case  of 
a  one  facility  minimax  location  problem  considered  in  Chapter  2. 

The  p-center  problem  may  be  written  equivalently  as  the  con- 
strained problem: 

min  z 

s-t.     ™JJ  w.j|x.  -P.J|+g.  <z     1  <  i  <  m.         (4.2.2) 

A  geometrical  motivation  for  the  p-center  problem  may  be  given 
using  the  following  definition.  A  diamond  centered  at  P.  and  of  value 
z   denoted  by  D(P.,z  )  and  is  defined  by 

D(P.  ,z*)  =  {X  e  E2  :  w.jj  X  -  PiJ|1  +  gj,  <  z*}  . 

Then,  in  reference  to  a  diamond  as  defined  in  [15],  D(P.,z  )  has 

* 
radius  (z  -g  )/w  . 
l    i 

The  p-center  problem  is  to  determine  the  smallest  value  z  for 

the  diamonds  D(P.,z)  and  to  locate  p  points  X.  within  these  diamonds, 

so  that  for  each  i.  1  <  i  <  m,  X   £  D(P  ,z)  for  some  j. 

J       i 

* 
Suppose  that  a  value  z   is  given.   Then  a  problem  related  to 

the  p-center  problem  is  to  determine  a  set  X  composed  of  a  minimum 

number  of  centers  that  satisfy  the  constraints  of  the  p-center  problem 
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(4.2.2)  with  respect  to  z*,  that  is,  determine  X  such  that  for  each  i, 


1  <  i  <  m  and  for  some  X.  €  X, 

J 


X.  e  D(P. ,z*). 


This  problem  is  now  formulated  as  a  set  covering  problem  and  its  rela- 
tion to  problem  (4.2.1)  is  established  in  Property  4.2.1. 

Consider  the  intersections  of  the  diamonds  D(P.,z*)  that  are 
determined  by  all  combinations  of  the  points  P.  taken  one  at  a  time, 
two  at  a  time,  and  so  forth,  up  to  m  at  a  time.   Let  M(z*)  be  the 
number  of  different  intersections  that  are  nonempty.   For  each  j, 
1  <  j  <  M(z*) ,  denote  the  j    nonempty  intersection  by  A.(z*).   Let 
I.(z*)  be  the  index  set  of  i  such  that  the  intersection  of  diamonds 
D(P.,z*)  equals  A. (z*).   Then 

A.(z*)  =   JJ    *   D(P.,z*)  . 
j       iei  .(z  )    i 

*  th 

For  any  X.  in  the  intersection  A.(z  ),  X   satisfies  the  i 
J  J       J 

constraint  for  i  e  I.(z*)  of  problem  (4.2.2).   The  problem  of  interest 

is  to  determine  the  minimum  number  of  intersections  A.(z*)  so  that 

J 

locating  a  center  X.  in  A.(z*)  has  the  property  that  each  diamond 
contains  a  center. 

Define  a  variable  x.  for  each  intersection  A  (z  )  by 

gfc 

f\  if  a  center  is  to  be  located  in  A.(z  ) 

I  0    otherwise  . 
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*  th 

Also,  an  m  component  vector  v . (z  )  is  defined  for  each  j,    where  the  i 

component  is  denoted  by  v..(z*)  and  is  defined  by 

C  1  if  i  e  I .  (z*) 

J 

v..(z*)  = 

0    otherwise. 


The  vector  v.(z*)  is  said  to  characterize  the  intersection  of 
J  

the  diamonds  D(P.,z*)  for  i  e  I.(z*),  since  it  specifies  which  diamonds 
intersect  to  determine  A.(z*). 

■J 

Thus  if  v   (z*)  =  1,  a  center  X.  located  in  A.(z*)  would  be 

a  point  in  D(P.,z*)  for  all  i  e  I.(z*).   For  each  point  P.  there  must 

be  at  least  one  center  X.  such  that  X.  e  D(P. ,  z*);  this  constraint  is 

expressed  by 

M(z*) 

I  v..  (z*)  x.  :>  1 

for  1  <  i  <  m.   Then  the  set  covering  problem  determined  by  z*  may  be 
written  as 


M(z*) 

min 

I        x. 
J-l        J 

M(z*) 
s.t.     E  v.  .  (z*)  x.  S:  1    1  <  i  <  m      (4.2.3) 
J-l  J1      J 


x  =  0,1,  1  <  j  <  M(z*) 


Example  4.2.1.     T>1   =  (-.5,4),  P2  =  (-1,3),  P3  =  (0,3),  ?4  =  (2,3.5); 

w  _  i  g  _  o  for  1  <  i  <  4;  z*  =  1.   The  diamonds  D(P. ,z  )  are 
i       i  x 
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illustrated    in  Figure   4.1.      There   are   eight    different   nonempty   inter- 
sections  so   that   M(z*)    =    8,    with    I      =    [l],    \^   =    {2},    I      =    {3}, 
I4={4},       I5={1,2};       I6={1,3},       I7={2,3},      Ig=    {1,2,3}. 
The   corresponding    intersections   A.(z    )      are  given  by 
Vl(l)    =    [1,0,0,0],      v2(l)    =    [0,1,0,0],      v3(l)    =    [0,0,1,0], 
v4(l)    =    [0,0,0,1],       v5(l)    =    [1,1,0,0],       Vg(l)    =    [1,0,1,0],    v?    =    [0,1,1,0,0], 

v    (1)    =    [1,1,1,0].       Then   the   set    covering   problem  may   be  written   as 
8 


mm   xi  +    x2   +    x3   +    x4   +    x5   +    xg   +    x?   +    xg 


10001101 
01001011 
00100111 
00010000 


h\ 


I  U    w 


jr  =  0,1     1  <  i  <  8 

l 


The  optimal  solution  is  x,  =  x„  =  1  and  all  other  x.  =  0.   Thus 
K  4    8  l 

any  point  in  A,(l)  and  any  point  in  A  (1)  determine  a  2-center  for  the 
4  o 

diamonds  D(P.  , 1) . 

Property  4.2.1.   Suppose  that  for  any  fixed  value  of  p,  the  optimal 

objective  function  value  of  problem  (4.2.1)  is  z  .   Let  z  =  z   and 

suppose  x  -    (x  , . . . ,x    *  )  is  an  optimal  solution  to  problem  (4.2.3) 
—       1       M  ( z  ) 

with 


M(z*) 

I  x°  =  p. 
j  =  l   J 
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■>•   x 


Figure  4.1      Diamonds   D(P.,1) 
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Define  the  p-center  X  =  fx X  }   as  follows:   if  x.  =  1,  X. 

-P   l  1      PJ  J       J 

*  o 

is  any  point  in  A.(z  ).   Then  X   is  an  optimal  solution  to  (4.2.1). 

J 

Proof :   Since  x   is  a  solution  to  the  set  covering  problem,  for  each 

i  there  is  at  least  one  x.  -    1  such  that  v..(z*)  x.  =  1.   That  is, 

J  Ji      J 

for  each  i,  there  is  at  least  one  A.(z*)  such  that  for  X.  in  A.(z  ) 

w  ||  X°  -  P.  ||  +  g.  <  z*  =  z°. 
i    J     i      i 

Therefore,  for  each  i,  1  <  i  <  m 

min     II  o   _  ||       _.  *    o 
j^pWjIXj  -Pjl  +g.  <z  =z 

and  the  set  X   and  the  value  z   is  a  feasible  solution  to  (4.2.2)  with 

o  o 

p  centers.   Since  z   is  the  minimum  value  of  z ,  X   is  an  optimal  solu- 
tion to  (4.2.2)  and  hence  to  problem  (4.2.1). 


Corollary  4.2.1.   Suppose  the  hypotheses  of  Property  4.2.1  hold  but 

o  *    o 

that  the  optimal  solution  x   to  (4.2.3)  for  z   =  z  has 

M(z*) 

l.     x.  =  r  <  p. 
j=l  J 

Then  the  set  X°  =  [X°   . .  ,X°]  U  fX^ .X^    }  .where  X°  is  defined 

—     -  1      r-*    -  Jl      Jp-rJ  J 

in  the  statement  of  Property  4.2.1  and  where  Xj   is  any  point  in  any 

diamond  D(P.,z  ),  is  an  optimal  solution  to  problem  (4.2.1)  for  a  fixed 

value  of  p. 

Corollary  4.2.1  considers  the  case  where  z   is  an  optimal 

solution  to  a  p-center  problem,  for  a  fixed  value  of  p,  but  only 


75 


r  centers,  r  <  p,  are  needed  for  an  optimal  solution.   In  this  case  the 
set  covering  problem  determines  only  r  centers.   To  obtain  a  p-center 
solution,  any  p-r  additional  points  X.,  as  specified  above,  may  be 
used  in  the  p-center  solution. 

Hence,  for  a  fixed  value  of  p  and  an  optimal  objective  function 
value  z   to  problem  (4.2.1),  Property  4.2.1  implies  that  an  optimal  solu- 
tion to  problem  (4.2.1)  may  be  determined  using  the  optimal  solution  to 
problem  (4.2.3).   The  procedure  for  solving  (4.2.1)  will  be  to  deter- 
mine candidates  for  the  optimal  objective  function  values,  and  for  each 
candidate  to  solve  the  corresponding  set  covering  problem. 

4. 3  Applications 

Applications  of  the  p-center  problem  may  occur  in  contexts 
similar  to  the  minimax  location  problems  studied  in  the  previous  chap- 
ters, that  is,  in  locating  emergency  service  facilities  in  settings 
where  travel  distance  is  approximated  by  the  rectilinear  distance 
between  facilities.   However,  in  the  location  problems  of  previous 
chapters,  a  new  facility  is  located  so  that  it  can  respond  to  any 
existing  facility;  in  a  p-center  problem  a  center  is  located  so  that 
it  can  respond  only  to  those  facilities  nearest  to  it. 

As  an  example  [18],  suppose  there  are  m  existing  facilities 
and  there  are  p  policemen  available  for  emergency  call  to  the  facil- 
ities, (p  <  m).   It  is  desired  to  locate  the  policemen  so  that  for  any 
emergency  occurring  at  an  existing  facility,  the  maximum  distance  that 
the  nearest  policeman  has  to  travel  is  minimized.   Similar  applications 
may  occur  in  locating  fire  stations,  ambulance  centers,  or  other 
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emergency  services  so  that  the  maximum  distance  from  any  existing  facil- 
ity to  the  nearest  service  facility  is  minimized. 

Applications  may  also  be  found  in  locating  centers  that  empha- 
size accessibility.   Suppose  it  is  desired  to  locate  p  mail  boxes  in 
a  city  with  respect  to  m  population  groups  idealized  as  points  P. , 
1  <  i  <  m.   The  mail  box  locations  are  to  be  determined  so  that  the 
maximum  distance  anyone  at  some  point  P.  has  to  travel  to  the  nearest 
mail  box  is  minimized. 

In  p-center  problems,  the  constants  w.  and  g.  may  have  similar 
interpretations  to  those  discussed  in  Chapter  2.   That  is,  the  w.  may 
represent  cost  or  some  other  constant  per  unit  of  distance,  and  the  g 
may  represent  fixed  costs  or  constants. 

4.4   Related  Work 

Hakimi  [18]  appears  to  have  originally  formulated  the  p-center 
minimax  problem.   His  formulation  was  a  p-center  location  problem  on 
a  graph,  and  no  solution  procedure  was  offered. 

Christofides  and  Viola  [3]  consider  the  p-center  problem  for 
a  graph  and  develop  a  solution  procedure.   The  approach  developed  in 
this  chapter  for  the  p-center  problem  in  the  plane  is  similar  to  the 
approach  taken  on  the  graph  problem  by  Christofides  and  Viola.   On  the 
edges  of  the  graph,  they  determine  regions  that  are  defined  by  the 
intersections  of  "rays"  of  fixed  length  along  the  edges  emanating  from 
each  node.   A  set  covering  problem  is  then  solved  to  determine  the 
minimum  number  of  regions  needed  to  cover  the  nodes  of  the  graph. 
By  incrementing  the  length  of  the  rays,  a  sequence  of  set  covering 
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problems  is  solved.   This  approach  may  not  determine  the  optimal  objec- 
tive function  value  for  some  p-center,  but  can  iterate  to  an  objective 
function  value  that  is  arbitrarily  close  to  the  optimal  value.   A  set 
covering  problem  is  solved  at  each  iteration. 

The  rays  that  are  used  in  [3]  on  the  graph  p-center  problem 
are  analogous  to  the  diamonds  that  are  used  in  Section  4.2.   One  result 
of  this  chapter  is  that  for  the  p-center  problem  using  rectilinear  dis- 
tances in  the  plane,  the  optimal  objective  function  value  is  determined 
for  each  value  of  p. 

Toregas  et  al.  [27]  consider  a  problem  of  locating  centers  in 
the  plane.   They  assume  that  centers  can  only  be  located -at  existing 
facility  locations.   Sets  of  existing  facilities  are  constructed  that 
contain  potential  center  locations  within  a  specified  time  or  distance 
of  each  existing  facility.   The  problem  of  locating  centers  within  these 
sets  to  respond  to  existing  facilities  is  formulated  as  a  set  covering 
problem. 

A  p-center  problem  using  general  distances  is  considered  by 
Hearn  [19].   He  develops  an  enumeration  algorithm  for  solving  the 
p-center  problem  for  a  specified  p  and  gives  extensive  computational 
experience. 

Total  cost  analogues  to  the  p-center  problem  have  been  referred 
to  as  location-allocation  problems.   Cooper  [  4 ]  develops  heuristic 
procedures  for  these  problems,  assuming  Euclidean  distances,  and  gives 
additional  references. 
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4.5   Problem  Development 

In  this  section,  properties  of  problems  (4.2.1)  and  (4.2.3) 
are  developed  that  form  a  basis  for  a  solution  procedure  for  problem 
(4.2.1). 

For  each  pair  of  points  P.  and  P.  for  1  <  i  <  j  <  m,  compute 
the  constants  B. .  by 


fe^r  [vj"  pi  "  pjHi +  giwj  +  gjwJ  -  gi-gj} 


B   =  max 
ij 

i   J 

The  number  of  B. .  computed  is  f  J  +  m,  which  is  denoted  by  m  . 

Next,  order  the  B  .  and  set  them  equal  to  t,  so  that  t,  <  t,  ,  for 
'  ij  k  k    k+1 

1  <  k  <  nu  -  1. 

o 
Property  4.5.1.   Suppose  p  is  a  given  integer  value,  p  <  m.   If  z   is 

the  optimal  objective  function  value  for  the  p-center  problem,  (4.2.1), 

then  z  =  t,  for  some  1  <  k  <  ni  . 
k  2 

Proof:   Let  X  =  fx, ,...,X  }  be  an  optimal  solution  to  problem  (4.2.1). 
-p     1      pJ 

Define  the  index  sets  N  for  1  <  v  <  p  by 

Nv=  f1  :  xv  e  D(pi>z°>}  • 

Consider  the  p  one  facility  minimax  location  problems,  written  as 
constrained  problems, 


mm   z 

v 

(4.5.1) 
s.t.         w.||Xv   -   P±\\1    +   g.    <   zv  i    C   Nv 
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for  1  <  V  <  p.   By  Property  2.8.2,  the  optimal  objective  function  value 

th  v 

f  the  V   problem  equals  some  B. .,  say  B. .,  which  equals  some  tfc  for 


o 

1  <  k  <  m  . 

By  definition  of  N'v,  X   and  z  are  a  feasible  solution  to  the 
V   problem  (4.5.1)  so  that     B. .  <  z   for  1<  V  <  p.   Let 


max   ,  v  , 
]<  v<p  l  Bijj 


Then  z'  <  z  . 


Let  X'  be  an  optimal  solution  to  the  V   problem  (4.5.1)  and 

consider  the  set  X'  =  [ X' X' } .   For  each  i,  1  <  i  <  m,  i  e  N 

—plP  v 

for   some   V,    that    is,    for   each    i,    there   is   some   X'    that    satisfies    the 
expression 

w.||   X^   -   P.||1    +    g.    <   z'     . 
Therefore,    for   each    i,    1  <   i  <  m, 

"**         W     J|X'      -    P     II  +     g.     <    2' 

1<  v<p      i"     \>  i"l         &i 


and  the  set  X   and  z  is  a  feasible  solution  to  problem  (4.2.2). 
-P 

Hence  z  <  z' ,  so  that  z   =  z'  =  t   for  some  1  <  k  <  m  . 

K  £ 


By  Property  4.5.1,  some  value  t  ,  for  1  <  k  <  m  ,  will  be  the 

optimal  objective  function  value  to  (4.2.1)  for  any  value  of  p. 

By  Property  4.2.1,  for  each  value  of  t   that  is  an  optimal  objective 

function  value  to  problem  (4.2.1)  for  a  fixed  value  of  p,  the  set 

covering  problem  (4.2.3)  generated  by  t   determines  an  optimal  p-center 

k 
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for  problem  (4.2.1).   By  solving  problem  (4.2.3)  for  each  value  of  t 

K 

in  increasing  order  of  t  ,  all  p-center  solutions,  for  p  =  m,m-l,...,2, 
to  (4.2.1)  are  obtained.   For  each  value  of  p,  the  smallest  value  of 
t   that  gives  a  p-center  solution  to  (4.2.2)  is  the  optimal  objective 
function  value  of  (4.2.1). 


4. 6   Determining  the  Set  Covering  Problems 

Attention  is  now  directed  toward  deriving  the  column  vectors 

for  the  set  covering  problem  determined  by  the  parameter  t  .   A  proce- 

dure  is  developed  that  determines  a  set  of  vectors,  denoted  by  V(t  )  , 

that  contains  vectors  v.(t  )  which  characterize  all  the  nonempty  inter- 

J   k 

sections  of  diamonds  D(P. ,t  ),  1  <  i  <  m,  that  are  needed  for  the  set 
covering  problem  generated  by  t  .   The  vectors  in  V(t  )  are  used  as 
the  column  vectors  in  the  set  covering  problem  determined  by  t  . 

K 

The  problem  of   determining   the   set   V(t    )    of   vectors   v . (t    )  ,    l<j<M(t    ), 

K  J    K  K 

is  decomposed  into  two  subproblems;  each  subproblem  determines  a  set  of 

vectors,  denoted  by  V  (t,  )  and  V  (t,  ),  respectively,  that  characterize 

x  k       y   k 

all  the  intersections  of  intervals  along  the  x  and  y  axes. 

For  a  given  value  of  t  ,  consider  the  diamonds  D(P. , t  )  for 

K  1    K 

1  <  i  <  m.   Applying  the  transformation  T  rotates  each  diamond  onto 

a  square  with  sides  parallel  to  the  rotated  axes  x  and  y,  respectively. 

Define  a  square,  in  the  rotated  space,  centered  at  P|  =  (a!  ,b'.  )  and  of 

value  t,'  as  the  set 
k 

S(P!  ,t^)  =  {X'  e  E2  :  w.|!  X'  -  P'l!^  +  g'  <  t^} 
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where  g1  =  g  /-J2   and  t'  =  t  A/2  ,  corresponding  to  the  results  in 
i    i         k    k 

Chapter  2.   Note  that  a  square  S(P*ft')  has  side  length  2<t^-g^)/wi< 
Then  each  nonempty  intersection  of  diamonds  D(pi»tk^»  for 

i  e  I  (t  ) ,  corresponds  uniquely  to  the  nonempty  intersection  of 
J   k 

squares  S(P',t')  for  i  !I,(t,),   Consequently,  the  procedure  actually 
M         i'  k  J   k 

determines  vectors  that  characterize  the  nonempty  intersection  of 
squares  in  the  rotated  space;  these  vectors  are  denoted  by  ^(t^)  for 
1  <  j  <  M(t^). 

The  problem  of  determining  the  vectors  v.(t')  may  be  decom- 

J   k 

posed  into  one-dimensional  problems  of  determining  vectors  that 
characterize  the  intersection  of  intervals  on  the  x1  and  y'  axes, 
respectively,  where  the  intervals  are  obtained  by  projecting  the  sides 
of  the  squares  onto  the  respective  axes.   Explicitly,  an  interval  cen- 

X  X 

tered  at  a!  of  value  t   is  denoted  by  I  (a'.  ,t  )  and  is  defined  by 
I(a^,tX)  =  fx'  e  E1    :  w  |x'  -  a^|  +  g^  <  tX]  . 

X  X 

Notice  that  an  interval  I(a!,t  )  has  length  2(t  -g!)/w..   Next,  the 
problem  of  determining  intersections  of  intervals  is  studied  and  a 
computational  procedure  is  developed. 

It  will  be  shown  that  the  vectors  in  V(t  )  are  determined  from 
the  vectors  that  characterize  intersections  of  intervals  on  both  axes. 

Example  4.6.1.   Consider  Example  4.2.1  and  apply  the  transformation  T 
to  obtain  the  squares  as  illustrated  in  Figure  4.2.   The  unique  corre- 
spondence of  intersections  among  diamonds  and  among  squares  may  be 
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Figure   4.2      Squares    S(P'.  ,1/*J2) 
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observed  from  Figures  4.1  and  4.2.   Figure  4.2  also  illustrates  the 
intervals  and  their  intersections,  corresponding  to  sides  of  squares, 
along  the  x'  and  y'  axes. 

Now  consider  the  problem  of  characterizing  the  intersections 
of  intervals  I(a.,t  )  by  zero-one  vectors.   In  order  to  define  these 
vectors  precisely,  it  is  necessary  to  establish  the  following  three 
lemmas. 

For  each  pair  of  points  a'  and  a!,  1  <  i  <  j  <  m,  compute  the 

m  constants  b   defined  by 
2  lj 

b.  .  =  max  i w.w.  a'.  -  a'.   +  w.g'.  +  w.g!  ,  g!  ,g'.  f  . 

ij       Lw.+w.  L  i  J   i    J1     i  J    J  ij   i   jJ 


X  X 

Lemma  4.6.1.   Suppose  i  4   j,    b.  .  >  g!  and  b  .  >   g'.  with  the  points 
ij     i      iJ     J 

a!  and  a'  named  so  that  a'  <  a' .   Define  the  point  x'   by 
i      j  i    J  ij 


x!  .  =  w.a'  +  w  a'  +g'  -g' 

ij   w  +w  L  i  i    J  J    J    ij 


i  y 


Then,  x!  .  is  the  unique  point  such  that 


w.|a!  -  x!  .1  +  g'.  =  w.l  a'.  -  x'  I  +  g'.  =  bX . , 


X  X 

Proof:   The  inequality  b  .  >  g'  and  the  definition  of  b    imply 
ij     i  ij 

w 


bX.  -  g!  =  —  |w .la!  -  a'.|  +  g'.  -  g!  |>  0  .     (4.6.1) 

ij    i   wx+wy  L  J   i    j'    J    ij 

A  similar  expression  with  i  and  j  interchanged  and  denoted  by  (4.6.2) 
results  from  the  inequality  b.  .  >  g'. .   The  expressions  (4.6.1)  and 

■J         J 

(4.6.2)  may  be  used  to  show  that  a'  <  x' .  <  a1  by  showing  that 


and 
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x!   -  a'  =  (bx.  -  g!)  w.  >  0  (4.6.3) 

ij     i     ij     i    i 


a'  -  x'  .  =  (bX.  -  g'.)/w.  >  0  .  (4.6.4) 

i     ij      iJ     J    J 


But  expressions  (4.6.3)  and  (4.6.4)  imply  that 


w  x'  .  -  a1  +  g'  =  w.x'.  -  a'.  +  g'.  =  b... 

i   ij    i  i    J   iJ  J  J    iJ 

Then,  these  equations  may  be  used  to  show  that  x'.  .  is  the  unique 
such  point. 


X  X 

Lemma  4.6.2.   Suppose  i  *  j,  and  b  J  -    g'  so  that  b.  .  ^   g'..   Define 
ij     i  ij     J 

the  point  x!  .  =  a'..   Then, 

w  Jx'  .  -  a*  I  +  g'  =  bX  and    w.lx!.  -  a'.|+  g'.  <  bX . . 
i   ij    i1    Bi           ij  j1  iJ    jl    J    iJ 


Proof:   The  assumptions  imply  that 


g '  ^  w  w  .  a '  -  a '.   +  w .  g '.  +  w  .  g ! 

6i    w.+w.  L  i  j1  i    J      i  J     J  iJ 


i   J 

This  inequality  implies  w.  a!  -  a'.   +  g'.  <  g!  =  b.  ..   Clearly, 

J   i    J     J     i     iJ 

i       i  x 

w  I  a'  -  a'  +  g'  =  s1  =  b   .   Hence  the  lemma  is  proved, 
i   i    i1    i    i    ij 


x  x 

Lemma  4.6.3.   Suppose  i  *. r  j  ,  and  b.  .  =  g'.  ,  so  that  b.  .  ^  g!  . 
iJ     J  iJ     i 

Define  the  point  x!  .  =  a'..   Then, 

ij    i 

w  |x*   -  a'  |  +  gf  <  bX.    and    w.lx'.  .  -  a'.  |  +  g'.  =  bX . , 
i   ij     i'     i     ij  J   iJ     J     3  iJ 
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By  the  above  three  lemmas,  if  i  4   j,  then  the  intersection 

X  X 

I  (a'  ,b   )  H  I  (a'  ,  b'  )  is  nonemDty   since  it  contains  x'.  ..   Note  also 
i   ij       j   ij  "  !J 

that  if  i  -  i,  then  b'.  .  =  g'.  ,  and  x'.  .  =  a!  is  in  the  interval 

11     1       hi 

I (a' ,b   ).   Hence,  for  any  i  and  j,  1  <  i  <  j  <  m, 
i   ii 

I  (a!  ,bX.)  n  I  (a'.  ,bX.)  4   0. 

x  x 

For  each  b   ,  define  a  vector,  denoted  by  v(b. .),  that  charac- 
ij  ij 

X  X 

terizes  the  intersection  I(a!,b'.  .)  D  I(a'.,b.  .),  where  the  characteriza- 

i   ij       3      iJ 

tion  is  provided  by  having  component  k  equal  one  if  and  only  if 

I(a'  ,bX  )  fl  I(a',bX  )  n  I(a',bX.)  *  0.   Clearly  components  i  and  j  of 
k   ij       i   ij       3   ij 

v(b   )  equal  one.   The  next  property  provides  a  procedure  for  deter- 
ij 

x 

mining  all  the  components  of  v(b. .). 

Property  4.6.1.   The  k    component  of  v(b. .)  equals  one  if  and  only  if 

XX  XX 

b    .    s  b.,       and     b.  .    >  b.,  . 
ij  ik  ij  jk 

XX  XX 

Proof:   Suppose  that  b   s  b'    and   b   ^  b  ,  .   It  must  be  shown  that 
ij     ik      .   ij     Jk 

x 

for  some  point  x'  in  the  intersection  characterized  by  v(b. .), 

ij 

x'  e  I  (a,  ,bX.). 
k   ij 

x  x 

Case  A.   i  =  j.   Then,  b'.  .  =  g'.  and  the  assumption  gives  g'  ^  ^v- 

ill  l  ik 

Thus   the  following   inequality  holds: 

1        fw.wja!    -   a'  |    +   w.g*    +   w  g!    I  <  g!  ,  (4.6.5) 

*'  •  +w,    liki  k  ik  k    l  I  l 

i      k 

which  implies  w,  a'  -  a'  I  +  g,'  <  g'   and   a!  e  I  (a,1  ,g! ),  so  that  the 
kik     ki        i      ki 

th 
k   component  is  one. 
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Case  B.   i  4   j.   Then  by  Lemmas  (4.6.1),  (4.6.2)  and  (4.6.3), 

XX  XX 

I  (a'  ,b'  )  PI  I  (a'  ,b   )  4   0.   Four  cases  in  terms  of  the  b' ,  and  b" 
i   ij       j   ij  ik      jk 

need  to  be  considered. 


x  x 

la:   b'   >  g'   and   b   >  g' .   Then  i  4   k,  and  Lemma  4.6.1 
ik     i        ik    k 

implies   I(a'.  ,bX  )  "    I(a,' ,b._  )  x   0;    since  b.,  <  b*.  .  , 
l   ik       k   ik  ik    ij 

I(a!  ,bX.)  n  I  (a'  ,bX.)  x   0. 
i  ij       k   ij 

lb:      b'       >  g'       and      b'      >  g,1  .      Then,    similarly, 
Jk        "j  Jk  k 

I(a'.,bX.)    0   I(a'  ,bX.)    x    0    . 
J       ij  k      ij 

2a:   bX   a  g! ,  bX  =  g' .   If  i  x    k,  Lemma  (4.6.3)  implies 

IK       1     IK       K 

that   I(a*    b*    )    PI  I(a*  ,b*    )   ji  0      so   that      I  (a!  ,bX .)    fl  I  (a*    bX .)    d  0 , 

i   ik       k   ik  i   ij       k   ij 

XX  X 

since  bv  assumption  b"   <  b   .   The  same  conclusion  holds  if  b.,  =  g! 

ik     ij  ik    l 

and   bX   ^  g'   by  Lemma  (4.6.2).   If  i  =  k,  b.  .  =  g'  and 
ik    k  n     l 

I(a! ,bX.)  0  I(a! ,bX  )  i   0. 

i   n       i   n 

2b:   bX,  >  g'.  ,  bX  =  g*    or   bX  =  g'.  ,  bX  2  g'    Then 
jk   &j    jk   bk'       jk    J    Jk    k 

x  x 

analogous  to  Case  2a,  I(a'.,b.  .)  fl  I(a'  b.  .)  =  0. 

j   ij        k'  ij 

Therefore  in  any  pairing  of  cases  a  and  b  the  resulting  non- 
empty intersections  of  intervals,  combined  with 

I  (a!  ,bX.)  fl  Ka'.  ,bX.)  4   G, 
i   ij       J   iJ 


imply  that 


I  (a*.  ,bX.)  D  I(a'.,bX.)  H  I(a'  bX.)  x   0. 
i   ij       J   ij       k  ij 
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XX  X 

Hence   there   is    an  x'    in   I(a!,b..)    fi   I(a'.  ,b*.  .)    that    is    also    in   I(a'    b    .) 

i      ij  J      t-J  k      ij 

th  x 

Conversely,    if   the  k        component   of   v(b.  .)    is   one,    then 

IIX  X  P)  X 

x'    -   a'      +   g'    <  b    .    for   some  x'    in    I(a'.  ,b    .)    M   I(a'.  ,b    .).       Suppose 
k  k'  k  ij  1      ij  j      ij 

b„    =   w.w    |  a'.    -   a'  |    +  w.g'    +   w  g'. 

lk        W-+WV  Liki  k  lk  kij 

Then  the  triangle  inequality  gives 

b       <  w . w      a '    -  x'      +   v.-  g !    +   w  w      a '    -   x1       +  w .  g' 

lk        w.+w,    Lik'i  '  ki  ik'k  l   k J 

l      k 


i r  *,*      »*  ~i 

w  b .  .    +   w .  b .  . 

w.+w.    L  k    ij  i    ijJ 


< 

V"k 


=   bX.. 


X  XX 

If   b.,  =  g! ,  then  by  definition  of  b  .,  g'  <  b  .. 
ik    i  ij   °i    ij 

Finally,  if  b*fe  =  g^ , 

g^wjx'  -aj  +g'<bX.. 

X       X 

A  similar  argument  proves  the  second  inequality,  b   <  b.  .,  and  the 

Jk    ij 

proof  is  complete. 

For  a  given  value  of  b'  ,  the  next  property  gives  a  method  of 
determining  all  the  vectors  that  characterize  the  collection  of  non- 
empty intersections  of  intervals  I (a' , b*  )  for  1  <  i  <  m. 

l   pq 
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Property  4.6.2.   Let  V  (b'  )  =  fv(t')  :  t,  <  bX  }  for  a  given  value  of 
x  pq     "•    k     k    pqJ 

X  X       1    X  X 

b   .   If  v(t,  )  e  V  (b'  ),  than  v(t"  )  characterizes  a  nonempty  inter- 
pq         k     x  pq  k  J 

x 

section  of  the  intervals  I(a'.  ,b   ).   Conversely,  each  nonempty  inter- 

x  x 

section  of  the  intervals  I(a! ,b"  )  is  characterized  by  a  vector  v(t  ), 

x  _  ,  x 

for  some  t.  S  b 

k    pq 

Proof:   Suppose  v(t,  )  e  V  (b   ).   Since   t,  <  b'X  ,  I  (a*  ,tX)  C  I  (a'  ,bX  ) 
k      x   pq  k     pq      l   k        l   pq 

for  all  i,  so  that  any  intersection  of  the  intervals  I(a!,t  )  for  some 

x 

subset  of  the  i  will  be  contained  in  the  intersection  of  I (a! ,b   )  for 

i  pq 

the  same  i.   Therefore,  v(t  )  characterizes  a  nonempty  intersection  of 

X  XX 

the  intervals  I (a' ,b   ) ,  for  t'  <  b'  . 
i   pq        k    pq 

Now  suppose  G  is  a  subset  of  the  indices  i,  1  <  i  <  m,  such 
that 


,2-  I (a! ,bX  )  *  0  . 
ieG    i   pq 

Consider  the  points  a'.  ,  i  e  G,  as  fixed  points  in  a  one  facility  minimax 
location  problem.   Then  the  optimal  objective  function  value  is 


»x,  =  ,7r  {bx.}  , 

St  1.J5G     I    ljJ 


and  b'   is  the  minimal  value  such  that  the  intervals  I(a! ,b  ,)  have 
st  1   st 

XX  XX 

nonempty  intersection.   Thus  b   <  b'    and   b  .  =  t,   for  some  k. 

st    pq        st    k 

By  Property  4.6.1,  v(b   )  characterizes  the  nonempty  intersection  of 

I (a' ,b   )  for  i  in  G. 
l   pq 

Some  of  the  vectors  characterizing  intersections  of  intervals 


may  be  redundont.   For  example,  if  the  intersection  corresponding  to 

-v  XX 

v(tj|   )    contains    the   intersection   corresponding   to   v(ti<0)  ,    then  vC^.,] 
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is  not  needed.   In  this  case,  the  vector  v(t^1)  has  a  one-component 
for  each  one-component  of  v(tk  ). 

Redundant  vectors  may  be  identified  using  the  following  defin- 

X  X 

itions.   A  vector  v(tj.  )  is  said  to  dominate  a  vector  v(tjj  )  if 

X  X 

v(t^  )  ^  v^ko^  *n  "^e  vector  sense  [17].   Then  a  dominated  vector 

v(t  )  is  redundant.   An  equivalent  expression  of  dominance,  and  one 
K 

that  is  useful  computationally  is  given  by  the  following  definition. 

The  Boolean  product  [3]  of  two  zero-one  vectors  v   and  v  ,  denoted  by 

th 
v   °  v  ,  is  the  vector  whose  k   component  is  the  product  (logical  and) 

th  x 

of  the  k   components  of  v  and  v  ,  respectively.   Then  v(tk  ) 

x         £  ± 

X  XXX 

dominates  v(t^  )  if  and  only  if  v(tk  )  o  v(tk  )  =  v(t^  ). 

Let  V  (b  .)  be  a  subset  of  the  vectors  V  (b  . )  =  fv(t  )  : 
x  ij  x  ij        k 

XX  X 

t  <  b  },  with  the  property  that  no  vector  in  V  (b.  .)  is  dominated  by 
k    ij  x  ij 

x  x  1  x 

any  other  vector  in  V  (b. .).   Then  V  (b. .)  is  obtained  from  V  (b  .)  by 

x  ij  x  ij  x  ij 

eliminating  all  the  dominated  vectors.   The  vectors  in  V  (b   )  still 

x  ij 

characterize  the  intersections  of  intervals  I(a.,b  .). 

An  extension  of  the  set  V  (b  .)  is  now  considered  that  has 

x   ij 

application  later  in  the  sequel. 

Remark  4.6.1.   Suppose  t  is  any  real  number  and  consider  the  set  V  (t  ) 

x     p 

X  XX 

where  t      is   the  largest  value  of   t .  ,    1  <  i  <  ni      such   that   t     <  t,    then 
P  l  2  P 

V  (tX)  =  V  (t). 
x  p     x 

X  X 

Proof:   Since  t  <  t ,  I (a! , t  )  c  I (a! ,t)  for  1  <  i  <  m,  and  any  inter- 

p         i  p       l 

section  of  I(a!,t  )  is  contained  in  the  respective  intersection  of 
i  P 

I (a! ,t);   thus  V  (tX)  C  V  (t). 
i  x  p     x 
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Suppose  v(t)  is  a  vector  in  V  (t) ;  let  G  be  the  set  of  indices 
i  such  that  component  i  in  v(t)  equals  one.   Suppose 


bx  =   max   fx  | 

pq    i,j£G  L  ljJ 


Then,  b   <  t  by  the  logic  in  the  proof  of  Property  4.6.2.   Therefore, 

v(t)  is  in  V  (tX). 
x   p 

The  procedure  developed  above  may  also  be  used  to  determine  the 
set  V  (t)  of  vectors  characterizing  the  intersections  of  intervals 

y 

I(b! ,t)  on  the  y1  axis  for  any  real  value  of  t. 

The  above  development  provides  a  method  for  characterizing  the 
intersections  of  intervals  of  a  specified  value  on  either  axis.   Atten- 
tion is  now  directed  toward  characterizing  the  intersections  of  squares 

x  y 

using  the  sets  of  vectors  V  (t  )  and  V  (t,  )  that  characterize  interval 

x  k       y   k 

intersections  in  each  axis. 


Property  4.6.3.   v.(t')  is  a  vector  characterization  of  an  intersection 

j   k 

of  squares  in  the  rotated  space  if  and  only  if  there  exist   a  vector 
v(t£  )  in  V  (t')  and  a  vector  v(tk  )  in  V  (t')  such  that 


V_.(t,'_)  ::  V(tk  )   o  V(tk,) 


x  y 

Proof:   Suppose  there  exist  vectors  v(tk  )  in  V  (t')   and   v(tk,p)  in 

V  (t,' )  whose  Boolean  product  is  v  .  (t,1  ) .   Let  G  be  the  set  of  indices 
y  k  J  k 

x  y 

i  such  that  component  i  equals  one  in  both  v(tk  )  and  v(t"k9)  ;  consider 

the  sauares  with  center  (a' ,b')  and  with  value  t'  for  i  in  G.   Let  x' 

a   l  k 

be  a  point  in  the  intersection  characterized  by  v(tk  )  and  let  y'  be  a 

y 

point  in  the  intersection  characterized  by  v(tk_).   Then  for  i  in  G, 
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w.Ja-  _  x.|  +  g-  <  t^<  t£ 

and  w.lb!  -y'|  +  g^<  tj2<  t£. 

Thus  w.N  (a!,b!)  -  (x',y')|L  +  g!  <  t' 

1     11  "■     l     k 

and  the  point  (x' ,y' )  is  in  the  intersection  of  the  squares  S( (a! ,b' ) , t' ) 

for  i  in  G.   That  is,  v.(t')  characterizes  an  intersection  of  squares 

J   k 

in  the  rotated  space. 

Conversely,  suppose  v  .  (t ' )  is  a  characterization  of  an  inter- 

J   k 

section  of  squares.   Let  the  set  of  components  i  of  v.(t')  that  equal 

J   k 

x  v 

one  be  denoted  by  G.   Consider  the  intervals  I(a! ,t  )  and  I(b' ,t  ),  for 

l   p         i   q 

x      y  xv 

i  in  G  and  where  t   and  t   are  the  largest  values  of  t   and  t  ,  respec- 

p      q  k      k 

x  y 

tively,  such  that   t  <  t'  and  tJ   <   t'         Suppose  that  the  intervals 

p    k      q    k 

x 

I(a! ,t  )  have  empty  intersection  for  i  in  G,  then  an  argument  similar 

to  the  last  half  of  the  proof  of  Remark  4.6.1  shows  that  the  intervals 

I(a!,t' )  have  empty  intersection  for  i  in  G.   But  this  implies  that  the 
1   K 

squares  S(  (a!  ,b!  )  ,  t'  )  have  empty  intersection.   Therefore     l(a't")^'0. 
ilk  i€G    i   p 

A  similar  argument  shows   that   H   l(b',t  )  ^  0. 

icG    i   q 

x  y 

Let  v(t  )  and  v(t  )  be  the  respective  vector  characterizations 
P        q 

of  these  nonempty  intersections  of  intervals.   Then,  since  the  two 

intersections  immediately  above  are  nonempty,  the  i    component  of 

each  equals  one  for  i  in  G,  which  implies  v(t  )  o  v(t  )  -  v  ft'). 

P       q     J  k 

Thus  for  a  fixed  value  of  t ' ,  the  Boolean  product  of  a  pair  of 

vectors,  one  from  V  (t' )  and  one  from  V  (t' )  ,  generates  a  vector  v  (f ) 

x  k  y   k  j   k 

that  characterizes  an  intersection  of  squares  of  value  t' ,  and  hence  of 
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diamonds  of  value  t  .   The  set  of  all  vectors  generated  by  all  such 

JEC 

pairs  characterize  all  intersections  of  diamonds  with  radius  t  . 

k 

Suppose  that  all  the  dominated  vectors  in  this  set  are  eliminated. 
Then  the  resulting  set  of  vectors,  denoted  by  V(t  ),  is  the  smallest 

K 

number  of  vectors  that  completely  characterize  all  nonempty  inter- 
sections of  diamonds  of  value  t,  .   These  vectors  are  the  column  vectors 

k 

for  the  set  covering  problem  generated  by  t  . 

4.7   A  Solution  Procedure 

A  procedure  is  now  given  for  obtaining  a  solution  to  the  p-center 
problem  for  any  p,  p  =  m,m-l,... ,2.   For  p  =  1,  the  problem  is  a  one 
facility  location  problem  which  is  studied  in  Chapter  2. 

Note  also  that  an  m-center  solution  is  immediately  given  by 
locating  a  center  at  each  of  the  fixed  facility  locations  P  for 

1  <  i  <  m  and  has  an  objective  function  value  z  where  z  =  _   fg.f. 

The  following  algorithm  determines  a  solution  to  the  p-center 
problem  for  a  given  value  of  p,  say  p  .   It  also  determines  all  p-centers 
for  p  <  p  <  m. 


Algorithm 


x  x 

Compute  and  order  the  b  . . ;  set  them  equal  to  t  with 

ij  k 

XX  V 

t     <  t'       .       Compute   and  order   the   b. .:    set    them  equal    to 
k  k+1  ij 

t?  with   t^  <  tj   ,.      Set  b. .   =   max   (bX.,by.).      Order   the 
k  k  k+1  ij  ij      ij 

b    .    and   set   them  equal    to   t,    with   t,    <  t, 
ij  k  k  k+1 

Set   k  =    1,    p     =   m. 
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2.  Compute  the  sets  V  (t1 ) ,  V  (f )  and  the  set  V(t  )  from 

x  k    y  k  k 

the  Boolean  product  of  vectors  in  V  (t' )  and  V  (f), 

x  k       y   k 

respectively.   Solve  the  set  covering  problem  using 

vectors  in  V(t')  as  column  vectors.   Suppose  the 

optimal  objective  function  value  is  p  .   If  p  <  p  , 

o       ok 

go  to  3.   If  p  =  p,  ,  go  to  4. 
o    k 

3.  The  value  t   is  the  optimal  objective  function  value  for 

K. 

the  p-center  problem  with  p  <  p  <  p,  .   If  p  <  p„  stop. 

o        k       o    1 

Otherwise  set  p,  =  p  ,  k  =  k+1  and  go  to  2. 
k    o 

4.  Set  k  =  k+1,  go  to  2. 

The  algorithm  is  simply  a  restatement  of  the  procedure  developed 
in  Sections  4.5  and  4.6  and  therefore  determines  a  p-center  solution  for 
the  desired  value  p  . 

A  set  covering  problem  is  solved  for  each  value  of  t  ,  unless 
V(t  )  =  V(t   ).   The  maximum  number  of  set  covering  problems  needed 
is  m  where  m  =  f  J  +  m.   Thus  as  m  gets  large,  large  numbers  of  set 
covering  problems  are  generated  and  the  speed  of  the  algorithm  is 
largely  dependent  on  the  speed  of  solving  set  covering  problems. 

The  algorithm  was  programmed  for  an  IBM/65  at  the  University 
of  Florida  Computing  Center.   The  set  covering  algorithm  that  is  used 
in  the  program  is  a  cutting  plane  algorithm  developed  and  written  by 
Ratliff  T25]. 

Each  vector  is  stored  as  a  binary  word.   The  generation  of  the 
vectors  and  the  determination  of  dominating  vectors  is  done  using 
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binary  operations  on  words.   Also  some  computational  efficiency  is 

gained  by  generating  the  vectors  in  the  set  V(t,  . )  using  the  vectors 

k+1 

in  the  set  V(t,  ). 
k 

Examples  were  computed  using  the  programmed  algorithm  for 

m  =  5,  10,  15,  and  20  with  the  points  P.  sampled  from  random  numbers 

with  uniform  distribution  between  0  and  100,  and  the  constants  w 

i 

between  0  and  10;  but  with  g.  =  0,  for  1  <  i  <  m.   From  three  to  five 
examples  of  each  size  problem  were  computed.   The  cumulative  execution 
time  to  generate  all  the  set  covering  problems,  and  the  cumulative 
execution  time  to  solve  all  the  set  covering  problems  that  are  required 
to  obtain  a  p-center  solution  were  measured  and  recorded  for  each  p, 
2  <  p  <  m-1.   The  average  of  these  times  was  taken  over  the  example 
problems  for  each  m  and  are  presented  as  times  A  and  B,  respectively, 
in  Tables  4.1  and  4.2.   The  total  time  required  to  obtain  all  p-centers 
is  presented  as  time  C  . 

In  Tables  4.1  and  4.2,  notice  the  initial  time  required  to 
generate  the  first  set  covering  problem  for  each  m;  but  then  notice 
that  the  accumulation  of  time  to  generate  the  succeeding  set  covering 
problems  is  fairly  constant  and  a  small  proportion  of  time  used  to 
generate  the  initial  problem. 

For  a  given  value  of  m,  the  time  required  to  solve  each  set 
covering  problem  is  fairly  constant.   The  increase  in  the  accumulation 
of  time  used  in  solving  set  covering  problems  to  obtain  succeeding 
p-centers  is  due  to  the  increase  in  the  number  of  set  covering  problems 
requiring  a  solution. 
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These  computational  results  indicate  that  the  algorithm  does 
provide  a  computationally  feasible  method  to  solve  p-center  problems 
for  all  values  of  p  when  m  is  small. 

Table  4.3  gives  the  input  data  and  output  data  for  the 
Example  4.2.1.   Figures  4.3  and  4.4  illustrate  the  optimal  p-center 
and  the  corresponding  diamonds  for  p  =  2  and  3,  respectively. 

4.8  Additional  Observations 

The  algorithm  of  Section  4.7  may  be  initialized  for  any  value 
of  t  ,  say  t  ,  by  setting  k=r  and  computing  V(t  )  in  step  2.   Then 
to  obtain  the  p-center,  where  p   is  specified,  computational  savings 
may  possibly  be  gained,  particularly  when  p   is  small,  by  starting  the 
algorithm  at  a  larger  value  of  t  .   If  the  solution  to  the  set  covering 
problem  generated  by  t   is  p  and  p  ^  p    then  the  algorithm  would  con- 
tinue to  iterate  to  the  p, -solution.   If  p  <  p  ,  then  t   was  chosen  too 

1  1        r 

large,  and  the  algorithm  would  have  to  be  restarted  at  some  t   <  t  . 

k    m 

This  discussion  points  out  that  considerable  computational 

savings  could  be  gained  if  bounds  were  known  for  the  objective  function 

value  of  the  optimal  p-center.   Using  these  bounds,  the  procedure 

would  require  solutions  only  to  set  covering  problems  generated  by  t 

k 

that  are  within  these  bounds.   No  such  bounds  appear  to  be  available 
presently,  and  additional  research  is  needed  to  establish  them. 

Using  the  transformation  T,  a  decomposition  of  the  p-center 
problem  into  one-dimensional  p-center  subproblems  may  be  obtained. 
Although  the  subproblems  are  not  independent,  solutions  to  the  original 
problem  may  be  obtained  from  solutions  to  some  subproblems.   Currently, 
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TABLE  4.  1 

AVERAGE  EXECUTION  TIMES  FOR  THE  p-CE.\TER  PROBLEM 

m  =  20,  15 

m  =  20  m  =  15 

p       A  B  C         A  B  C 

19  20.074  .677 

18  20.441  1.273 

17  20.815  2.209 

16  21.151  2.677 

15  21.514  3.577 

14  21.854  4.108  4.530  .333 

13  22.120  4.601  4.802  .790 

12  22.564  5.489  5.082  1.438 

11  22.910  6.256  5.353  1.940 

10  23.365  8.067  5.633  2.696 

9  23.780  9.278  5.945  3.753 

8  24.262  11.278  6.238  4.605 

7  24.841  14.944  6.596  6.390 

6  25.418  17.802  6.890  7.023 

5  26.060  22.201  7.274  10.425 

4  26.848  29.254  7.601  12.449 

3  27.792  38.831  7.871  13.974 

2  28.808  52.457      81.265  8.156  16.835      24.991 


A:  average  cumulative  time  in  seconds  to  generate  all  ¥"(*,)  UP  to 
the  p-center  solution. 

B:  average  cumulative  time  in  seconds  to  solve  all  set  covering 
problems  up  to  the  p-center  solution. 

C:  total  average  time  in  seconds  required  to  obtain  all  p-center 
solutions. 
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TABLE  4.2 

AVERAGE  EXECUTION  TIMES  FOR  THE  p-CENTER  PROBLEM 

m  =  10,  5 


m  =  10 


m 


B 


B 


1.085 


244 


8 

1.185 

.462 

7 

1.271 

.752 

6 

1.387 

.984 

1.691 


1.908 


4  2.031 
3  2.345 
2     3.002 


2.997 
4.209 

6.849 


9.851 


.070 

.070 

.176 

.189 

.278 

.480 

.758 


A:   average  cumulative  time  in  seconds  to  generate  all  V(t  )  up  to 
the  p-center  solution. 


B:   average  cumulative  time  in  seconds  to  solve  all  set  covering 
problems  up  to  the  p-center  solution. 

C:   total  average  time  in  seconds  to  obtain  all  p-center  solutions. 
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TABLE  4.3 
DATA  FOR  EXAMPLE  4.2.1 


Input 

i 

a 
i 

b. 

l 

w 

l 

gi 

1 

-0.5 

4.0 

1.0 

0.0 

2 

-1.0 

3.0 

1.0 

0.0 

3 

0.0 

3.0 

1.0 

0.0 

4 

2.0 

3.5 

1.0 

0.0 

Output 

k 

X 

v(t*] 
k 

I 

< 

*<tj) 

*; 

1 

0.  177 

[0.1,0, 

,0] 

0.  177 

[0,0,1,1] 

0.354 

2 

0.354 

[0,1,0, 

,0] 

0.354 

[0,1,1,0] 

0.530 

3 

0.530 

[0,1,1 

,1] 

0.530 

[1,0,1,0] 

0.530 

4 

0.707 

[1.0,0, 

,1] 

0.530 

[0,1,1,1] 

0.  884 

5 

0.  884 

[1,1,0, 

1] 

0.884 

[1,1,1,0] 

1.061 

6 

1.237 

[1,1,1, 

,1] 

1.061 

[1,1,1,1] 

1.237 

p=3       V(tJ_)  =  {[0,0,0,1],  [0,1,1,0],  [1,0,0,0,]} 

jg  =  {p  =  (-0.5,4.0),  P4  =  (2.0,3.5),  X3  =  (-0.5,3.0)} 

P  =  2        V(t^)  -  {[0,0,0.1],  [1,1,1,0]} 

-2  =  {P4  =  (2-°.3-5),  X9  =  (-0.5,3.25)} 
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£  =  f?l'P4'X3^ 


->■  x 


Figure  4.3   Optimal  3-Center  for  ExamDle  4.2.1 
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x°  =  fp  x  1 

-2    '•4'  2  J 


Figure  4.4   Optimal  2-Center  for  Example  4.2.1 


101 


solutions  to  the  one-dimensional  problems  require  the  same  procedure 
and  computational  effort  as  the  original  problem,  that  is,  solving 
a  sequence  of  set  covering  problems  to  determine  the  minimum  number 
of  intersections  of  intervals  along  an  axis.   However,  the  one- 
dimensional  problems  appear  to  have  special  structure,  as  do  the  one- 
dimensional  minimax  location  problems  studied  in  the  previous  chapters. 
Computational  efficiencies  might  be  gained  if  the  subproblem  structures 
allowed  an  improved  solution  procedure  to  the  set  covering  problems. 


CHAPTER  5 


A  ONE  FACILITY  PROBLEM  WITH  RESPECT 
TO  A  COMPACT  SET 


5.1   Introduction 

In  this  chapter  an  extension  of  the  one  facility  location 
problem  studied  in  Chapter  2  is  considered.   The  problem  is  to  locate 
in  the  plane  one  new  facility  with  respect  to  a  compact  set  of  exist- 
ing facilities  so  as  to  minimize  the  maximum  cost,  where  costs  are 
linear  functions  of  the  rectilinear  distance  between  new  and  exist- 
ing facilities.   The  problem  also  includes  upper  bounds  on  the  recti- 
linear distance  from  the  new  facility  to  each  existing  facility  in 
the  compact  set. 

The  transformation  T  is  applied  to  this  problem  and  gives  a 
decomposition  into  independent  subproblems. 

Some  results  are  established  which  extend  the  work  of  Francis 

[13]  for  a  one  facility  minimax  location  problem  in  E   with  respect 

n 

to  a  composite  set  and  using  more  general  distances  then  the  recti- 
linear distance.   These  results  are  used  to  establish  properties  of 
the  subproblems,  which  are  analogous  to  those  properties  established 
in  Chapter  2,  and  to  obtain  a  solution  procedure  for  the  subproblems. 

The  optimal  solutions  to  the  subproblems  determine  the  optimal 
solutions  tc  the  original  problem. 
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5.2   Formulation 

Let  X  denote  a  point  in  E   to  be  determined  for  the  location  of 
the  new  facility.   Let  S  be  a  compact  set  in  E   representing  the  loca- 
tion of  existing  facilities.   For  a  point  P  =  (a,b)  in  S,  define  w(P)  to 
be  a  positive  continuous  function  on  S  and  define  g(P)  to  be  a  contin- 
uous function  on  S.   Define  a  function  f(X)  on  E  by 


f(«  =  res{w(p)Hx-p"i +g(p)) 


Since  w(P)  ,  g(P)  and  ||  X  -  P||   are  continuous  functions  of 
P  on  S,  the  term  inside  the  brackets  is  continuous.   Thus  the  maximum 
over  the  compact  set  S  exist  and  f(X)  is  well  defined. 

Let  c(P)  be  a  positive  continuous  function  on  S  that  is  an 
upper  bound  on  the  rectilinear  distance  between  X  and  P  for  each  P  in  S. 

The  one  facility  problem  to  be  studied  in  this  chapter  may  be 
written  as  follows: 

min   max 


|w(P)|j  X  -  p||   +  g(P)j  (5.2.1) 

s.t.      ||  X  -  p||   <  c(P)  ,   Pes  . 


XeE   Pes 


Notice  that  for  the  special  case  with  S  finite,  problem  (2.2.1) 
is  obtained.   Also,  problem  (5.2.1)  without  constraints,  referred  to  as 
problem  (5.2.2),  is  a  special  case  of  (5.2.1)  obtained  by  letting  c(P) 
be  infinite  on  S. 

5. 3  Applications 

The  one  facility  location  problem  (5.2.1)  may  find  applications 
in  contexts  of  locating  an  emergency  facility  or  a  facility  that  is 


104 

accessible  to  existing  facilities,  and  in  settings  where  the  rectilin- 
ear distance  approximates  travel  distance.   The  assumption  that  S  is 
compact  allows  consideration  of  applications  that  may  require  the  new 
facility  to  respond  to  an  emergency  at  any  point  in  some  set  or  region. 
For  example,  an  ambulance  center  might  be  expected  to  respond  to  emer- 
gencies that  could  occur  at  any  point  in  a  city  district,  not  just  at 
a  finite  number  of  points. 

The  function  w(P)  may  have  interpretations  analogous  to  those 
considered  in  the  previous  chapter,  for  example,  a  cost  function  per 
unit  of  distance;  the  function  g(P)  also  may  have  analogous  interpreta- 
tions as  a  cost  function  on  the  set  S. 

5.4  Related  Work 

The  literature  on  problems  of  locating  new  facilities  with 
respect  to  compact  sets  of  existing  facilities  is  not  as  extensive  as 
the  literature  on  location  problems  for  a  finite  set  of  existing  facil- 
ities.  Francis  [15]  considers  a  special  case  of  problem  (5.2.2)  with 
w(P)  =  1  and  g(P)  =  0;  this  work  appears  to  be  the  only  work  explic- 
itly using  rectilinear  distances.   Hearn  [19]  considers  a  one  facility 
minimax  problem  using  Euclidean  distances  with  respect  to  a  convex 
polyhedron  and  references  related  work  for  Euclidean  distances. 

The  developments  in  Section  5.5  are  extensions,  to  the  case 
with  S  compact  and  g(P)  *'  0 ,  of  corresponding  results  obtained  by 
Francis  |"13]  for  the  one  facility  minimax  problem  with  S  finite. 

5. 5  Properties  of  a  Minimax  Problem 

Results  are  established  in  this  section  that  will  be  used  in 
developing  a  solution  procedure  for  problem  (5.2.1),  but  the  results 
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also  hold  for  more  general  one  facility  minimax  problems.   Hence  for 
completeness,  the  results  are  established  for  the  general  case;  special 
cases  of  these  results  will  be  referenced  explicitly  when  they  are 
needed  in  later  sections. 

For  any  norm  defined  on  E   and  for  any  point  X  in  E  ,  let  ||  X|| 
denote  the  norm  of  X.   Then  a  distance  function,  relative  to  a  norm, 
may  be  defined  on  E   by  writing  the  distance  between  any  two  points  X 
and  Y  in  E   as  ||  X  -  Y||  .   It  can  be  shown  [12]   that  any  norm  is  con- 
vex, and  hence  continuous  on  E  ,  and  that  \\   X  -  y||  is  convex  and  con- 
tinuous on  E0  . 

To  formulate  the  problem  under  consideration,  let  (J  X||  be 

a  norm  on  E  ,  let  S  be  a  compact  set  in  E   and  let  w(P)  and  g(P)  be 
n  n 

functions  on  E  ,  with  w(P)  positive  and  continuous  on  S  and  g(P)  contin- 
uous on  S.   Then  the  problem  of  interest  may  be  written  as 


mi 
Xe 


in   max  f    ,  n      n        "1 


Notice  that  if  the  distance  in  (5.5.1)  is  the  rectilinear 
distance  and  if  n  =  2,  problem  (5.2.2)  is  obtained. 
Define  the  function  F(X)  by 


F(X)    »   J£{w<P)||x-p||     +   g(P)} 


The  function  F(X)  is  well  defined  by  an  argument  similar  to  that  in 
Section  5.2. 

In  each  of  the  following  properties,  the  special  case  obtained 
by  assuming  S  is  finite  and  g(P)  -  0  on  S  corresponds  to  a  similar  prop- 
erty obtained  by  Francis  [12]. 
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First  a  lower  bound  for  F(X)  is  obtained.   Define  a  real  valued 

function,  B(P  ,P-),  for  all  (P  ,P„)  in  S  x  S  C  E„   by 

a     P  a     :  2n 

B(P  ,PS)  =  max  {__!_-—  W(p  )w(Pft)  ||  P   -  pj| 
a    p       U(P  )+w(p.)  L  a         P  "  a        p" 

a  z 

+    w(P  )g(P  )  +  w(P  )g(P  )   ,  g(P  ),  g(P-)}  .    (5.5.2) 

Then   B(P  , P~ )  is  continuous  on  S  x  S,  since  each  component  is  contin- 

a     f 

uous  on  S  x  S.   Also  S  x  S  is  compact  in  E,,   since  S  is  compact  in  E  . 

2n  n 

Therefore   there   exist    a   point,    say    (P    ,P    )    in   S    x  S   and   a   number,    B      , 

st  st 

defined   by 

Bst=    (P    ,pa)cS.S   {B(IW}   =    B(Ps'Pt)    • 
a     - 

Property    5.5.1.       For    all    X   in   E    ,    B       <  F(X). 
n        st 

Property    5.5.2.       B        =    F(X)    if    and   only    if 

(i)   ||p.  -Ptll-  II  p.  -x||+  ||x-Pt||. 

(2)    w(P    )||  X    -   P    ||     +    g(P    )    =    w(P    )  \\  X    -   P    ||     +   g(P    ), 

9  O  O  b  L  I* 

and      (3)  w(P)  ||  X  -  pj|  +  g(P)  <  B    for  all  P  in  S,  P  i  P  ,P  . 

The  proofs  of  Properties  5.5.1  and  5.5.2  are  direct  extensions 
of  the  proofs  of  the  analogous  results  in  [12]. 


Property  5.5.3.   There  exists  an  X   in  E   such  that  F(X  )  <  F(X) 
on  o 

for  all  X  in  E  . 

n 
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Proof:   For  each  P  in  S,  w(P)  \\   X  —  P  |(  +  g(P)  is  a  convex  function  in 

X  on  E  .   Since  the  maximum  of  convex  functions  is  convex, 
n 


max 
pes 


|w(P)  ||  X-p||  +  g(P)}  =  F(X) 


is  a  convex  function  on  E  .   Since  F(X)  is  a  convex  function  on  E  , 

n  n 

and  since  E   is  an  open,  convex  set,  F(X)  is  continuous  on  E  . 
n  n 

Therefore  the  set 

K  =  {X  e  E   :  F(X)  <  k] 
n  J 

for  any  real  k  is  a  closed  set  on  E  .   K  is  also  bounded  so  that  it  is 

compact.   Thus,  F(X)  attains  a  minimum  value  for  some  X  in  K,  and, 

as  shown  in  [12],  there  exists  an  X   in  E   such  that  F(X  )  <  F(X)  for 

on  o 

all  X  in  E  . 

n 

If  X   is  a  minimum  point  of  F(X) ,  then  by  the  definition  of 
o 

B   ,  F(X  )  s  g(P)  for  all  P  in  S.   The  following  property  concerns  the 

case  when  F(X  )  =  g(P  )  for  some  P   in  S. 
os  s 

Property  5.5.4.   Suppose  that  X   is  any  point  in  E   such  that  F(X  )  =g(P  ) 

o  n  os 

for   some  point   P      in  S.      Then  X     =   P      is   the  unique  minimum  value  of   F(X). 
s  os 

Proof:       Suppose  F(X   )    =   g(P    )    for   some   P      in   S,    then 

os  s 

g(P    )    =    F(X   ) 
s  o 


max 

Pes 


{w(P)J|  Xq   -  P||     +   g(P)} 
>  w(P    )||  X      -   P J|     +    g(P    ) 
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implies  0  ^  ||  x  -  p  |j  so  that  X  =  P  .   Now  suppose  X  is  any  point 


o    s 

in  E  with  X  ^  P  ,  then 
n  s 

F(X)  ^  w(Pg)||  X  -  pj|  +  g(Pg) 

>  g(P  )  =  F(P  )  . 

s       s 

Therefore  P   is  the  unique  minimum  of  F(X). 
s 

Under  the  assumption  that  F(X  )  =  g(P  ),  it  is  possible  that 

o        & 

there   is   only   one  point   in   S,    namely,    P    ,    such   that 

F(X   )    =   w(P   )||  X      -  P    ||     +   g(P   ). 

O  SOS  s 

However,  if  F(X  )  >  g(P)  for  all  P  in  S,  then  there  are  at  least  two 
o 

points  in  S,  namely,  P   and  P  ,  such  that 

F(X  )  =  w(P  )||  Xq  -  P  ||  +  g(Pg)  =  w(Pt)||  Xq  -  Pt||  +  g(Pt). 

This  assertion  is  stated  as  a  necessary  condition  for  X  to  be  a  min- 
imum in  the  following  property. 


Property  5.5.5.   Suppose  S  contains  at  least  two  points;  let  X  be 

such  that  F(X  )  <  F(X)  for  all  X  in  En ,  and  suppose  F(X  )  >  g(P)  for 
o  2  o 

all  P  in  S.   Then  there  exist  distinct  points  P  and  P   such  that 

F(Xq)  =  w(Pg)||  Xq  -  Pg||  +  g(Pg)  =  w(Pt)||  Xq  -  Pt||  +  g(Pt). 


Proof:   By  definition  of  F(X) ,  there  exists  some  point  P  in  S  such 
that 

F(X  )  =  w(P  )||  X   -  P  ||  +  g(Pg)  ^  w(P)||  Xq  -  P,'|  +  g(P) 

for  all  P  in  S,  P  4   P  .   Assume  that  the  assertion  is  not  true,  that  is, 
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(P  )J|  X   -  P  ||  +  g(P  )  >  w(P)  ||  X   -  p|J   +  g(P) 


w 
'  s"  o    s 


for  all  P  in  S,  P  4   P  • 

s 


Part  I  of  the  proof  shows  that  under  the  above  assumption  there 
exists  a  neighborhood,  N(X  , a) ,  such  that 

w(P)  ||  X  -  P  ||  +  g(P  )  >  w(P)  \\   X  -  p||  +  g(P) 

for  all  X  in  N(X  ,cr)  and  for  all  P  in  S  with  P  ^  P  .   Part  II  shows 
o  s 

that  there  exists  an  X'  in  N(X  , a)  such  that  F(X')  <  F(X  )  which 

o  o 

contradicts  the  assumption  and  proves  the  theorem. 
Part  I:   For  each  P  in  S,  with  P  4   P  ,  define  6(p)  by 

S 

6(P)  =  w(P  )  \\   X  -P  ||  +  g(P  )  -w(P)||  X  -P||  -  g(P).      (5.5.3) 

Then  6  (p)  is  positive  for  every  P  j£   P   in  S  by  assumption. 
For  every  P  jL   P  in  S,  define 

£1(P)  =  4||  XQ6-PP||  +    1  (5'5-4) 

so  that  e  (P)  is  positive  and 

ifl-  ei(P)||xo-p||  >0. 

Then,  since  w(P)  is  continuous  on  S,  for  each  P  ^  P   in  S,  there 
exists  a  p  (P)  >  0  such  that 


|w(p)  -  w(p  )|  <  e  rp) 
a  1 


for  all  P   in  S  1*1  N(P,0  (P) ) . 
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Also  for  every  P  d   P  in  S,  define 

e2(p)  -§(nr-  Vp)l'xo-pl'  >  (5-5-5) 


so  that  e  (P)  is  positive  and 


6(P) 


-  c  (P)||  x  -  p||  -  e0(P)  >  o 


Then,  since  g(P)  is  continuous  on  S,  for  each  P  ^  P   in  S  there  exists 
a  p  (P)  >  0  such  that 

|g(P)  -  g(PQ)|  <  e2(P) 

for  all  P   in  S   N(P,p0(P)). 
a  2 

Since  w(P)  is  continuous  on  the  compact  set  S,  there  exists 
a  real  number  M  such  that  w(P)  <  M  for  all  P  in  S. 
Define  p (P)  by 

p  (P)  =  2^1  min  {P1(P)  '  D2(P)  '  6(P)/2 

-  gi(p)||  Xq-p||  -e2(P)}  •       (5.5.6) 

Then  p (P)  is  positive  and 

6(P)/2  -  e,  (P)||  x  -p||  -e  (p)-md(P)  >  0. 
1      o        ^ 

It  also  follows  from  the  definition  of  p (P)  that 

|w(P)  -  w(P  ) I  <  e  (P)  for  all  P   in  S  n  N(P,p(P)),  and  that 

a  1  q 

|g(P)  -  g(P  )|  <  e.(P)  for  all  P   in  S  ^  N(P,p(P)). 
'  °         or  '     2  a 

Now  the  collection  of  neighborhoods  {N(P,p(P))}  for  P  *  Pg 
is  an  open  cover  for  S,  and  since  S  is  compact,  there  exists  a  finite 
open  cover  of  S,  denoted  as  {N(P . ,p (P. ) ) }  for  1  <  i  <  n.   Then,  for 


Ill 


each  i,  the  terms  5 (P.),  e, (P.),  eo(P.)  and  o(P.)  are  defined  by 

i  1     i  2     i  l 

(5.5.3),     (5.5.4),    (5.5.5),    and    (5.5.6),    respectively.       Furthermore, 
(5.5.6)    implies   that 


6(P.)/2-  e   (P.)  II  X    -P.  ||    -  e    (P.)  -Mc  (P.)   >  0 

X  1         X  O  1  <bl  X 


for  1  <  i  <  n.   For  each  i,  define  the  function  f . (X)  by 

l 


f.(X)    =   w(Ps)||  X-Pjl    +g(Pg)  -w(P.)||  Xo-P.||    -g(P.). 


Then     f    (X)    is  continuous,    and  f.(X  )   =   6(P.)    for   each   i.      Therefore, 

1  1       o  1 

for  each  i,  there  exists  a  a!  >  0  such  that  f.(X)  >  6(P.)/2  >  0  for 

all  X  in  N(X  ,a\). 
o   1 

Define  the  positive  number   a.    by 

l 


a.    =  —J—  min  ia\  ,    6(P.)/2  -  e    (P.)||  X    -P.| 
l        2M+1  L   l  l  1      l    "     o        l 


-    e2(Pi)    ~  Mp   (P.)j    .  (5.5.7) 

Define  P.    for   1  <  i  <  n  by 
i 

B.   =  6(p.)/2  -  e   (p.)||  x    -P.  ||    -  e_(P.)   -  Mo  (P.)    -Ma., 
i  i  1     i   "     o       i"         2     i  i  i 

and  note  that    (5.5.7)    implies    that   0.    >  0   for   each    i.       Then   it    follows 

that   f.(X)    2:  5(P    )/2  >  0   for   all    X   in  N(X    ,a.). 
11  o       1 

Now  for   all   P   in  N(P.,p(P.)),    for   all   X  in  N(X    ,c),    and  for 

11  o       1 

each   i,    1  <   i  <  m,    it   follows   that 

w(P    )J|  X-P    ||    +g(P    )  -w(P)||  X-P||      -   g(P) 
s  s  s 

*w(P)||x-P    !|    +g(Po)  -w(P)  \\  X-x\\    -w(P)||x    -P.  ||    -w(P)||  P.  -p||  -g(P) 

5  S  S  O  OX  1 

(by    the   triangle   inequality) 
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>   w(P    )||  X-P  J|  +    g(P)    -   M0.  -   Mp   (P.)  -w(P)||  X     -P.  |!    -g(P) 


o        1 


(by   the  definitions   of   M,    a.    and  p(P.)) 

11 


=   w(P)  ||  X-P    ||  +   g(P   )  -w(P   )[|  X     -P.  ||  -g(P.)  +w(P.)||  X     -P.||  +   g(P.) 

S  o  S  X  O  1  1  1  O  1  1 

-   w(P)||  X     -P.  ||  -g(P)  -Ma.  -  Mp   (P.) 


o         1 


(by   subtracting    and   adding  w(P.)      X    -P.      +   g(P.)) 

l   "     o       i"  i 

=    f .  (X)  +  (w(P.)    -   w(P))||  X     -P.  J|    +g(P.)  -g(P)  -Ma.  -   Mp   (P.) 


o  1 


(by   the   definition  of   f.(X)) 

l 


;>  6(p.)/2  -  e,(P.)||  X  -P. ||  -eo(p.)  -Ma.  -  Mo  (P.) 
l       1   l  '   o    l     2   i      i        i 


(by  the  definitions  of  6 (P  ) ,  e  (P  )  and  e  (P  )) 

l     1   l       2   i 


£  g   >  0  . 

l 

-i         min         .     a    min  ._  . 

Now  let  a  =  ,_ (a.)  and  let  S  =  „„.^  (3.) 

l<i<n   l  l<i<n   l 

Then  for  all  X  in  N(X  ,a)  and  for  all  P  *  P   in  S, 

o  s 

w(P    )  ||  X-P    ||  +    g(P    )    -   w(P)  ||  X-P||  -   g(P)    s   3   >   0. 

S  S  S 

This   proves   Part    I. 


Part   II:       It   will   be   shown   that    there   exist    some  X'    in  N(X    ,a)    with 

o 

X'    4   X      amd  X'    ^   P      such    that 
o  s 

w(P   ){|  X     -P    ||  +g(P    )   >  w(P   )  \\  X'  -P    II    +g(P   ). 

S  O  S  5  S  S  S 

Let 

X  =  -  mm   ] ,    If- 
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so   that    \\\  X     -  P   jj       <  o     and      X  <   1,    0  <  1   -    X  <  1. 

O  S       *- 

Choose  X'    =    (1  -  X)X     +    XP    . 

o  s 

Since      f(X   )   >  g(P   ),      w(P   )\\  X     -  P    J|  +   g(Pg)    >  g(Pg)      so   that   Xq  4  Pg. 

Then        X'    4  P        and     X'    4  X   .      Furthermore, 
s  o 

!|x'-Xo||2    =    l|XPs-XXoH2    =    x||Ps-XoJ|2    <o 

so   that  X'    is   in  N(X    ,  a).      Also 

o 

||  X'  -Psl|=   ||  (1-«)XQ  -  (l-c)P8||  =   (l-e)||  Xo-Ps||<  ||  Xo-pJ|  . 

Thus 

w(P    )||  X'  -P    \\  +   g(P   )    <  w(P   )||  X     -P    ||  +   g(P   ).  (5.5.8) 

s  s  s  S  0        s  s 

The  result   for  Part   I    implies 

w(P    )||  X'   -P    ||  +   g(P    )    =    ™*  {w(P)||  X'   -P||  +    g(P)}   =    F(X'). 

S  S  S  r  Go      \.  J 

Therefore   (5.5.8)  implies  F(X' )  <  F(X  )  and  is  a  contradiction. 

o 

This  proves  Property  5.5.4. 

5. 6   The  Transformation  and  Decomposition 

The  transformation  T  is  applied  to  problem  (5.2.1)  to  obtain 
a  rotated  problem  and  a  decomposition  that  is  analogous  to  the  develop- 
ments in  Section  2.5.   Let  S'  be  the  image  of  the  set  S  under  T.   The 
following  functions  are  defined  on  S' :   W(P' )  =  w(T~  (P')), 

G(P')  =  —  g(T_1(P')),  and   C(P* )  =  -i-  c(T_1(P')).   Then  for  a  point 

a/2-  v/2~ 

P  in  S  and  P'  =  T(P)  =  (a',b')  in  S'  ,  it  follows  that  W(P' )  =  w(P) , 
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G(P')  =  — —  g(P)  and  C(P' )  =  c(P).   Also   each  of  these  functions 

is  continuous  since  each  is  the  composition  of  the  continuous  functions 
T    and  w(P) ,  g(P)  or  c(P),  respectively.   With  these  definitions,  the 


rotated  problem  is  written  as: 


min 
X'  CE2 


"  {w<p,>H x'  -p'IL  +  g(p')} 


(5.6.1) 
s.t.  \\  X*  -P'  \\a  <  C(P'),      P'    e  S'  . 

The  equivalence  of  problems  (5.2.1)  and  (5.6.1)   is  established  by  the 
following  property  which  is  analogous  to  Property  2.5.2. 

Property  5.6.1.   The  point  X   is  an  optimal  solution  to  problem  (5.2.1) 

if  and  only  if  T(X  )  is  an  optimal  solution  to  problem  (5.6.1).   If  z 

and  z'  are  the  optimal  objective  function  values  for  (5.2.1)  and  (5.6.1), 
o 

respectively,  then  z   =  z'  J2~   . 

o    o 

Two  subproblems  are  now  defined  from  the  rotated  problem. 
Recall  that  a  point  in  S'  is  written  P'  =(a',b'). 


mi 
x 


■\\    p™  {■*>!«• -'I  ♦«<»•>} 

s.t.  Jx'  -a'  |  <  C(P')  ,   P'  e  S'  , 


(5.6.2) 


s.t.  |y'  -b'  1  <  C(P')  ,   P'  e  S' 


(5.6.3) 
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Each  subproblem  is  a  location  problem  in  one  dimension,  but 

with  respect  to  a  compact  set  in  E  .   The  following  property  relates 
optimal  solutions  of  the  subproblems  to  optimal  solutions  to  the 
rotated  problem  (5.6.1). 

Property  5.6.2.   If  the  point  x'  is  an  optimal  solution  to  problem  (5.6.2) 

o 

with  optimal  objective  function  value  z'  and  if  y'  is  an  optimal  solu- 
tion to  problem  (5.6.3)  with  optimal  objective  function  value  z'  then 
(x',y')is  an  optimal  solution  to  problem  (5.6.1)  with  optimal  objective 

function  value  z'  =  max  (z' ,z'). 
o         1'  2 

Thus  the  optimal  solutions  to  the  one  facility  problem  (5.2.1) 
are  obtained  from  the  independent  solutions  of  the  subproblems  (5.6.2) 
and  (5.6.3). 

The  constraints  of  each  of  the  problems  (5.6.2)  and  (5.6.3)  may 
be  simplified.   For  problem  (5.6.2)  define 

CI  =  P^  {*'-C<P'>} 
Then  problem  (5.6.2)  is  equivalent  to 


mm   max 
x'  eE  P'  es 


,  |w(P')|a'  -  x' |  +  G(P)j 


(5.6.4) 


s.  t.  C'  <  x'  <  C*  . 

J.  ^ 


Similarly,  for  problem  (5.6.3)  define 
so  that  problem  (5.6.3)  is  equivalent  to 
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"ce   Pmcs-  {w(p')'b'  -y'l  +  G(p,>} 


mi 

y 


1 

(5.6.5) 


are: 


min   max 
x'  eE   P'  e  S' 

mil 

y' 


s.t.         D^  <  y'  <  D'2. 

Problems  (5.6.4)  and  (5.6.5)  have  identical  structure  so  that  it  suffices 
to  develop  properties  and  solution  procedures  only  for  problem  (5.6.4). 

5. 7   Properties  of  the  Subproblems 

The  subproblems  (5.6.4)  and  (5.6.5)  written  without  constraints 

(w(P')la'  -  x'  |  +  G(P' )j  ,  (5.7.1) 

™n   P™es-    {w<p*>'b'  -  y'l  +  G(p,)}  •  (5-7-2) 

Analogous  to  the  development  in  Chapter  2,  the  optimal  solutions  to  the 
constrained  subproblems  may  be  obtained  from  the  optimal  solutions  to 
the  unconstrained  subproblems.   Since  problems  (5.7.1)  and  (5.7.2)  have 
identical  structure,  only  problem  (5.7.1)  is  considered  explicitly. 

Properties  are  now  developed  for  problem  (5.7.1).   These  prop- 
erties are  analogous  to  those  properties  presented  in  Section  5.5  for 
problem  (5.5.1).   However,  problem  (5.7.1)  cannot  be  considered  as 
a  special  case  of  problem  (5.5.1),  since  (5.7.1)  is  a  problem  in  one 
dimension  with  respect  to  a  set  S'  in  E9.   The  properties  of  (5.7.1) 
are  presented  below  as  remarks;  they  may  be  proved  in  an  analogous 
fashion  to  the  corresponding  properties  in  Section  5.5. 

Define  the  function  f(x')  on  E  by 

f(-V)  =  P'l€S'  {W(P,)IX'  "  3'l  +  G(P,)}  • 
It  is  assumed  that  there  exist  two  points  (a' ,b')  and  (a^.b^)  in  S' 
with  a'  -  <*l  ;    otherwise  x'  =  a'  is  the  optimal  solution  to  (5.7.1)  with 

W>1 1~  0(P.). 
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For      P1    =    (a*  ,b')    and  Plfc'.d')    in   S1    define 
Of  P 

*%>*?  -  -x  {w(P-)+w(Pi)[w(p;)w(pp)|a'  -c,i  +  ^p;)g(p3> 

a  P 

+  W(PI)G(P')1  ,  G(P'),  G(P')} 

p    a  J     a      P  J 

and   brt.b<P-.P->.    .  _^esxs -{blP;.»J>}  •  (5-7.3) 

a  P 


Remark  5.7. 1.   For  all  x'  inE.,  b   <g(x'). 
.  1   st 

Remark  5.7.2.   b  ^  =  f(x')  if  and  only  if  (1)  | a'  -  c ' |  =  | a'  -  x' | 

st  sis 

+  |x'  -  c'l,  (2)  W(P')|x'  -  a'|  +  G(P')  =  W(P')|x'  -  c'|  +  G(P' ) , 
I        t  s  s         s         t  t'        t 

(3)  W(P')|x'  -  a' |  +  G(P')  <  b  .  for  all  P'  in  S' ,  P'  4   Pi, Pi- 

St  S    L 

Remark  5.7.3.   There  exists  an  x'   in  E„  such  that  f(x')  <  f(x') 

o     1  o 

for  all  x'  in  E  . 

Remark  5.7.4.   Let  x'  be  such  that  f(x')  <  f(x)  for  all  x'  in  E  , 

o  o  J- 

and  f(x')  >  G(P')  for  all  P'  in  S' .  Then  there  exist   distinct  points 
o 

P'    and   P'    such   that 
s  t 

f (x' )  =  W(P'  ) |x'   -  a'  |  +  G(P' )  =  W(P' ) |x'   -  c'  |  +  G(P' )    . 
o  sos  s  tot  t 

Additional  properties  are  now  developed  for  problem  (5.7.1). 

Analogous  to  the  unconstrained  one  facility  subproblem  in  Chapter  2,  the 

lower  bound  b    is  always  attained  at  optimality  for  problem  (5.7.1). 
st 

Property  5.7.1.   Let  b   be  the  lower  bound  defined  by  (5.7.3)  and  let 

1 J pq 

z'  be  the  optimal  objective  function  value  of  (5.7.1).   Then  z'  =  b 
o         P         J  o    pq 
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Proof:   Suppose  x'  is  an  optimal  solution  to  problem  (5.7.1).   If 

f(x')  =  z'  =  G(P')  for  some  P'  in  S* ,  then  by  Remark  5.7.1 
o     o      s  s  J 

b   <  G(P'),  but  by  taking  s  =  t,  b   =  G(P').   Since  G(P' )  Sb   =  G(P' ) , 
pq      s  ss      s  s    ss     s  ' 

b   =  b   . 
pq    ss 

If   f(x')  >  G(P')  for  all  P'  in  S* ,  then  by  Remark  5.7.4,  there 

exist  distinct  points  P'  and  P'  in  S'  such  that 

s  t 

z'  =  Hr(P')|x'   -a'  |  +  G(P')  =  W(P')  Ix'   -cl|  +  G(P') 
o  sos  s  tot  t 

which    implies    that    z'     is   determined   by   one  of   four  possibilities. 

(1)  x'  -a'  =  (z'  -G(P'))/W(P')      and      x*  -  c!  =  (zf  -  G(P'))/W(P! ) 

oso  s  s  otot  t 

(2)  x'  -  a'  =  (z'   -G(P'))    W(P')      and   -x'  +  c'  =  (z'  -G(P'))  'W(P') 

osos  s  otot  t 

(3)  -x'  +  a'  =  (z'   -G(P*  ))    W(P')      and      x'   -  c'  =  (z'  -  G(P'  ) )  /W(P'  ) 

osos  s  otot  t 

(4)  -x'  +  a'  =  (z'   -G(P'  ))/W(P'  )      and   -x'  +  c'  =  (z'   -G(P'  ))/W(P!  )  . 

oso  s  s  otot  t 

In   each    of    the   four   cases    eliminate  x1    and   solve   for   z' .       Cases    (1) 

o  o 

and  (4)  imply  that  z'  is  either  undefined,  with  denominator  W(P' )  -  W(P' ) , 

o  st 

or  gives  a  value  of  z'  greater  than  the  value  resulting  from  cases  (2) 

o 

and  (3).   Since  z'  is  the  minimum  value,  only  cases  (2)  and  (3)  need 
be  considered.   Case  (2)  or  (3)  yields 

,;-w<P')iw<P')[w(p:)w(pl)  (±1)  (a;-c;)+w(Ps)G(P;)  +  w(p;)G(P;)], 

s     t 

respectively.   Since  |a'  -  c'   =  max  (±(a*  -c')), 

s    t  s    t 

z;-w(P')!w(P-)  [w(p;)w(p;> |a;  -  c;|  +  w(p;)o(p;)  *  w(pyG<p;>]  . 

s      t 

(5.7.4) 
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Recall  that  z1  >  G(P' )  for  all  P'  in  S'  so  that  b    is  the  maximum  of 
o  pq 

expressions  like  the  right-hand  side  of  (5.7.4);  therefore,  b   ^  z' , 

pq     o 

but  Remark  5.7.1  gives  b   <  z'  so  that  b   =  z' . 

pq     o  pq     o 

The  next  property  of  problem  (5.7.1)  corresponds  to  a  similar 
result  for  problem  (2.6.1)  in  Chapter  2. 


Property  5.7.2.   Problem  (5.7.1)  has  a  unique  optimal  solution. 


Proof:   Suppose  problem  (5.7.1)  has  optimal  objective  function  value 

z' .   If  z'  =  G(P' )  for  some  P'  in  S'  ,  then  the  discussion  of  this  case 
o       o 

in  Section  5.5  shows  that  problem  (5.7.1)  has  a  unique  optimal  solution. 

If  z'  >  G(P' )  for  all  P'  in  S' ,  suppose  that  x'  and  x'  are  both 
O  12 

optimal  to  (5.7.1)  with  x'  4   x'    Then,  by  Remark  5.7.4,  there  exist 
P'  and  P'  such  that 

J.  «-> 

Zo=  W(Pi)ixi  -a[l  +G<Pi)  =  ^pl3^  ~a'?\  +  G(P2)" 
Remark    5.7.4    also   allows    the   assumption   that   P'    4   P~-       Since   f(x')    is 

J.      £ 

a  convex  function,  x'  =  \  x'  +  (I-a)x'  for  0  <  \  <  1,  is  also  an  optimal 

solution  to  problem  (5.7.1).   By  the  definition  of  f(x'),  there  exists 

some  point  P'  in  S'  such  that 
o 

z'  =  W(P')|x'  -a'J  +  G(P') 
o       3        3        3 

which    implies 

Zo  =  W(P3)  I  X(a3  "  Xi}  +  (1~X)  (S3  "  X2}  I    +   G(P3} 

<    >j  W(P^)|a^-x^J  +G(P^)    |+  (l-\)^W(P^)|a^-x;|  +  G(P^)J   .        (5.7.5) 
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Now  observe  that  by  definition  of  f(x')  and  f(x') 

J.  w 


W(Pp|a^-x'|  +G(Pp  <  W(Pp|a^-x^|  +G(Pp    =    z^  (5.7.6) 


and 


W(P3)'a3  ~X2^  +G(P3}   -    WCP25  I  a2  "^   +G(P2}  =  Zo    '  (5.7.7) 

If    strict    inequality   holds    in   either    (5.7.6)    or    (5.7.7),    then 
substituting    into   expression    (5.7.5)    gives 

z'    <   \z'     +    (l-X)z'    =    z' 
o  o  o  o 

so  that  equality  must  hold  in  both  (5.7.6)  and  (5.7.7).    But  equality 
in  both  expressions  implies  that 

W(P3)'a3  "Xll  +  G(P3}  =  w(p3>la3-x'  I  +  G<p3) 
and 

W(P3}la3  ~*fcl  +  G(P3)  =  W(P3}'a3  "X'  '  +G(P3)- 

These  equations  imply  that  either  x  =  x9 ,  a  contradiction  to 
the  assumption,  or  that  x'  =  x'  or  x'  =  x'   a  contradiction,  since 
0  <  \  <  1.   Therefore  problem  (5.7.1)  has  a  unique  optimal  solution. 

The  next  result  gives  an  explicit  expression  for  the  unique 
optimal  solution  to  problem  (5.7.1)  in  ternis  of  the  bound  b   . 

Property  5.7.3.   Suppose  b  ^  =  b(P' ,P!).   If  b  .  =  G(P' )  define 

E : St        St  St        S 

X*   =  a'  ;  if  b   =  G(P'  )  ,  define  x'   =  a'  ,  and  if  b   >  G(P'  )   and 
St     s      st       t  st     t  st      s 

b  _  >   G(P! ) ,  and  for  a' <  a'  ,  define 
st      t  s    t 
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,-^i^  [w(p;)a;  +  w(P')a;  +  g(p;)  -  g<p;)] 


Xst   W(P  , 

s     t 

Then   in   any   of  the  above  cases,  x'   is  the  optimal  solution  to 

st 

(5.6.1). 

Proof:   Suppose  b  .  =  G(P').   Then  by  Property  5.7.1,  if  x'  is  an 
st      s  o 

optimal  solution  to  (5.6.1),  f(x' )  =  b  .  =  G(P' )  and  by  Property  5.5.2, 

o     st      s 

x'  =  a'  is  the  unique  optimal  solution  to  (5.6.1).   Hence  x'   =  a'  is 
os  st     s 

the  unique  optimal  solution  to  (5.7.1). 

Suppose  b   =  G(P' ) ,  then  by  a  similar  argument,  x'   =  a'  is 
st      t  st    t 

the  unique  optimal  solution  to  (5.7.1). 

Suppose  b  >  G(P' )  and  b   >  G(P').   Then  a  development  anal- 
st     s       st      t 

ogous  to  the  proof  of  Lemma  4.6.1  shows  that 

W(P*)|a'  -x'J  +G(P')  =  W(P')|a!  -x'  |  +  G(P' )  =  b  ., 
s  '  s   st      s      t  '  t   st'      t     st 

and  that  x'   is  the  unique  point  satisfying  these  equalities.   Also, 

analogous  to  the  proof  of  Lemma  4.6.1,  a'  <  x'   <  a' 

s    st    t 

It  remains  to  be  shown  that  f(x'  )  =  b   .   For  any  other  point, 

st     st 

say  P1  in  S' ,  b  .  ^  b   and  b  .  £  b   ,  since  b   is  the  maximum  of  all 
r         st    sr      st    rt         st 

such  expressions.   Then  a  proof,  analogous  to  the  proof  of  Property  4.5.1 

shows  that 

W(P')|x*   -  a' I  +  G(P')  <  b  ^. 
r  '  st    r'      r     st 

Therefore  the  conditions  of  Remark  5.7.2  are  satisfied  so  that 

f(x'  )  =  b  .  and  x'   is  the  optimal  solution  to  (5.7.1). 
st     st      st 

The  preceding  remarks  and  properties  in  this  section  indicate 
a  procedure  for  determining  the  unique  optimal  solution  to  problem  (5.7.1) 
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Recall  however,  that  it  is  desired  to  obtain  an  optimal  solution  to 

the  subproblem  (5.6.2),  which  is  equivalent  to  problem  (5.6.4). 

If  x1  is  the  optimal  solution  to  problem  (5.7.1)  and  if  x'  is  feasible 
o  o 

to  (5.6.4),  then  x'  is  the  optimal  solution  to  (5.6.4).   The  case  where 

o 

x'  is  inf easible  to  (5.6.4)  is  considered  in  the  next  property. 

Property  5.7.4.   If  x1  is  an  optimal  solution  to  (5.7.1)  and  x' 
o  o 

is  infeasible  to  (5.6.4),  then  the  optimal  solution  to  (5.6.4)  is 
either  C'  or  C'    Specifically,  if  x'  <  C'  ,  C'  is  the  unique  optimal 

1       £  Oil 

solution  to  (5.6.4)  or  if  x'  >  C '  ,  C'  is  the  unique  optimal  solution 

O     Z  2. 

to  (5.6.4). 

The  proof  of  Property  5.7.4  follows  directly  from  the  proof  of 
Property  2. 6. 1. 

5.8   Solution  Procedure 

The  properties  developed  in  the  preceding  section  for  the 
subproblem  (5.5.4)  and  the  analogous  properties  that  exist  for  sub- 
problem  (5.5.5)  are  used  to  determine  a  solution  procedure  for  both 
subproblems.   A  procedure  is  also  given  for  finding  all  the  alternative 
solutions  to  the  one  facility  problem  in  two  dimensions,  (5.2.1),  using 
the  solutions  to  the  subproblems.   This  procedure  is  completely  anal- 
ogous to  the  procedure  for  the  one  facility  problem  in  Chapter  2. 

Subproblem  (5.6.2)  is  solved  by  computing  the  lower  bound  from 

expression  (5.7.3)  and  determining  the  value  x'  using  Property  5.7.3. 

If  x'  is  feasible  to  problem  (5.6.4),  then  x'  is  the  unique  optimal 
o  o 

solution  and  b   =  z'  is  the  optimal  objective  function  value. 
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If  x'  is  not  feasible,  the  unique  optimal  solution  is  x'  =  C'  or 
o  o     l 

x'  =  C'  as  determined  by  Property  5.7.4,  and 
o    2 

Z'l  =   P^eS'  {W(p,)'a'  "  Cj'  +  G(p,)} 


for  j  =  1  or  2 ,  respectively. 

Similarly,  subproblem  (5.6.3)  is  solved  by  computing  the  lower 

bound  b    from  expression  (5.7.3)  but  with  b'  and  d'  replacing  a'  and  c' , 
st 

respectively.   Then  y'  is  determined  from  b   using  analogous  results 

to  those  given  in  Property  5.7.3.   If  y'  is  feasible  to  (5.6.3),  then 

it  is  the  unique  optimal  solution  and  b   =  z'  is  the  optimal  objective 

function  value;  otherwise,  y'  =  D'  or  y'  =  D'  and 

o     1     o    z 


Z2 


=  p^s'  {w(p,)'D;  -  b'i +  G(p,)} 


for  j  =  1  or  2 ,  respectively. 

By  Property  5.6.2,  the  optimal  objective  function  value,  z' , 

to  the  rotated  problem  (5.6.1)  is  given  by  z'  =  max  (z'  z')  where 

o         1   2 

z'  and  z'  are  the  optimal  objective  function  values  to  (5.6.2)  and 

(5.6.3),  respectively. 

The  point  fr' ,y')  is  one  optimal  solution  to  problem  (5.6.1)  by 
o   o 

Property  5.6.2.   Three  cases  need  to  be  considered  to  determine  all 

alternative  solutions. 

Case  1:   z'  =  z'  =  z' .   Then  (x' ,y')  is  the  unique  optimal  solution 
o     1    z  o   o 

to  (5.6.1). 
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Case  2:   z1  =  z*  >  z' .   Then  with  a  development  similar  to  that 
o     1    2 

in  Section  2.7,  the  points  (x' ,y' )  on  the  line  segment  determined  by 

o 

(x*  ,y')  and  (x' ,y' )  are  all  optimal  solutions  to  (5.6.1)  where 
o   1        o   2 

yi  =  max[Di-p"s-  lb'  - (z;  -g(p'»/w(p')}] 

y'2=   -n  [d- ,  p^,  {b'  ♦  («;  -G<P'»  TO')}]. 

Case  3:   z'  =  zl  >  z' .   Then  determine 
o    2     1 

x^  =  max  [r  ,  p™,  {a1  -  (z^  -  G(P'))  W(P')|  J 


x;,  =  mi 


Then  the  points  (x' ,y' )  on  the  line  segment  determined  by  (x' ,y') 

o  1   o 

and  (x' ,y')  are  all  optimal  solutions  to  problem  (5.6.1). 
2   o 

By  Property  5.6.1,  the  optimal  solution  to  the  problem  (5.2.1) 

is  obtained  by  taking  the  inverse  rotation,  T   ,  of  the  solutions  to 

(5.6.1).   In  Case  1,  (x  ,y  )  =  T   (x' ,y' )  is  the  unique  solution  to 

o   o         o   o 

problem  (5.2.1).   For  Case  2,  the  points  (x  ,y  )  on  the  line  segment 

o   o 

determined  by  (x  y  )  =  T   (x1 ,y')  and  (x  y  )  =  T   (x' ,y' )  are  all 
11  ol  22  o2 

optimal  solutions  to  problem  (5.2.1).   Similarly  for  Case  3,  the 
optimal  solutions  are  on  the  line  segment  determined  by  the  points 
O^.yj)  =  T_1(x^.y^)  and  (x^y^  =  T_1(x^,yM. 

In  the  above  solution  procedure  for  problem  (5.2.1),  there 
are   several   optimization  problems  to  minimize  or  maximize  some  func- 
tion over  S  or  S  x  S,  for  example,  determining  b   .   Thus,  obtaining 
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the  solution  to  problem  (5.2.1)  depends  on  being  able  to  solve  some 
optimization  problems  as  subproblems.   Assumptions  have  been  made  so 
that  an  optimal  point  to  these  functions  exists; it  is  recognized  that 
obtaining  these  optimal  solutions  may  be  difficult.   The  level  of  dif- 
ficulty depends,  of  course,  on  the  properties  of  the  functions  involved, 
that  is,  w(P) ,  g(P),  and  on  the  properties  of  S.   Chapter  2  shows  that 
for  the  special  case  of  S  finite,  these  optimization  problems  are 
finite  maximization  or  minimization  problems  and  are  readily  solved. 


CHAPTER  6 
SUMMARY  AND  FUTURE  RESEARCH 

6 .  1   Summary 

A  study  of  several  minimax  location  problems  is  presented  in 
this  dissertation.   Each  problem  is  to  locate  one  or  several  new  facil- 
ities in  the  plane  with  respect  to  existing  facilities  so  as  to  min- 
imize the  maximum  cost,  where  costs  are  specified  functions  of  the 
rectilinear  distance  between  new  and  existing  facilities  or  among  new 
facilities.   For  all  but  one  of  the  problems  considered,  constraints 
are  assumed  that  are  upper  bounds  on  the  rectilinear  distances  between 
new  and  existing  facilities. 

For  each  problem  considered,  a  formulation  is  given,  possible 
applications  are  discussed,  related  literature  is  surveyed,  properties 
of  the  problem  are  developed,  and  a  solution  procedure  is  provided. 

For  the  problems  considered  in  Chapters  2,  3  and  5,  the  develop- 
ment shows  that  a  transformation,  T,  introduced  in  a  study  of  location 
problems  by  Francis  [15]  ,  applies  to  the  respective  problem  and  allows 
a  decomposition  into  one-dimensional  subproblems.   The  optimal  solutions 
to  the  subproblems  determine  the  optimal  solution  to  the  original  problem. 

Solution  procedures  are  established  for  each  subproblem  that 
take  advantage  of  the  properties  of  rectilinear  distance,  the  objective 
function,  and  other  structures  of  the  subproblem.   In  Chapters  2  and  5 
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a  solution  procedure  is  given  for  the  one  facility  subproblem  that 
expresses  the  unique  optimal  solution  in  terms  of  specified  parameters 
of  the  subproblem;  the  parameters  are  specified  by  an  expression  that 
is  determined  for  the  optimal  objective  function  value  of  each  sub- 
problem.   In  Chapter  3,  the  multif acility  subproblem  is  written  as  an 
equivalent  weighted  network  flow  problem  which  is  solved  by  deter- 
mining flows  of  optimal  "cost-to-weight"  ratio  in  the  network. 

For  the  p-center  problem,  considered  in  Chapter  4,  properties 
of  the  rectilinear  distance  enable  the  specification  of  a  finite  number 
of  candidates  for  the  optimal  objective  function  values  for  all  p-center 
solutions.   An  iterative  solution  procedure  is  provided  that  deter- 
mines the  optimal  p-centers  for  all  values  of  p.   The  transformation  T 
is  used  in  the  computational  development  of  the  iterative  procedure. 

In  Chapter  5,  properties  are  presented  of  a  one  facility  minimax 
location  problem  with  respect  to  a  compact  set  of  existing  facilities 
and  using  general  distances.   These  properties  extend  similar  results 
previously  reported,  [13],  for  a  finite  set  of  existing  facilities. 

6. 2   Future  Research 

One  possible  direction  for  future  work  might  be  toward  includ- 
ing more  constraints  on  the  problems  considered.   The  solution  proce- 
dures developed  for  the  one  facility  and  multif acility  problems  in 
Chapters  2  and  3  suggest   that  linear  constraints  on  the  location  of 
the  new  facilities  could  be  handled.   Assuming  linear  constraints,  the 
linear  programming  formulation  of  the  transformed  problem  does  not 
decompose  into  independent  subproblems,  but  the  constraints  have  block 
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diagonal  structure  with  coupling  constraints  corresponding  to  the 
transformed  linear  constraints.   The  subproblems  in  linear  programming 
decomposition  appear  to  share  many  of  the  properties  of  the  subproblems 
obtained  in  Chapters  2  and  3.   An  investigation  into  the  subproblems 
might  lead  to  an  efficient  decomposition  procedure  for  the  linearly 
constrained  facility  location  problem. 

As  mentioned  in  Chapter  4,  the  solution  procedure  for  the 
p-center  problem  could  be  made  more  efficient  if  good  bounds  were  known 
for  the  optimal  objective  function  value  of  the  p-center  solution  for 
all  values  of  p.   It  would  appear  that  the  properties  of  the  rectilinear 
distance  would  be  useful  in  any  investigation  into  these  bounds. 

Another  area  of  research  mentioned  in  Chapter  4  is  an  investi- 
gation into  the  p-center  problem  in  one  dimension.   An  efficient  solu- 
tion procedure,  capitalizing  on  the  structure  of  the  one-dimensional 
problems,  that  improves  on  or  bypasses  the  set  covering  problems,  might 
lead  to  a  more  efficient  algorithm  for  the  p-center  problem. 
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